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ABS ERACT 


A assessment OF Some approximations in the 
dynamical theory of electron diffraction has been 
eae out. The approximations which have been 
considered are firstly, the s-value approximation 
fOr calculating the diagonal elements or the A matrix, 
secondly, the Kinematical approximation, and thirdly, 
the approximations which can be used to take the 
imaginary part of the lattice potential into account 
when absorption is considered. 

The s-value approximation has been found to 
Dena eleausOlagle approximation uncer the ditiraction 
conditions which were considered. 

The kinematical approximation has been assessed 
for large déviations of the first order reflection g | 
from the Bragg condition by comparing the kinematical 
theory to the many-beam theory. It has been found that 
the kinematical theory is not as good an approximation 
aS comparison with the two-beam theory had indicated. 
In addition, it has been found that the kinematical 
LHe OLY vals) as POO is approximation when the deviation para- 


meter, s of a non-systematic reflection is equal to 


ee 
Zero OL ane It has also been. found that the kinematical 
theory becomes a poorer approximation as atomic number 


or accelerating voltage is incréased or if a lower order 


systematic row is considered. 
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In order to assess the approximations for taking 
absorption into account, a general formulation of the 
dynamical theory of electron diffraction including 
absorption has been developed. This theory is based 
on generalized perturbation theory. It has been shown 
that the introduction of the imaginary potential leads 
to a mixing of the unperturbed Bloch waves characteristic 
of the system without absorption. The assumptions 
involving the importance of Bloch wave mixing which 
are required to reduce the exact theory to the low 
order perturbation theories previously published in the 
literature have been examined. This is of interest since 
an understanding of these assumptions provides a physical 
basis for comparing these theories as well as giving the 
insight required to develop new approximations. The 
details of a new approximate theory, referred to as the 
sub-matrix method, have also been given. Numerical cal- 
culations have been presented which show that it is 
important to take the mixing of Bloch waves into account 
in high atomic number materials or when quasi-degeneracy 
or degeneracy of the kinetic energies of the unperturbed 
‘Bloch waves is obtained. A comparison of the approximate 
and exact methods for taking absorption into account has 
been carried out and it has been found that the sub- 
Matrix approach has a number of advantages over other 


approximate methods. Firstly, it is applicable under 
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conditions which lead to either non-degenerate, quasi- 
degenerate or degenerate Bloch waves. As well, the 
method can give very close agreement with the exact 
theory while at the same time providing computational 
advantages. 

As an illustration of a phenomenon where Bloch 
wave mixing is important, the critical voltage effect 
has been examined. It has been shown that it is 
necessary to take Bloch wave mixing into account in 
order to obtain a PhYysicalsunderstanding of the con- 
trast mechanism involved. For small deviations from 
the exact second order Bragg condition, this mechanism 
has been found to be different from the standard theory 


mechanism. 
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CHAPTER 1 


THE HISTORY OF ELECTRON DIFFRACTION THEORY AND ~ 


INERODUCTION TO THE’ PRESENT “STUDY 


deo b Introduction 


The electron microscope is an important tool 
in the study of the structure of materials. In order 
for the microscopist to interpret the experimental 
micrographs ag tee ey from his instrument, it is 
necessary that reference be made to some theory which 
relates image contrast to specimen structure. In 
some cases, such aS in the interpretation of micro- 
graphs of Piological specimens, useful information 
can be obtained by using the mass thickness method 
for interpreting image contrast, This method -.assumes 
that the intensity reaching a particular aréa of the 
Micrograph is proportional to the product of fife dude 
times the thickness in the corresponding area of the 
specimen. The interpretation of electron micrographs 
OL cheery ite materials, however, requires the use of 
a more complicated theory which takes the diffraction 
= electrons into account. In this case, useful 
information can usually only be obtained by comparing 
experimental micrographs to the results of theoretical 


calculations based on electron diffraction theory. 
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In using electron diffraction theory in practice, 
approximations of various kinds must be made. The 
validity of some of these approximations will be 


examined in this thesis. 


1) “elite torical Gintroducliongcosc nen, neory sores ectron 


Diffraction 


The phenomenon of electron diffraction was 
first discovered experimentally by Davisson and Germer 
(1927) and almost simultaneously by Thomson and Reid 
(1927). These authors confirmed the de Broglie postulate 
of wave particle duality in the case of electrons 
(de Broglie, 1924) and provided the basis for the 
development of the electron microscope. In the year 
following these first experimental discoveries of 
electron diffraction, Bethe (1928) published a quantum 
mechanical theory describing the interaction of a high 
energy electron with a periodic crystal potential. 
This theory, known as Bethe's formulation of the dynami- 
cal theory of electron diffraction, now plays a very 
important role in the analysis of electron diffraction 
effects. 

Initially, however, after the development of 
the electron microscope by Knoll and Ruska (1932), 
the interpretation of the images from crystalline 


specimens was carried out in terms of the simpler 
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kinematical theory of Pastor diffraction. * This 
theory is similar to the kinematical theory of X-ray 
diffraction which had been used successfully to 
explain many aspects of the interaction of X-rays 

and crystals. However, in the case of electrons, 

the kinematical assumptions of single scattering and 
of negligible diffracted beam intensity in comparison 
to that of the directly transmitted beam are not 
usually justified. It is not surprising, therefore, 
that a number of authors (for example, Blackman, 1939; 
von Borries and Ruska, 1940; Heidenreich, 1942; 

Hillier and Baker, 1942; Boersch, 1942, 1943; Kinder, 
1943; Heidenreich and Sturkey, 1945) reported experi- 
mental results which could only be partially explained 
in terms of the kinematical theory. ASa result, there 
was a renewed interest in Bethe's dynamical theory. 

In Bethe's theory any number of strong diffrac- 
ted beams could, in principle, be taken into account. 
In practice, however, it was found that analytical 
expressions for the intensity of the different difiirac= 
ted beams could be obtained easily only for the case 
where two reflections, the directly transmitted beam 
and one diffracted beam, were taken into account. 

These analytical solutions gave good agreement with 
experiment in strong beam imaging situations where only 


one low order reflection was close to its Bragg 
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condition. As a result, the use BE this two-beam 
approximation of the dynamical theory became wide- 
spread. 
In order to take the effects of additional 
weak beams into account and yet retain the simplicity 
of the two beam analytical solutions, it is possible 
to use Bethe's second RA OPOC INA On (see for example, 
Bethe, 1928; MacGillavry, 1940; Kato, 1952; Pinsker, 
1953). In this approximation, the effects of additional 
weak beams are included by adding correction terms to 
che lattice potentials corresponding to the two strong 
beams considered. As would be expected, however, | 
Bethe's second approximation was found to be inadequate 
when three or more strong beams were present, and 
interest pve athee turned to theories which could tks 
the dynamical interactions of many beams into account. 
There were two principal approaches which were 
taken in the late 1950's for calculating intensities 
in a many-beam situation. The first was quantum 
mechanical in nature while the second involved the 
use of optical models to obtain expressions for the 
diffracted beam intensity. Although Fujiwara (1959) 
showed that the Schrodinger wave equation could be 
solved by the method of higher order Born approximation, 


the quantum mechanical approach which was generally 
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adopted was Bethe's approach. In order to make Bethe's 
Eon redeli more suitable for many-beam numerical calcu- 
lations, af eigenvalue formulation of this theory was 
developed (see Niehrs, 1959; Fujimoto, 1959; Howie and 
Whelan, 1960). The eigenvalue formulation of Bethe's 
theory has gained wide acceptance for calculating many- 
beam effects and is the method used in this thesis. 
Methods based on obtaining the eigenvalues and eigen- 
vectors by approximate methods have also peal published 
byhGevers:(1970) andeSerneels “and: Gevers (1973A) .'“-In 
these methods, the diffracted beams are classified into 
strong and weak beams. The eigenvalue equation is then 
solved using partition theory. In applying this theory, 
it is assumed that the strong beams interact dynamically 
and the weak beams kinematically. The latter method, 
termed a kinematical approach by Gevers, should not be 
confused with the kinematical theory discussed in this 
‘thesis. 

The second of the two principal approaches 
which were developed for taking the dynamical inter- 
action of many strong beams into account involves the 
use of optical models to obtain expressions for the. 
diffracted beam intensity. Examples of this approach 
are the cil geantceus (1957) theory, which was deve- 


loped from scattering principles in light optics, and. 
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the Howie-Whelan (1961) differential equation Aa 
which was based on the plane wave reflection approach 
of Darwin (1914). These methods have been justified 

by showing that they are equivalent to the more basic 
quantum mechanical approach (Fujiwara, 1959; Hirsch et 
al.,1965; Goodman and Moodie, 1974). The relationships 
between the different many-beam dynamical theories have 
been summarized by Goodman and Moodie (1974). 

It can be seen from this brief historical sketch 
ofathe,development of electron diffraction theory that 
the kinematical, the two-beam, and the many-beam 
theories came into use as theories were required to 
better describe the interaction of the fast electron 
with the crystal. Although the many-beam theories 
have proven to be very useful for calculating diffracted 
beam intensities, these theories are numerical in form 
and results are generally forthcoming only with the aid 
of a computer. For this reason, the two-beam dynamical 
and kinematical theories have not fallen into disuse. 
Rather, they have remained as important qualitative 
a1dssrOr Givang ansight into avitraction. phenomena, 
These theories are also used in a more quantitative 
sense under diffraction conditions where they are 
thought to be good approximations. The two-beam and 
kinematical theories are therefore discussed at length 


in the textbooks of electron diffraction (see for 
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ex ample Pinsker, 1953; Heidenreich, 1964; Hirsch et 
al., 1965; Amelinckx et al., 1970) and have been used 
extensively in the literature. Recently, for example, 
the kinematical, and the two-beam and many-beam 
dynamical theories have been used by Cockayne (1972) 
in the analysis of his weak-beam results. Under weak- 
beam diffraction conditions, the intensity of the 
diffracted beam being imaged is small and the appro- 
Ximations of the kinematic theory are more WEY! to 


apply. 


1.3. Introduction to the Present Study 


ineaneroducing the research to a carried Out 
in this thesis, it is first convenient to discuss some 
of the basic concepts involved in Bethe's approach to 
the dynamical theory of electron diffraction. (The 
details of Bethe's approach will be considered in 
Chapters inary 3.) Bethe's approach is a quantum 
mechanical approach in which the fnteraction of the 
fast electrons with the crystal is described by the 
Schrodinger wave equation. The wave function solutions 


to this equation are Bloch waves of the form 
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Bloch wave vector and ee is the amplitude of the 

| : : : =e Or) eae 2 : 
Blocn wave in the direction k “i+g. .in- solving the 
Schrodinger equation, the lattice potential is also 
expressed as a Fourier series. If only the elastic 


Scattering Of Electrons is’ considered, then this 


potential is real and can be written 
v(x) =) Vy exp (2Tig-f) (dl 22) 
g 


where Mie is the Fourier coefficient of the lattice 
potential. However, inelastic scattering processes 
also occur. Some of these elecerons will be scattered 
outside the objective aperture and therefore do not 
contribute to the image. These Serene are effec- 
tively absorbed and it is in this context that the 
term absorption is used in this thesis. As will be 
shown in Section 2.2, absorption can be taken ae 
account in the dynamical theory by the introduction 

of an imaginary part to the lattice potential. The 
lattice potential is then considered to be complex and 


can be written 


v(t) + iv' (4) =) (v 


+ iv!) exp (27ig-r) é ips (alice) 
g g 


g 


The expression for the Bloch wave amplitude in 
equation 1.1 involves the sum over an infinite number 


4 > . 3 
of reciprocal lattice vectors g. However, when this. 
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Bloch wave function is incorporated into Bethe's 
theory, an infinite set .of equations results (see 
Gnapter. 2). in order. to obtain, numerical soltttions, 
this infinite set of equations must be reduced to a 
finite set by taking only a finite number of reciprocal 
lattice SiSceee into account. Since every point in the 
reciprocal lattice corresponds to rossi ore seraeayeta 
from the real lattice, this means that the sere tete OL 
diffracted beams taken into account must be finite. 

By reducing the number of beams considered, the many- 
beam and two-beam dynamical theories can be developed. 
If only the directly transmitted beam and one diffrac- 
ted beam of small intensity are considered, then the 
kinematical theory. can be obtained. These theories: 
will be considered in the first part of this thesis 
(seenocection. 1.3.1). 

Another important aspect of the dynamical 
theory is the manner in which the imaginary Fourier 
coefficient iV, of equation 1.3 is incorporated into 
the theory. The method widely adopted inthe literature 
is based on the fact that th is generally much smaller 
than Vg: This permits the imaginary potential to be 
taken into account by using first order perturbation 
theory for the non-degenerate state. (This method is 
referred to as the standard method in the literature.) 


In the second part of this thesis,the validity of this 
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method is examined and an alternate method is developed 


(see Section 1.3.2). 


1.3.1 An Examination of the Approximations Involved in 


the Many~Beam, Two-Beam and Kinematical Theories 


There are a number of approximations which can 
be used in taking the different diffracted beams into 
account. The “£irst is..the many-beam theory, which takes 
a Shae number N of the beams into account. The actual 
‘number of beams included is generally determined by 
USinG es Criterion co order the different diffracted 
beams according EORCNeLEIeTLSCe Onitiic Intensity beams 
are then creed tuvorder of their decreasing importance 
until the inclusion of additional beams does Merce e 
in any significant change in the diffracted beam inten- 
Sity. In the many-beam formulation of the dynamical 
theory, ten to one hundred beams are generally required 
and the theory is usually numerical in form. Because — 
all the important diffracted beams are included in a 
properly executed many-beam calculation, the approxima- 
lens Tnvolved in developing the theory are very small. 
Recotdingiys the many-beam theory has been found to be 
in good agreement with experiment under all diffrac- 
tion conditions and has become the Secebli shed method 
for interpreting experimental pean (see for example, 


Cann and Sheinin, 1974A, 1974B; Sheinin and Botros, 
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1971; Ayroles, Mazel and Ajustron, 1973). 

ve The second approximation obtained by taking 
only a finite number of diffracted beams into eaters 
is the two~-beam theory. The theory considers only the 
directly transmitted beam, o, and one diffracted pear 
g. This approximation results in expressions for the 
diffracted beam intensity which are analytical in form. 
The two-beam theory is been a questionable approxima- 
tion, however, since higher order reflections of ene 
type ... -2g, -g, 2g, 3g, ... are always excited when - 
the reflection g is excited. (These reflections are 
referred to as systematic reflections.) The effects 
OL these Systematic reflections have been studied by a 
number of atnuors:| For example,.Sheinin (1967). found : 
in examining the (110) dark field intensity in 
molybdenum that the two-beam theory was a good appro- 
ximation for crystal orientations approximately between 
OU AG eo ea CLS thesis, the deviation 
of the systematic reflections from their Bragg con- . 
ditions will be indicated by a deviation parameter Oa 
See Fig. 1). For values of 6 approximately between 2.4 
and 3.2, Sheinin found that the two-beam theory was a 
very poor Ease tine ort For Values. of) 0 73.2 .5Gne 
agreement between the two-beam theory and experimental 


results improved. 
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Fig. 1. The method used in this thesis to indicate the 
deviationwor the systematic reflections: ,£rom 
their Bragg conditions. » The deviation’ of: these 
reflections from their Bragg conditions is 
indicated by a deviation parameter 6 which is 
obtained by dividing the distance OA by the 
distance OG, For example, in the above figure 
GO =64.5. "if 6 = n where mn iS a positive ox, 
negative integer or zero, then the reflection 
RGuLS ihwucs, exact Bragg condition. —ihespomre 
L is the perpendicular projection of the tie-— 
point T onto the systematic row. 
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The third theory which will be considered is 
the kinematical theory. In this theory no further 
assumptions about the number of important diffracted 
beams are made. Rather, it is assumed that the inten- 
sity Remeron e diffracted beam taken into account is | 
small compared to the directly transmitted beam. As 
will be shown in Chapter 3, the kinematical theory is 
analytical in form and can be derived from the two- 
‘beam theory if certain approximations are made about 
the size of the Fourier coefficient Vg: Because the 
intensity of the diffracted beam is assumed to yh small, 
it has been found eepente for very thin crystals) that 
the kinematical theory is a poor approximation for the 
Petar eTaioen g close to its een condition, i.e. for. 

5 Cel s0eeesOowe Vel, Lie wwinemaci cal theoryohas been 
assumed to be a good approximation, even for thick 
dated if the veflection g ¥s~ far from 1ts Bragg 
condition (see for example, Whelan, 1970; Gevers, 1970) 
because he intensity of the diffracted beam will then 
be small compared with the directly transmitted beam. 
This assumption has been justified by the fact that the 
two theories become equivalent for large deviations of 
the reflection g from-its Bragg condition, i.e. for 

6 >> 1.0. It has therefore become widely accepted in 
theerliteraturesthatiitthe kinematical ‘theory iswvalivdrfor 


values of the deviation parameter 6 >> 1,0. 
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Although there is good agreement between the 
kinematical and@eworbeam theories for 06 >> 1.0,'i1t is 
important to note that the kinematical theory still 
takes only two beams into account. However, it is 
possible that other beams will also be important and 
should be included. In eaer to properly assess the 
kinematical theory for 6 >> 1.0, it is therefore not 
sufficient to compare the Minenaeical and two-beam 
theories as has been the practice in the past. Rather, 
it is necessary to compare the kinematical theory to 
the more Poprece mMany-beam theory. Such a comparison 
has not yet been presented enehe literature. A com- 
parison of the eneneties ee cwaebean and many-beam 
theories sce be large values of 6 will therefore be carried 


out in the first part of this thesis. 


J25.2- An Examination. of the Approximations Which Can 
be Used to Take the Imaginary Part of the Lattice 


Potential into Account 


There are two approaches which have been adopted 
in the literature to the solution of the Schrodinger 
wave equation including the imaginary potential iv'(r). 
The first assumes that the imaginary part is much smaller 
than the real part. Low order perturbation theory is 


then used to calculate corrections to the Bloch wave 
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functions found using the real potential. The standard 
method has been to use first order perturbation theory 
for the non-degenerate state (Hirsch et al., 1965). 
The use of this theory has been brought into question, 
however, when the imaginary potential becomes large, 
as in high atomic number materials (Sheinin and Andrew, 
1974), or when the assumption of non-degeneracy does 
not apply. Under certain diffraction conditions this 
second assumption is questionable, as for example, 
near critical voltage diffraction conditions (Andrew 
and Sheinin, 1974A) or when a non-systematic reflection 
is excited (Cann and Sheinin, 1974A, 1974B). Although 
higher order perturbation theories have been developed 
(Serneels and Gevers, 1972, 1973B) and the theory has 
been modified to properly take into account Bloch waves 
which are precisely degenerate (Sheinin and Cente T9755 
Serneels and Gevers, 1973B), this perturbation approach 
is incomplete in that no method is available to handle 
quasi-degenerate Bloch waves. (In this thesis the 
term Bloch wave degeneracy will be taken to mean an 
equality in the kinetic energies of the Bloch waves, 
while quasi-degeneracy will be taken to mean nearly 
degenerate Bloch waves or close branches of the dis- 
persion surface.) 

The second and alternative approach is to solve 


the Schrodinger wave equation including the complex 
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potential without making any approximations about the 
size of the perturbation or about how close the 
kinetic energies of the Bloch waves are.to one another 
(Hewat reported by Lally et al., 1972; Thomas, 1972; 
- Andrew and Sheinin, 1974A). eecer this method, 
referred to as the exact method in this thesis, does 
not give any insight into the interaction between 
Bloch waves which results from the introduction of the 
complex lattice potential. Furthermore, the relation- 
ships between the exact method and the perturbation 
approaches are not ee ente | 
Because of the inadequacies of the two avail- 
able approaches,: a study is arcercncen Le) one second 
part of this thesis of the approximations which can be 
ween to take the imaginary part of the periodic lattice. 
potential into account. In this research, generalized 
perturbation theory is used to develop a general for- 
mulation of the Seeeetieeer theory of electron diffrac- 
tion, including absorption. This new-.formulation of 
Sore theory is very useful since it provides a physical 
basis for comparing the diffraction conditions under 
which the various low order perturbation approximations 
are applicable as well as giving the insight required 
to develop new sere w ime niet ee The details of a new 
approximation, referred to as the sub-matrix method, 


will be presented. This approach has.a number of 
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advantages over other approximate methods. Firstly, 

it is applicable under conditions which lead to either 
non-degenerate, ques Gdedenebats or degenerate Bloch 
waves. As well, the method can give very close agree- 
ment with the exact theory while at the same time 
providing advantages from Seecerietin Meine point of view. 

In order to List rate the usefulness of the 

new theory, it is then applied to the situations where 
the standard theory does not apply, namely for high 
atomic number materials and for diffraction conditions 
which lead to degenerate or quasi-degenerate Bloch 
waves. Because of the ability of the new theory to 
handle quasSi-degenerate and degenerate states, this 
theory is also used in a study of the mechanism leading 
to the reduction of the second order dark field inten- 
sity which is associated with the critical voltage 


effect. 
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CHAPTER 2 


ASPECTS OF THE DYNAMICAL THEORY OF 


ELECTRON DIFFRACTION 


2.1 Basic Outline of the Dynamical Theory 


2.1.1 The Solution of the Schrodinger Wave Equation 


The dynamical theory of ocr rsnucireract ion: 
proposed by Bethe (1928), starts with the Schrodinger 
wave equation for an electron in a crystal potential 
| V(r). This equation can be written as 


81-m e 


weve) ee (E+ V(r)) v(t) =O. |. (oe) 
where (rt) is the wave function of the electron, and E 
is the Sten tat through which the electron was accele- 
rer harere entering the crystal. The quantities Mo 
and e denote the rest mass of the electron and the 
electronic charge respectively, and h is Planck's 
constant. Since the crystal potential, V(r), is periodic 


in nature, it can be expressed as a Fourier series of 


the form 
> > > h2 > > 
Vilrjo =e) ov exp (27ig-r) = 5=—— ) U_ exp (2nig-r) (252) 
x ame og 


where the summation is over all reciprocal lattice vector 


h 


g and Veg = (n*/(2me) )U, is the g* Fourier coefficient of 
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the lattice potential. If only elastic scattering is 
considered, the crystal potential is real, V(r) = Vv*(r), 


and therefore 
U == lic 2 a es) 


in addition, if the crystal is centrosymmetric, which 
is the case for the crystals considered in this thesis, 
then an origin can be chosen such that vit) Vier): 

In this case, it can be seen from equations BR einel oe 


that 


U. = U =U 7 (olage 


and the Uays are allireal. 

Hing the case of electrons “with energies of the 
order of those used in the electron microscope, E > V(r) 
and the Schrodinger equation can be solved using the 


nearly free electron approximation. The solutions out- 


side the crystal, where V(r) = 0, are plane waves of the 


form 


W(t) = exp(27ix-r) - (2.5) 
where the magnitude of the wave vector X is given by 


5 
x= (AS S226) 
h | 


In the crystal, the effect of the lattice potential, 


V(r), results in solutions to the Schrodinger equation 
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which are Bloch waves, that is plane waves modulated 
byra unctdon which has the periodicity of the lattice. 


Accordingly, solutions to equationw2...are ofethe, form 


y(t) = [bp eo) r)> ape oh tl git ) veegenaqade +g) -r} 
g 
Cy) 
where g (1) is the Bloch wave vector and ¢ () va )y is 


the amplitude of the Bloch wave in the direction 
pee If equations 2.2 ahde ni GeesUubS Ll tuted iil co 
thes sch rodinger equations 2. UL, the fol Lowinggexpression 
is obtained: 

ee (2) &) 74x? 40, Ios?) (1) )4 ca cf) BO) y 
g Toa Be 


ee aye) = 0) e588) 


where the prime on the second summation indicates that 
the term h = 0 is excluded. Since each of the terms 


Ss it 


involving exp{2mi (kK g)-r} is linearly independent, 
the coefficients of these terms must all be identically 
Cqudleto 2eLlo. Swlilis Condition Lresults ina Set or 
equations of the form 

ee) ) yt) ) ey ig ela 
Keim (kK ) me k ae (k ) =O ( 2929)) 
g qptin! ee) h'g-h 
where 


Kee yeu | C255 18.08) 
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and 

Rey = RT yt S : | | (2211) 
The quantity K is the magnitude of the Siesrren ive 
vector in the crystal after correction for the wave- 
length change due to the mean crystal potential Us: 

The homogeneous set of equations 2.9 gives the 
general relationships between the erritddes ei cee 
the Fourier coefficients Ugi and the Bloch wave vectors 


oe This equation will have solutions only if the 


determinant formed by the coefficients vanishes, i.e. 


y? U elem s 
uf ere (aed Wee ab eaeyr (era): 


It can be seen from equations 2.10 and 2.6 that this 


condition relates the Bloch wave vector Rae 


to the 
energy eE. The locus of the end points of the wave 


vectors k Satiastying this) relacrvon traces out branches 


aay) 
Ee 
> 

of a surface of constant energy eE in k space which is 

known as the dispersion surface. For this reason, 

equation 2.9 is known as Bethe's dispersion equation. 
The method which will be adopted in this thesis 

to solve Bethe's equation 2.9 is based on the eigenvalue 


method. The manner in which the eigenvalue equation 


is derived can be understood by considering inks an PG 
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Sphere of radius 
K centered at g 


Sphere of radius 
K centered at O 


Brillouin zone Ewald sphere of radius 
boundary - -K centered at T 


Fig.2. The dispersion surface and the Ewald sphere 
construction for high energy electrons. 
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‘In this figure the incident wave vector, after allowing 
for the ertect ot refraction; 1s) given by K. For 
Simplicity only the two epee ar cee lattice points, o 
and g, and the two dispersion surfaces, (1) and (2), 
have been included. In Bee to determine the allowed 


Bloch wave vectors er 


within Ehescrvystal, it; as 
necessary. to apply the boundary conditions at the top 
SuUrLLaCemOre tic crystal, s:lnsorderstor the wave functions 
to join smoothly at the boundary, it is necessary that 
the wave function and its derivative normal to the 
Peres should be continuous at the surface. The fact 
that these conditions must hold ye any point on the 
surface requires that the wave vectors kia, have the 
Same tangential components along the surface as Ke 
Generally, reflected waves from the surface can be 
ignored up to angles of incidence of approximately 80° 
(see Hirsch et al., 1965) because of the large difference 
between the energy of the incident electron and the 

_ lattice potential. 

In order to apply the above boundary conditions, 
it is necessary to consider the inclination of the 
crystal surface with respect to the incident beam. 
Limorder £66 simplify this process, it will be assumed 
in this thesis that the symmetrical Laue case holds. 

In the symmetrical Laue case, the different Laue zones 


(see Fig. 3) are assumed to be. parallel to the crystal 
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Fig. 3. Ewald sphere construction showing reflections 
in different Laue zones. Because |K| is much 
greater than the reciprocal lattice spacing, 
the Laue zones are nearly perpendicular to the 
incident beam. 
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surptace.’ }Forkthe small Bragg angles found in high 
energy electron diffraction, this means that the 

crystal surface is assumed to be nearly perpendicular 

to the incident beam. . It has been shown by Whelan and. 
Hirsch (1957) that this symmetrical Laue seo eeimabTon 
is a good approximation in the two-beam case for angles 
of incidence up to approximately 70° to 80°. Spencer 
and Humphreys (1971) have noted that in nearly all cases, 
the symmetrical Laue approximation is a Hood iene, but 
that there ate cases where it is not, as for example in 
the description of scanning electron microscope channel- 
ing patterns from a tilted specimen. 

The boundary conditions can now be applied in 
bil aan by drawing, through point T, a line perpendicular 
to the crystal surface. The points where this line ) 
intersects the branches of the dispersion surface are 
termed wave points and determine the allowed Bloch wave 
vectors within the bevatal, In the symmetrical Laue 
Case forythe reflection g, this line will be parallel 
to the Brillouin zone boundary as shown. Now let y (4) 
be the di weaios of a possible wave point from T and let 
ce be the distance of the reciprocal lattice point g 
from the Ewald sphere, measured in the Zz direction 
(i) 


(Normal tO.g). MoOtestawl, cy and Sq! it can be seen 


from Fig, 2 that 
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ra) O71") 2 (i) 2 
K. - (k ~ wee (Ke + 
( g ) ( Sy cos Pe Y cos 8) 
= K - (K+ =5 Cc ore 
( (y g? os 8,) i 13) 
seen Be ay ao 3 
== od Ss ; ° 
Y Sg co | aK cos By 
ey cee is assumed thate~Kn>>ag, and Kys> (yO mg) then 
COS io. os aL 26 
Pees (2.14) 
(y Bac 
and aK Sal) eae (2315, 
Equation 2.13 then reduces to 
ss aes Ue te Seok ty 84) (2-16) 


and Bethe's equation 2.9 can be written in the eigen- 


value form 
ac‘) mi y fg) : (oie 


Hene ¢ (4) is a column vector whose elements cate abeachie 


components of the Bloch wave of equation 2.7 and A is a 


matrix with elements ee =O, A s_, and epee eh OCS 


Hopes 


where g #h.- The eigenvalue yt) is related to the z 
component of the Bloch wave vector ae through. the 
relation 
Le ee (i) | 7 
(ko es me ao 5 | (2). 18) 


ss | * 


6 


bn 5» << AH seers bemyess 


(#4 
i — \47 r) 
es ee 
7 fia 
c subsi dst tL£.$ noltsaspa 
c) SC (b)« ,. 
lo Cj ( - 4 J y) AS- c { ») im” y ; 
: - . A 
7 = 
~nepis oft at motsiaw ad aso €.%) noijsupe eB” ented Sas 


mtet sulev 


7 2 
; 7 


eee: (ton 


(i) | at ey | Sir 
Sid Sis se etinemels seonw rotoSay apuloo 6 Bf 2 soe 


7 


‘os ei A bun ¥.S celLasups- to svaw doolf,. 34 to : etmenogmos 


| 7 
wee PhS y= 2A - “aa ednemals idkw wide i 
; SAW ran te %, Site 


7: - 
sit? ody buderon pb ‘* ly. 28 Covneeio ot 
: | ; | = sf ik vty he * 7 ari. 
| a —e oy a) 


as) = ” fy = 
“_ as 


* 


| 


In order to be able to refer to the different Bloch 

| waves present, it is convenient to label the Bloch 
waves and their corresponding dispersion surfaces in 
accordance with some accepted convention. The method 
generally adopted in the literature is to use the 
untceger Paberls) i= 727397 e. in order of pees eins 


(i) 


values of (ko 


= (Humphreys and Fisher, 1971). 

The formulation of the dynamical theory just 
presented has been éxpressed in terms of the eigenvalues 
and eigenvectors, of an infinite A matrix. The applica- 
tion of the theory in practice, however, clearly will 
require that this matrix be finite. The approximations 
Caen are involved in reducing the matrix to a finite 
size will be GOR OLeE: in the next chapter. In the 
remainder of fis chau: erjatie methods which can be 
used to calculate the diffracted beam intensity, and to 


take inelastic scattering and relativistic effects into 


account will be considered. 


2.1.2 The Calculation of Diffracted Beam Intensities 


inga PerfectCrystal 


The total wave function p(t) of the high energy 
electron in the crystal can be written as a linear com- 
bination of all the Bloch waves excited, i.e. 


Wit) RE peo pe! Jol, explaniik +)er). . 
er 7 ie CMa Rae (2.19) 
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The coefficients x (4) in this expression are known as 
excitation Pee egenee and determine the extent to 
which the associated Bloch wave is excited in the 
crystal., In order to calculate the intensity I of 

a directly transmitted or diffracted beam g, equation 
ee Lois) Sorted ee components in the directions of 
these beams. If these components are then multiplied 


by the phase term exp (-20ik-r), then the intensity of 


beam g at. a depth -z in the crystal is found to be 


2 


(1) 5) | (2.20) 


2 bs el Gs 
PZ) = 105 °(2 = X Xe exp (Zia 
g (2) =1%, (2) | p 5 p(2miy 
where oa) is the amplitude of the beam g. 

The excitation coefficients required in equation 


2.20 can be calculated from the boundary conditions at 


.the top surface Ore the crystal, <1. 
POM aw Pe (0). 0 6g 40). = (2.21) 
bone tvonr2e20. together with 2.21 gives 


CX =u "425229 


where. C, 1s a matrix having the elements oe in the eee 
row and pet column, X is a column vector containing the 


excucLallon cOGt ii cuen ts x) in the rt row, and Uys 
a column vector containing the diffracted beam amplitudes 
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xh) required in equation 2.20, can then be found by 


obtaining numerical solutions to the nonhomogeneous 
set of linear equations 2.22. The effort required to 
solve this set of equations can be greatly reduced, 
however, by noting that Aid A-matrix of equation ‘2.17 


ToeLealpandesyMMetrics. Hence tie matrix sc or normalized 


eigenvectors is orthogonal and cht “s ats the transpose 


of C. If equation 2.22 is multiplied from the left by 
ee then 
x 2.0) u (23) 


andsitican immediately be seen that the excitation 


coefficients ae se given by the elements in the first 


BOweoheG pj. Le. 


xf) 2 Gl) | (2.24) 
fe) : 


The expression for the diffracted beam intensity there- 


fore becomes 


2 28 (i). (4) 
T(z) = [0 (2) | 12 Shaner 


exp (2niy ‘t)2){7°... (2.25) 


2.2 Effects of Inelastic Scattering 


The high energy electrons incident on the crys- 
tal can undergo both elastic and inelastic scattering. 
The dynamical theory as formulated by Bethe takes only 


elastic scattering into account. It has been found, 
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i however, that the inelastic scattering of electrons 
outside of the objective aperture leads to the very 
important anomalous absorption effect. This section 
will consider the important types of inelastic scatter- 
ing processes and the modification which is usually 
made in Bethe's theory to take these processes into 
account. , This will be followed, by. a. discussion: of a 


anomalous absorption effect. 


2i2e ok wlneLastac Scattering Processes 


The high energy electrons travelling through the 
crystal can be involved enisetene different inelastic 
scattering processes. The three most important of these 
processes are plasmon scattering, short range single 
electron scattering, and phonon. or thermal diffuse 
scattering. 

Plasmon scattering results from the long range 
Coulomb interaction between the high energy electron and 
the conduction electron gas aS a whole. The incident 
electron loses energy by exciting long wavelength AuAarey 
oscillations, called plasma oscillations or plasmons, 
in the electron gas. These oscillations have a charac- 


teristic frequency Ww, and an excitation energy fw, which 


if 


is usually in the range of 3 eV to 30 eV. The mean free 


paths for plasmon excitation lie typically in the range 
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of 1000-5000 A (Pines, 1963), so that in many cases 
almost all the electrons emerging from the len oa 
have lost energy by this fee isa B: Because the. plasmon 
wavelength is long compared to the lattice spacing 
(i.e. the wave vector q << g), it is not expected that 
Bragg reflection would occur during plasmon excitation 
(Howie, 1963, 1970). Only small angle scattering is 
therefore dupautwenrdnd most electrons suffering plasmon 
losses ia T pass through the objective aperture (see 
also Heidenreich, 1964). 

In addition to the long range Coulomb interac- 
FCLton -yust discussed, thecincident electrons ican interact 
be short range with the individual electrons in the 
crystal. Most of the electrons inelastically scattered 
by single electron excitation are thought to scatter we 
through angles less than a Bragg angle (Hirsch et al., 
1965). Although it is possible for some of these 
inelastically scattered electrons to pass through the 
Be Veeries aperture, most of. the inelastically scattered . 
electrons e6htrebecing to the image are thought to have 
been scattered by plasma excitation (Howie, LOW Opis 

The third of the important energy loss processes 
is phonon or thermal diffuse scattering,which involves 
the scattering of Bloch waves by thermal vibrations. | 
‘In the areWenart of Lit ckerertetdig process, the scat- 


tered intensity can be divided into two parts (see for 


ie 


‘ees 4at mox? oAkpyems enovtools efd Lhe prey 
dottesig edz saysesth .bontvom area vd yorene taol even 
prbkonye sotitvsi eft oF bateqmoo piol 2i Atpaelevew 
Jsd3 Bssosqxeo tom ei si ,(R >> P tojosv evaw ond .8,1t) | 
noliasioxs nomaniq oni tue 1HID0 bigow aoktosl tex ppexd 
ai prixsdisoe sipis lilsme vino .(0vel ,coRt ,oiwolt) =f 
fomesiq pnaixzstive anottoslis Jaom bie betusake gtotesreds 
992) Sivjiags svisospdo els Mpuords aesg filw easeol 
(hae, ,dorainsirel oalfs 
optetnt dmolued sprzex pnodt aia oF noirgibts nt 
tosxgini aso enortyels tnebiont sig .boeeuoesb taut nots 
sit aL enortoois [subivibat eft iisw apne Jiona te 
betedisse yileoigesleni enoitosis sid to taoM ,ieseyio Pr 
istjaoa o2 tdpwode sxe cot4stioxs fertosle sfipaka va 
,-45 38 noeitH) sipnas ppsid 5s msnd easel eslpns tepoxds 
sagnt to omos xot sidtesoq at di cpvonalA .(eeek 
ad4 dpyotds 2a6q o3 Anotdoels bexestesd yilesrsesions 
bexetdsv2 yllesidasloni aft to teom ,saddeege eyigsette ~~. 
Sve od tdpuotia sue spemt odd of paksgotliziaaon anonioels 


-(OVEL \siwoll) mottatioxs smasta yd Parse 


S2 


example Hail and Hirsch, 1965) .° The first "term is the 
Bragg scattered term. This term is found to be equal 
to exp(-2M) times the intensity Bragg scattered by a 
static lattice. Here M is the Debye temperature factor 
given by M = Bg*, and B = enue where u2 is the mean 
square displacement of the atoms in a direction perpen- 
dicular to the Bragg planes. The Bragg scattered term 
is taken into account in the dynamical theory by multi- 
plying the Fourier coefficients of the lattice potential 
by the Debye-Waller correction factor exp (-B|g7|/4). 
The second term is called the diffuse scattering term. 
This term is taken into account in the calculation of 
the imaginary part iu! of the lattice potential (see 
Section oe The energy changes involved in 
‘thermal diffuse scattering are small (AE ~x k,T ced O29 i PAS) 
at room temperature) and the scattering has been found 
to be greatest at moderate angles of the order of “10° 
(Howie, 1970). The contribution of thermal diffuse 
scattering is of prime importance to electron absorp- 
tion since most electrons suffering this type of loss 
will not enter the objective aperture and are therefore 
considered to be absorbed. Since thermal diffuse 
scattering is a function of the thermal motion of the 
‘crystal, it is not surprising that it has been found to 


decrease with decreasing temperature (Steeds and Valdré, 


1968; Howie and Valdré, 1964). 
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It should be noted that, for crystals. of the 
thickness commonly examined in the electron microscope, 
nearly all the electrons are expected to have undergone 
inelastic scattering. This point is of interest since 
good agreement with experiment can be found by using 
diffraction theories which assume that none of the 
electrons. reaching the image Nave been inelastically 
scattered. This apparent contradiction has been 
resolved by Howie (1963) who showed that the phase 
relationships among elastically scattered Bloch waves 
are the same as among Bloch waves which have been inelas- 
tically scattered through small angles. fThus, en 
inelastically scattered electrons give the same image 
contrast as those undergoing only Biased scattering. 
This was first shown experimentally by Kamiya and Uyeda 
(1961), who displaced the objective aperture in a manner 
such that only inelastically scattered electrons contri- 
buted to the image. The resulting images were very 
Similar to those obtained with elastically scattered 


electrons. 


2.2.2 Extension of Bethe’'s theory to include Inelastic 
Scattering Processes 
In Bethe's original formulation of the dynamical 


theory, the depletion of the elastically scattered beams 
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by inelastic scattering processes is not taken into 

account. In order to include these inelastic scatter- 
ing effects, Moliére (1939) introduced the phenomeno- 
logical concept of a complex lattice potential. If a 
small imaginary component iv' (r) is added to the real 


potential V(r), then the total crystal potential becomes 


ae 
2me— 


V(r) + iv' (4) = Y (U_+ iU')exp(2mig.r) . OL D49 26) 
g g g 

The addition of this imaginary component leads to 
imaginary components in the Bloch wave vector KR‘) of. 
equation 2.7 and in the eigenvalue y $2) of equation 2.25. 
It will be shown below that this leads to the attenuation 
of the diffracted beams with thickness and, in this 
manner, accounts for the loss of electrons from the 

beam by inélastic: scattering. 

The addition of an imaginary component to the 
lattice potential was substantiated theoretically by 
Yoshioka (1957). Rather than solve the Schrodinger 
equation using a static potential V(r) as Bethe had 
done, Yoshioka used a potential which was able to take 
inelastic scattering effects into account. He showed 
that these additional interactions give rise to an 
additional complex term in the Fourier potential. ‘The 


real part of this additional term is small and is 


usually neglected, while the imaginary part is more 


> 
' jc oe 
otai newted jon af satesnord snesnuss oleae aE 


-ysyie0a Sisesioni seen% sbafont 09 4ob40 AT ey 
~onenonsig sit besuborsni (22 eL) siSrLomM ,otostts pat 
6 22 .latinetod sordtal xelgmoo 5 to 3qsyno5 les ipol 
fsex eft ot bobbs et (1) "VE Pnenogmos yisnipsm. fleme 


; | wv. : 
zemoved leiinetoq lstayr1> fisted silty ners , (2)V Ls tonszog. 


(a8.8) ; (%-fkw 9) qx (Ue +t) Ge (2) "vt + (3) 
oj a@bscl gnsnoqmoo ytenipsmi aidt to noksibbs edit 
to (214 sojosyv ovaw doola ode ak adnenoqmoo yisnipsm? 
28.5 notteups 30 ith, avkeviagie edt at bas €.S mottsupe 
noltisugedds ort of ebsol eindd tends wolsd awode od {itw 32 
eidd ni .bns ekondotds dtiw amesd hetve1tt lb edt to 
eit moti encrsosls to eeol off tzo*t aetnvoose , tennsm 
-pritestsoe coistasleni yd maed 

sd4 ot Anonogmoo yranipsmi os to moistibbs siT 
yd yilsvisorood? batsisassadve esw Inisnetoqg eoldset 
tepniboxrise oft eviog nety asntteA =. (TTCL) ndoisdeo’ 
ped esdjo@ as (Ci)V (nidtsdog pidete s pabes nokdaupe 
saps of eit. asw dokdw Lettiuetoq s been sAoideoY ano’ 


bewoile si Mere ANs es etostie paizossso8 oitasle 


important and represents absorption effects associated 
with inelastic scattering. The contributions of the 
various inelastic scattering processes to the imaginary 
potential have been calculated by several authors 
(Heidenreich,:- 1962; Whelan, 1965A, 1965B; Hall and 
Hirsch, L965. Humphreys and Hirsch, 19607 Radi','1970) ; 
These calculations have indicated that the principal 
contribution to the inelastic scattering of electrons 
outside the objective aperture arises from thermal 
diffuse scattering. (A review of Yoshioka's theory and 
an assessment of the importance of the different inelas- 
tic scattering processes are given by Kambe and Moliére 
BLO 7 O ee) 

There are several approaches which can be taken 
in order to take the imaginary part of the periodic 
lattice potential into account. Since these aeprogetes 
will be discussed in detail in the second part of the 
thesis, only one approach, which has been widely adopted 
in the literature, will be introduced in this section. 
This standard method assumes that the v' (r) << V(r), 
and takes iv' (r) into account by using the methods of 
first order perturbation theory for the non-degenerate 
state. The effect of the perturbation iv' (r) is co 
_ change the energy of the Bloch wave by an amount eAE 
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This. energy change can be related to a change AK, 


in the z-component of the Bloch wave vector, i.e. 


ak) =- = ME : (2.28) 
Ned. 
ie oe 
where aioe is the absorption coefficient of Bloch wave 


Ib), This coefficient can be calculated from the 


integral in equation 2.27 to be 


(1) ee (a) C1) 8 
he soy L , Cy ay oUpy (2.29a) 


In matrix form this equation becomes 
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git) = (ct) Taw eG) 


(22a) 


where A' is a matrix containing the elements Ug-n/2k in 


the got row and ae column. It can be seen from the 
ia) in equation 2.18 that this change 
(i) 
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Accordingly, when the effects of inelastic scattering 
are included in the dynamical theory, the expression 
for the diffracted beam intensity given in equation 2.25 
becomes 
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By examining this relation it can be seen that the 
diffracted beam intensity is effectively attenuated 
or absorbed with depth in the crystal due to the 
presence of the exp (-21q ‘1 z) terms, | 
Calculations of the Bloch wave absorption 
parameters indicate that some Bloch waves are more 
highly absorbed than others. This observation has 
. been explained by several authors (for example, 
Hashimoto, Howie, and Whelan, 1962; Howie, 1966) in 
terms of channeling of Bloch waves. In general, the 
intensity <b {4) |p fh), of a Bloch wave as it propagates 
through the crystal is not evenly distributed, but 
instead, is localized or channeled into certain regions 
of the crystal either at or between the atomic planes. 
On the other hand, scattering processes such as thermal 
diffuse scattering occur with a relatively higher pro= 
bability near the atomic positions thansingzother 
regions of the crystal. Accordingly, Bloch waves 
whose intensity is concentrated at the atomic planes 
will be more highly absorbed and therefore have larger 
absorption coefficients than Bloch waves whose intensity 
is concentrated between the atomic planes. The fact 
that different Bloch waves generally have different 


absorption coefficients leads to the anomalous absorp- 


tion effect which is considered in the next section. 
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2.2.3 Anomalous Absorption Effects 


The principal reason for including inelastic 
scattering in Patera rena of diffracted beam intensity 
is to take anomalous absorpciongintoraccount. Pin order 
to miihisarate this effect, it is first convenient to 
examine the variation of the diffracted beam intensity 
with thickness in a perfect wedge-shaped crystal. 
Consider first a situation where aeeorption is neglected 


and where there are only two Bloch waves with large 
(i) 
g 


I, is found to vary sinusoidally with crystal thickness 


as shown in Fig. 4(a). Fig. 4(b) shows the same situa- 


excitation amplitudes Ceeac In this case, intensity 


tion except that both important Bloch waves have been 
assumed to have equal absorption coefficients. It can 
be seen of Fig. 4(b) that the mean intensity decreases 
with increasing thickness but the sinusoidal variations 
iLeno Lina pLescna. lhiG os termed normal absorption 

and simply decreases the scale of intensity .of the 
fringes mn thicky regions of the crystal... The visibility 


Vi 


OF Contrasts of the fringes, as measured by (i ae aan 


Glee laa isgnotuchanged.. (Pig. 4(c) son ‘thesoether 
hand shows the intensity profile obtained when one Bloch 
wave is absorbed much more strongly than the other. 


This profile illustrates the effect of anomalous 


absorption in a wedge-shaped crystal. In thicker 
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regions, the intensity variations eye pueeeren Thy 
dite deueee due to the strong attenuation of one of 
the Bloch'waves. However, considerable diffracted 
beam intensity is still present due to the weakly 
absorbed beam. Besides having an important effect on 
intensity profiles in perfect crystals, anomalous 
absorption has also been found to lead to important 


effects in the images of crystal defects. 


2.3 Relativistic Corrections to the Dynamical Theory 


Since the velocity of the imaging electrons in 
the eieeteen microscope is approaching the speed of 
light, it is necessary to take relativistic effects 
into account in the dynamical theory. These effects 
have been gonsidéréd by Fujiwara (1962), who developed 
a relativistic dynamical electron diffraction theety 
using the Dirac wave equation. His results showed that 
the non-relativistic theory developed using the 
Schrodinger wave equation can be corrected for rela-' 
tivistic effects by employing two simple sabe tieusions: 
These consist of replacing the non-relativistic wave- 
length by a relativistically corrected one and multiply- 
ing the Fourier coefficient terms, by bY a relativistic 


mass correction term 


me (2.31) 
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Here v and c are the velocities of the electrons and 
light respectively. The validity of these corrections 
has been confirmed experimentally by Hashimoto (1964), 


Dupouy et al. (1965) and Goringe et al. (1966). 
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CHAPTER 3 


APPROXIMATIONS INVOLVED IN THE MANY-BEAM, TWO-BEAM 


AND KINEMATICAL THEORIES 


Seleeelittr OGUC ULON 


The dynamical theory was formulated in the last 
chapter in terms of the eigenvalues and eigenvectors of 


an infinite matrix 
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- This matrix is infinite because an infinite number of 
diffracted beams were taken into account in the Bloch 
wave function of equation 2.7. In order to carry out 
calculations using this theory, it is necessary to 
approximate this Infinite Matrix by wa omatnixson asia ce 
size. This is accomplished in the dynamical theory by 
assuming that only a finite number of diffracted beams 
are important and need to be taken into account. In 
this way, the many-beam and two-beam dynamical theories 


can be developed. These two theories and the reduction 
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of the two-beam theory to the kinematical theory will 


be considered in this chapter. 


3.2 The Many-Beam Dynamical Theory 


In the many-beam dynamical theory the number 
of diffracted beams included is reduced to a finite 
number N. The question which immediately arises is 
which diffracted beams are of sufficient importance 
tenbe takensintdor account.: «Thisiisean important 
question because the computing time required to 
diagonalize the A matrix increases as the Square of 
the number of reflections included. The problem is 
generally solved by deriving some criterion for indica- 
ting the relative importance of the different reflec- 
tions. Reflections are then included in order of their 
decreasing importance until a point is reached where 
additional reflections make no significant change in 
the.diffracted beam intensity. Since the diffracted 
beam intensity has converged in this manner, the many- 
beam theory can be considered to be an exact theory 
because no Sa udidueant approximations are made by 
neglecting the remaining beams. 

The effect of a particular diffracted beam on 
‘the diffracted beam intensity is found to vary with 


crystal,onientation. »Forrthistreason; thenquestion of 
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which beams to include can be conveniently divided - 
into two situations: the systematic case and the 
non-systematic case. These two cases will ne con- 
Sidered in this section and typical many-beam 


intensity profiles will then be discussed. 


3.2.1) The ssystematic Case 


When a high energy electron beam is incident 
on a crystal, it is generally found that a number of 
different sets of lattice planes are close to satisfy- 
ing their Bragg: conditions. It is possible, however, to 
orient the crystal so that the only important diffracted 
beams result from different order Bragg reflections from 
the one set of planes. This situation is termed the 
| systematic case. An example of a diffraction pattern 
obtained under systematic diffraction conditions is 
shownsin tid eo. Ine slow or: rer lections 1s) called the 
systematic row and the reflections are termed systema- 
tic reflections. 

In the systematic case, the reflections are 
generally included by first including those reflections 
which are on or inside the Ewald sphere (see Fig. 5). 
Additional reflections are then included in pairs, 
one from each end of the row already included, until 


additional beams do not result in any significant 
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Fig 


An experimental diffraction pattern obtained 
under systematic diffraction conditions in 
molybdenum. The micrograph shows, from left 
to right, the reflections (000) through 
(550) "of the (110) systematic row. The 
retlegerwon (550) i5 in its Bragg condition 
and the reflections (000) through (550) are 
either on or inside the Ewald sphere. The 
higher order systematic reflections lying 

CO heeeeent and wert of the portion of eae 
diffraction pattern shown are outside the 
Ewald sphere. 


pcnangesin thesdittracted beam intensity. “The justifi- 
cation for this procedure is based on the fact that 

the effect of a systematic reflection h is found to 
decrease both with the increasing order (1.é€. increasing 
magnitude of the reciprocal lattice vector h) and 
increasing deviation from the Bragg condition. Gen- 
erally ten to twenty reflections are required in a. 


-Many-beam calculation for the systematic case. 


3.2.2 The Non-Systematic Case 


The second of the two situations which can occur 
in determining anion beams to include in a many-beam 
calculation is ease the non-systematic case. In this 
more complicated situation, the crystal is assumed to 
be at an orientation where a number of lattice planes 
are close to satisfying their Bragg conditions and, 
therefore, a large number of diffracted beams are impor- 
tant. ° A typtcal drifraction pattern (for the one 
systematic case is shown in Fig. 6. In such a pattern, 
one row of reflections iS considered to be the systema- 
tic row while those beams which are not colinear with 
this row are termed the non-systematic reflections 
(Hoerna, 1956)% 

In the non-systematic case, all of the important 


reflections will not necessarily lie in the same plane 
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A portion of a typical diffraction pattern 
in the non-systematic case. Many diffracted 
beams are important and must be taken into 
account in calculating the A matrix of the 
dynamical theory. 


or zone. The criterion for ordering reflections which 
was used in the systematic case is therefore not 
applicable. This is because the question arises of 
whether or not a low order reflection far from its 
Bragg condition is more important than a high order 
metlection which’) 1s in another zone end,close to its 
Bragg CONdicLoll. = Cann (19/5) Nas studied this problem 
and has developed a criterion based on Bethe's second 
order approximation (Bethe, 1928). In Cann's theory, 
the effect of a reflection h on the diffracted beam 
intensity of reflection g, where S5 is small, is given 
by 
Ue U 
cone Bela “fh ; (SZ) 


h 2KU 58h 


By using this criterion, reflections can be included in 
order of decreasing En until additional beams do not 

result in any significant change in the diffracted beam 
intensity. Generally thirty to one hundred reflections 


are required. 


36263, The Di fiimacted) Beam gdintensity sin, thes Many-Beam 
Theory 


48 


The expression for the diffracted beam intensity - 


in the many-beam theory follows directly from the inten- 
sity expression derived in Section 2.3 where an infinite 


number of diffracted beams were considered. When only 
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N reflections are considered, only N Bloch waves are 
excited in the crystal and equation 2.30 for the 
diffracted beam een ere iecones 
Detar Ni) 
= hR 0 (2) | a lies Fa 
(3 33) 
The form of the diffracted beam intensity profile 
calculated from he equation can be conveniently 
examined with reference to a plot of the dispersion 
surface in the many-beam case. The first seven branches 
of a dispersion surface calculated by taking thirteen 
Systematic rer lections into acGount are shown, 1nAaPlg. /. 
There is one Bloch wave associated with each branch 
of the dispersion surface, and, in general, each of 
these Bloch waves has a different value of ee It 
is therefore apparent from the intensity expression 
3.3 that the beating between the pee components of 
these Bloch waves can lead to a diffracted beam inten- 
Sity protile whichy\is a complicated. function of tchystall 
thickness. 
lives important sto note, however, that Braise 
not the number of Bloch waves, but rather the number 
of important Bloch waves, which is most significant 
in determining the form of the intensity profile... The 


only Bloch waves which will have a large effect on the 
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diffracted beam intensity at thickness z will be those 


(i) 


for which the terms Oarony 


. exp (-21qg 


z) are large. 
There ‘are two situations which generally occur. In 


a) 


O 


the first case, the terms C CF axpla2mq\ty.2) for 

two of the Bloch waves are significantly larger than 
the corresponding terms for the other Bloch waves. 
Theee two Bloch waves, termed primary Bloch waves in 
this thesis, have a significantly larger effect on the 
diffracted beam intensity than other Bloch waves. The 
form of the) intensity profile for this situation is 
shown in Fig. 8. The beating between the two primary 
Bloch waves leads to the large variation in diffracted 
beam intensity with thickness. The spacing between the 
intensity peaks is approximately given by ih ginal re 
which is the reciprocal of the difference between the 

z components of the wave vectors of the two eery 
Bloch waves. Othe: Bloch waves, however, can have some 
effect. If the effect is small, they introduce a 
variation in the peak to peak spacing as shown in Fig.8. 
The Bloch waves causing this variation will be termed 
secondary Bloch waves in this thesis. 

The second profile form which is commonly found 
in a many-beam calculation occurs when one or more of 
the secondary Bloch waves has such a large term 
Ca) 


exp (-21q Z) thet. tt Can be classified as a 


Be) a 
o “g 


primary Bloch wave. When this occurs, there will be 
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A diffracted beam intensity profile typical 
of the case where there are two primary 
Bloch waves and other secondary Bloch waves 
present. 
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three or more primary Bloch waves. A typical profile 
for this case is shown in Fig. 9. If more than two 
different wave vectors We ay are involved, then 

the beating between the three or more primary Bloch 
waves leads to an intensity profile which exhibits a 
much larger variation in peak to peak spacing with 
thickness than occurs when there are only two primary 
Bloch waves. Profiles of this second type will be 
referred to ingthis thesis ascomplex profiles. 

In order to indicate the manner in which the 
diffracted beam intensity varies with crystal thickness, 
it is convenient to introduce the terms extinction 
distance, average extinction distance, and variation 
in extinction distance. The extinction distance is 
defined in this thesis to be the peak to peak spacing 
in plots of the diffracted beam intensity versus 
crystal thickness. (i showl dibemnoted that this 
definition differs from that used by some authors who 
define the extinction distance to be the constant ao 
The average extinction distance is the average peak to 
peak spacing over some specified range of crystal thick- 
ness. The third parameter, the variation in extinction 
distance, is defined to be the difference between the 


largest and smallest values of the extinction distance 


over a specified thickness range, expressed as a 
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A diffracted beam intensity profile typical 
Oreille, Case wicre Chere are Uiresé or More 
primary Bloch waves present. 
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percentage of the average extinction distance. By 
using these three terms it is possible to convey 

some indication of the manner in which an intensity 
profile varies with thickness. In Fig. 9 for example, 
the percent variation is large (~100%). This indicates 
that the intensity profile is complex and that more 
than two primary Bloch waves are present. The percent 
VablLation~2ObL tne profude inygbaig. 8, where there are 
two primary and other secondary Bloch waves, is much 
smaller (~20%). In the two-beam theory considered in 
the following section, only two Bloch waves are present 


and the percentage variation will be found to be zero. 


3.3 The Two-Beam Theory 


The number of diffracted beams included in the 
Many-“beam theory generally léads tora Large A Seer 
which can only be diagonalized numerically. in order 
to Fae eate examining the dependence of the diffrac- 
Lsdabeatl inlensaty) on the Al matrix elements; 11S of 
interest to reduce the A matrix to a size which can 
be diagonalized analytically. Fukuhara (1966) has 
shown that,| for certain, specific crystal orientations, 
it is possible to develop analytical expressions when.- 
aS Many as twelve beams are included. Simple eigenvalue 


and eigenvector expressions, valid for all crystal 
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orientations, can only be found, however, for the case 
where two reflections are included. In this two-beam 
theory, the directly transmitted beam and the diffrac- 
ted) beam. g; are,.included, and«the eigenvalue equation 


ya becomes 


U 
(n) (n) 
Oy oe Co 
U Oh eS 5.40) 
g ee) (n) | 
2K roti Cg 


The eigenvalues of this equation can be found by 
requiring. that the determinant of the left hand matrix 


be zero. This leads to eigenvalues given by 


(ay-th n ty/med ev 2eqee 
Y 5 [s, Sgt Ue 
(379) 
CZ) ee Ve Diese 
and Y aaa [s, ae ue | : 


By substituting each of these eigenvalues back into 
equation 3.4, the matrix C- of normalized eigenvectors 


‘can be written 


f ab (l+ w')*- w | 
on 4 y (3.6) 
C eat 

w—- (1+ w) ? ae 
N N 


where 
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and the normalization factor 


Nee Vl e(H IU ES es fe a (3.8) 


This matrix can then be reduced to 


cos (8/2) sin(8/2) 
Cs (3.9) 
-sin(8/2) cos (6/2) 
by making the additional substitution 
w= cot B. (3.10) 


If these analytical expressions for the two-beam 
eigenvalues and eigenvectors are now substituted 
into the many-beam intensity expression 3.3, then the 
intensity of the diffracted beam g, neglecting absorp- 


tion, can be shown to be 


2 
U 

j= = sin (7s Z) (3641) 
g Ks ert 

eff 

where 
a 2 oe es P 

Sore 7 (S5 “- oc. Jute f= (3. dal) 


The form of the two-beam intensity profile can conven- 
iently be examined with reference to the dispersion 
SunLace plot -of Fig..L0.. Sance only twoirefections 
have been included in the A matrix, there are only two 


branches of the dispersion surface, each branch being 
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associated with one of the two Bloch waves present. 

The beating between eo. components of these two 
Bloch waves leads to an intensity profile in which 

the extinction distance is constant as shown in Fig.1l.. 
The extinction distance in the Pre shean theory is. 
equal to the reciprocal of the spacing between the 


dispersion surfaces and can be calculated analytically 


to be 
( ) ( ) U Ss 
Y -Y ae ie =o} eff 
Gx 


When =f is; zero!’ (i.e., 6= 1.0), this extinction distance 


is equal to BS while for large Pe (i.e. large 6), 


CoB 


from the dispersion surface plot of Fig. 10. since 


approaches 1/8,- This latter property is apparent 


the distance 55 can be seen from Fig. 12 to be equal 
to the distance between the spheres of radius K 


centered on o and g. 
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The kinematical treatment is the simplest 
approach to the problem of calculating the diffracted 
beam intensity but it involves two rather drastic 
assumptions. These assumptions are: 


as) An electron can only be scattered once; and = 
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A typical two-beam intensity profile 
showing a constant peak to peak spacing. 
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Fig.12. A geometrical construction showing that the 
distance d between the two K spheres centered 
Sie anaemic mequat Os Sias /O1nCG dais 
parallel’ tovs-Vanderc =Age K,vATCoy canrbe 
shown to be a parallelogram. and s_=d. 
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C2) The depletion of the incident beam, as it 
progresses through the crystal, can be 
ignored. 

The kinematical theory is usually derived in 
one of two ways. . The first is the geometrical method 
(see for example, Hirsch et al., 1965) in which the 
scattering by each individual scattering center on the 
lattice’ 1s considered. ~in this method, which is 
analogous to the X-ray case, the scattered intensity 
is calculated by taking the phase relationships 
between waves from different scattering centers into 
account. The second approach is a quantum mechanical 
approach in which the first order Born approximation 
PemUsed tomo Onve stile soci OdMmngerewavVewed Uagt1one Zui. 

A good description of this second approach is given, 

for example, by Vainshtein (1964) and also by Gevers 

CEs OE. 

The approach to be taken here, however, will 
be to examine the approximations which can be made in 
onder -toereduce ther very nearly exact many—beam theory 
to the kinematical theory), The approximations which 
are required have previously been given by Serneels and 
Gevers (1972). These approximations are first, to 
reduce the many-beam A matrix to a two-beam matrix, 


and second, to calculate the eigenvalues and eigen- 
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vectors of this two-beam maeioms diddy to first 
order. However, their approach is very involved 

and therefore a simpler approach has been adopted 

in this thesis. In the approach used here, the 
eigenvalues and eigenvectors of) thegj-ewo-beam A 
Matrix will’ bescalculated correct ‘to first. order 

by employing the approximate methods of matrix 
diagonalization presented in Appendix A, 

| If the procedure given in Appendix A is used 


to diagonalize the two-beam matrix 


“g 
: 2K 
A= rey 
Ug 
2K Sg 


then the first order eigenvalues can be found from 


equation Ac2 to be 


ee) ee EG OP ih Gee (3.14) 


Similarly, the matrix of eigenvectors, Correct to first 


Orders Can be Lound from Equations A.Y and ALi co be 
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LE-these,f£irst-ordersexpressions for the eigenvectors 
and eigenvalues are substituted into the INany-beam 
intensity expression 3.3, then, neglecting absorption, 
the intensity of the diffracted beam becomes 


Use 22 
J (4) sin? (=n (y 2) -y (2?) 2) 
g 


(ee) sin® (5,2) - 7 (36>) 


This expression is the kinematical expression for ene 
diffracted beam terensieye 

The kinematical theory will be a good appro- 
ximation to the many-beam theory if two conditions are 
satisfied. Firstly, the effects of systematic reflec- 
‘tions other than o and g and of non-systematic 
reflections must be small, and secondly, the Liest 
order eigenvalues and eigenvectors must be good 
approximations to the two-beam eigenvalues and eigen- 
vectors. The conditions under which this second 
approximation holds can be examined by considering 
the approximations used to Aegonalawe the two-beam 
Se 


A matrix, The fiest order approximation, (y 


ae 


used in the intensity expression 3.16 will be a good 
approxi Mat Ol wae ciemd1therence, between mec and 
CA te to second order is small. -It can be seen 


from the eigenvalue expression A.14 in Appendix A that, 
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correct to second order, 


U. 2 
cy (1), (2)) aaah eae = (53) ? (2 


g 
Accordingly, the first order eigenvalue difference 


will be a good approximation if 


fe Pee fie ee (3.18) 


This same approximation is also used in calculating 
the orthogonal matrix of first order eigenvectors in 
equation 3.15. The kinematical theory can therefore 
be seen to be a good approximation providing only two 
beams are important and providing the inequality 3.18 
holds. 

The form of the kinematical theory intensity 
PLrOolLvie can bewseen by examining the intensity _ 
expression 3.16. The diffracted beam intensity will 
vary sSinusoidally with thickness with a constant 
extinction distance of 1/55 The way in which the 
kinematical extinction distance varies with the 
deviation parameter 6 can be examined by referring to 
the kinematical dispersion surface which is obtained 
from the kinematical eigenvalues O and ey. CoOL GortG 
was shown in Fig. 12 that the distance d between spheres 
of ea sen centered on the reciprocal lattice points o 


and g is simply So° the branches of the kinematical 
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dispersion surface can be represented by these spheres. 
These spheres are indicated on the two-beam dispersion 
Ssuneace: DLOtL.OL M10. 10. Since the kinematical extinc- 
tion distance is 178 or the extinction distance can be 
seen to be infinite for 6 = 1.0 Se = 0) and to approach 
the two-beam extinction distance with increasing 6 


(increasing Fels 
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CHAPTER 4 


COMPUTATIONAL METHODS 


ie Lae nertoacuctLLon 


The research presented in this thesis relies 
extensively on numerical calculations based on the 
theory presented in the last two chapters. It is 
therefore important to indicate the procedure used 
) BOsGCalculatesand diagonalize. the Asmatrix, (see equation 
Zoab,,eoeCelonez.). 1) sandsto calculate the wittiracted 
beam intensity and extinction distance. These proce- 


dures will be discussed in this chapter. 


4.2 .Methods Used to Specify the Crystal Orientation 


In Order tovtcarry out calculations using the 
dynamical theory, it is first necessary to obtain the 
elements S}, on the diagonal of the A matrix (see 
@quatzron +2175 Seceron Jee It can be seen from 
Fig. 2, pagé’' 22, that the values of these elements 
depend on the coordinates of the tie-point, i.e. the 
center of the Ewald sphere. This in turn depends on 
orientation of the crystal with respect to the incident 
electron beam. In the research carried out for this 
thesis, two different methods were used to specify fhe 


CEYSualtOLreneation.  Inwther first method, -déveloped 
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by Canny (L973), an initial crystal,orientation was 
obtained by specifying two reflections which were 

in their Bragg conditions. This information was then 
used to calculate the coordinates of the tie-point. 
The orientation for which the calculation was to be 
performed was then obtained by changing the direction 
of the wave vector K where K is the .vector «from -the 
Cle sponte to the syorigin sof the reciprocal slattices 

In this manner, new coordinates for the tie-point were 
obtained. These coordinates were then used to enV 
late the s-=values.,...in using,this method it was Lound, 
however, that significant numerical error could arise 
in the numerical procedure used to calculate the 
coordinates of the tie-point. AS a result, it was 
very difficult to specify that a particular reflection 
was, in its exact Bragg condition. This was.required, 
for example, in order to obtain dispersion surface 
shapes accurately and to carry out the critical voltage 
calculations discussed in Chapter 8. 

In order to overcome the problems associated 
with Cann's method, a second method was developed for 
specifying the crystal orientation. In this method, 
referred £00 in this sthesiseas the tie-poinl mernod, 
the crystal orientation is’ specified by the coordinates 
CEBENe DON UG Wiich ws the perpendicular projection 


ot the tie-point onto the zero order Laue plane (see 
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Fig. 13). The s-values are then obtained from the 
eee Ofapointiiledsee Section 430) —This 
method was found to be very accurate because no 
numerical error was involved in specifying the co- 
ordinates of the point L. 

it should be noted that, if only systematic 
reflections are considered, complete specification 
Oletne croystalyocientation 18 not. required. “in the 
systematic case, the deviation of the reflections 
from their Bragg conditions can be specified by the 
perpendicular projection of the tie-point onto the 
systematic row (see Fig. 1, page 12). As indicated 
in Fig. 1, these deviations will be indicated in this 


thesis by the parameter 64. 


aoe The Calctiilation,of the Diagonal Elements of the 


A Matrix 


Once the crystal orientation has been specified, 
Che,sorientalion dependent diagonal velemen cS Ora thean 
matrix can then be calculated. The procedure used ZO 
calculate these elements when the position of the Ewald 
sphere is specified using Cann's method has already 
been documented (Cann 1973) and will not be discussed 
agai nere. wit aseor interest in this thesis, however, 


to consider in detail the methods which can be used to 
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Ewald Sphere 


The intersection of the Ewald sphere with 
the zero order Laue zone. The crystal 
orientation is specified in the tie-point 
method by the point L which is the perpen- 
dicular projection of the Center Gomeune 
Ewald sphere (tie-point).onto the zero order 
Laue plane. 


ot 


® 6 
a ‘ ~. 
ssenqe biows / 


/ 


_ { 
\ 
° & 


isiy stedge biswa sit) lo GL F984 TSE ant 
Pete ine iT .onos Syed aabaa Oras ont - 
Seis ae Seed ae er as gpeetaebibewsiond 
ae ache tol age rg pe aot te 
Ed reg ak Mgt eo 


2 . aa 
‘8 


43) i 
we? Mic Cone h = Bow . 
i a Prins iebn: 


7a: 


calculate the diagonal elements when the position of 
the Ewald sphere is specified using the tie-point 
method. Three possible methods, termed the s-value, 
t-value, and exact methods will be considered and 


compared in this section. 


4.3.1 The s-Value Method 


When the A matrix is derived by using the methods 
discussed in Section 2.1.1, then the elements on the 
diagonal of gthée A matrix ane, termed, ssvalues.~ The s- 
value for a particular reflection h is defined to be 
the distance in the z direction from the reciprocal 
lattice point h to the Ewald sphere. When the position 
of the Ewald sphere is specified using the tie-point 
method, then these s-values can be calculated with the 
ald %otremig. 14, in this Ligure, the z-axis has been 
taken to be perpendicular to the crystal surface and 
the symmetrical Laue case has been assumed. The point 
L is the projection of the tie-point T onto the zero 
order Laue plane. It can be seen that the Sh value 


for some reflection to the nuh Laue zone is 


Sh = LN - SNe nL Negetoge ei ee Ph 
2 pe 2 jue 
= (kK o> Th) - (K" - TT) TER), (cae) 
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Fig.14. The dispersion surface and Ewald sphere 
construction for calculating the diagonal 
elements of the A matrix. 
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where Pi is the distance measured in the z direction 
Lrom the 2ero Order Laue™ zone to the n‘? order Laue 
zone. Although this expression for Sy is analytically 
COLEecce, Signiticant numerical error can aie: Cen 
the two large square root terms are nearly equal. 
Under these conditions, increased numerical accuracy 
can be obtained by using a Taylor series to expand the 
Square root terms. If. terms correct to second order 
are retained, then Sy becomes 

— 


ate - i 
Sh 5e (2G Tt) + 


(4.2) 
re) 8K? 


ee 
h of a ee 
This method of calculating the diagonal elements of the. 


Plate ixe (i.e. equations4 1) and ferasmall s -equation 


lay 
4.2) will be referred to as the s-value method. 

fee USsOteincerestetounoue Chat ondy “tlie 
first order term of equation 4.2 is retained, then 


the s-values for reflections g in the zero order zone 


become 


2 
Ce i (4.3a) 


ii 
K g(t - ey) 


=~ g AQ (4.3b) 
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wnere AG’ vs ‘defined’ in Fig. 14. This first order 
approximation is used by many authors (see for 


example Head et al., se erRcth 


4.3.2 The t-Value Method 


It is next of interest to consider an alter- 
F 72 (1) 2 . t 
nate expression for the terms (K™- (Kp )~) in Bethe's 
equation 2.9, page 20, which in turn leads to a new 
expLession for the diagonal ,clements jf sLhe. Avmatrx. 
This alternate approach is similar to one developed by 
Miyake (1959) for reflections in the zero order Laue 


zone. The expression to be presented here, however, 


can also be applied to reflections not in the zero 


order Laue zone. If the eigenvalue ve is defined 
to be 
CPG Peer) Ke 
Y = (Ko ‘bg K (4.4) 
then it can be seen from Fig.. 14 that 
2 a(n 2 i (i) 2 eee 
fea EDS re ci heect pP,} ae 
. Wars (i) Bil be 
=K* - {K + (y + py) Th, 
BOHR rey irs pastes PPS ey 
% ¥ Ph h h 
Gltd nept ject: 
te (a5) Bene 
Sh ORs aaa) Lick a8 Sie Belt e)) 


ST ; vy 


' 
$ ; 


sebro gavkiR 1d? <*L -prt ab Beatteb et 84 oxeeiw ae 
+G2 soe) axodjus yaem yo bean ef so itamixouggs ; a 
EN@1 ,.f5 79 best silgmexe a 


bosseM evylseV-3 ont 4,&.% 7 


~yosis np xsbienos oF Jaststnt Fo Jxon ef St _ 
| _ jb be Se | = 
a’ oAargee arr t 4 4 \}) -"*) efxet edad tot noilsesrqss oJsen : 


wen 5 of 2bss! oxvt ni fold ,Of spag ,&.& soresupse 7 


_xiatem A ods to atnemeis Isnops to sns 103 aoLeesiqxs 
yd beqolsveb ono oF sHlimie ef Nosotdas oteasedis eiAT 
sug] 1sbi0o oes’ sit at anotsssiges sor (ede) SAByiM - 
,teavewod ,Sisd Bestnsesig Sd 02 noteesiqx2s Saif .Ssaox of 
oss oot ni ton enoidsveliex af bettqgs ad ogis nso 
benijab el (4) , sulévnspis of3 TI antes seat tsb10 


“bd oF | (a 


ib. b) a Ronee: = (by ' 7 as 


° _— 


tet bl pit most ayse od aso te useRs 
a - *t 4 + tw i= y= (ON 8a : .* 
_ a 7 a 


a= ae 


7 
<< 


4a oa at (6 eae} +H} = : 
A ao i Ts : ; —_ 


- es (oan — 


% = : : a om a 


Tf “LteiSsassumed sthat 


i, tT Pr 
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then 
Me (a) ee (i) 
K (ky ) 2K (y + t)) (407) 
where 
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Th. 
Cys Se ON) Ph ° : (4.8) 


Tf Equation 4.7 is substituted into Bethe's equation 


2.9, page 20, then the eigenvalue equation 


ag ) tdogdile d) falas 


can be obtained. In this equation, the off-diagonal 


elements are given by A = U siq Es the diagonal 


gh 
elements are the t-values of equation 4.8. The eigen- 
values and eigenvectors of equation 4.9 can then be 


used to Calculate the diffracted béam intensity. 


4.3.3 The Exact Method 


The third method which will be used to; calcu= 
late the eigenvalues and ergenvectors Of tne sAsmatrax 
is an exact approach when has not yet appeared in the 
literature. . This. méthod considers-reflections.in the 
zero Laue zone only, but it does not make the approxi- 


mations inherent in the s-value and t-value approaches. 
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(If reflections other than those in the zero order zone 
rane considered, PehenoBethe|ts,equation cannot be cast 
in eigenvalue. form.) In this approach, the term 
2 eri) 5 aoe 
K Poiks )~ is calculated from Fig. 14 in the follow- 


ing manner: 


2 h(i) eee ae hl (Vy ez 2 
K'- OSs ) ae i oe ((ko ie + ge 
= (K2 = (k (4) )7) 4 (2 - 72) 
Zz O Z O g 
x a (a) (1) Coe 
TAA (Ko PORES vor 6 Wea, uae bes (4.10) 
1Ese y 2) is defined in the same way as for the s-value 


method, i.e. 


pe (a) - 
Y = (ko \e K, (ee a) 
then 
a ets Rupee 
Bae ee eee ON 0 
2K zy K 2K 2K 
eli 
etary) O g 
= —T + oI Gacgl2)) 
where 
; : K (5) 
<(ea) Sg levee Ela) nae Y 
I = Vf (= + 3K) . (aes) 


If the terms (cola) 72k are used as the diagonal 
) 


elements of the Asmatrix and the eigenvaiues bee 
are calculated, then the exact eigenvalues y ) defLaned 


by equation 4.11 can be found from equation 4.13. The 
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correct eigenvectors are those found by diagonalizing 
Che 2 Matrix. For retiections in the zero order Laue 
zone, the exact eigenvalues of equation 4.4 can be 


calculated in a similar way from equation 4.5. 


4.3.4 A Comparison of the s-Value, t-Value and Exact 


Methods 


Because the exact method is a new approach which 
does not make the approximations inherent in the s-value 
and we-Value methods, itis possible for the £irst time 
to assess the accuracy of the two approximate methods. 
This was done by comparing the diffracted beam intensity 
calculated using each of the approximate methods with 
that found using the exact method. These calculations 
showed that for thin crystals the profiles were quite 
Similar. As the thickness increased, however, signi- 
ficant differences in the diffracted beam intensity 
were found although the profile shape remained very much 
the same. An example of these differences is shown in 
Fig. 15. These calculations’ were carried out for the 
(110) dark\ tied deinymolybdenum at “OUR KV. 1 can (be 
seen from this figure that the peaks in intensity 
obtained using the t-method are displaced slightly 
from those of the s-method and the exact method. fThis 


leads to significant differences in the diffracted beam 
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Fig. 15. {Asconparison oh the ditfiracted beam 
intensities found an thick crystals 
using the t-method, s-method, and exact 
method. Thirteen reflections of the 
(110) systematic row in molybdenum were 
considered and 6 = 3.5. 
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intensity ata particular crystal thickness. For 
example, at a thickness of 7600 A the intensities 
foundyby (che 4t-method and exact method drtiersby 8°... 

As can be seen from Fig..15, good agreement was obtained 
between pthe ts-methed hand xsexactymethod tefonethicker 
crystals. 

The reason that the t-value method is a poorer 
approximation than the exact noened can be seen by 
comparing the two methods analytically. This, LS most 
easily done by.first»showing that.the t-method and the 
first order s-method approximation of equation 4.3a 
are equivalent. This can be shown by noting that the 
differences between the eigenvalues of a matrix will 
notebeychangedeif asconstant is added to,each of the 
diagonal elements before diagonalization. Furthermore, 
thenaddittiont#ofisthis scconstantywwili not change the 
eigenvectors. If the constant 12/2K is added to all 


the ey elements, then the t's of equation 4.8 become 


(4.14) 


For the zero order Laue zone. This expression for the 
diagonal elements is identical to the expression 
obtained when the s-values are calculated correct to 
the first order (see equation 4.3a). Accordingly, 


the eigenvalues and eigenvectors calculated by using 
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the t-method are of the same accuracy as those calcu- 
lated by using the first order s-value approximation 
of equation 4.3a. 

The intensity profile differences which were 
found when the t-method was compared to the exact 
method can now be examined by comparing the first 
order s-value approximation to the exact method. 

Since the same diagonal elements are used jinieectie A 
matrix in each case (compare equation 4.3a to equation 
4.12), then the eigenvectors for each case will be 
identical. However, the eigenvalues obtained by 
diagonalizing the A matrix are used, without alteration, 
in calculating the intensity in the first order s-value 
case, while these same eigenvalues are corrected (see 
equation 4.13) before being used in the exact calcula- 
tion. wAccordingly, the’ eigenvalues of the first order 
s-method will differ slightly from the exact eigen- 
values and this results in the slight displacement 

in the peaks in the intensity profiles an. thick panes 

of the crystal (see Fig. 15). When a calculation is 
carried out using the more accurate s-values of equations 
4.1 or 4.2, then the eigenvectors will no longer be 
exact. The eigenvalues, however, are more accurate 

‘and generally, improved agreement with the exact method 


is obtained. 
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For the diffracted beam intensity calculations 
presented in this thesis, either the s-value method 
(equations 4.1 or 4.2) or the exact method was used. 

The exact method is preferable when only reflections 

in anew ed order Laue zone are considered because it 

is the most accurate and yet it does not take any more 
computing time than the s-value method. (The extra 

time required to correct the eigenvalues is compen- 
sated for by the simpler form of the diagonal elements 
GE CiewA fmatenices. Wslowever the et=methodawds Sised for 
plotting dispersion surface plots. The t-method is 
useful in this case because the eigenvalues and diagonal 
elements are defined with reference to a plane at dis- 
tance K ingtherzrdiarection fromthe zerosorder lauewzone. 
This allows the dispersion surface and the K spheres 
about the reciprocal lattice points to be obtained by 
plotting the eigenvalues and the Ey elements respec-— 


tively. 


"4.4 The Calculation of the Off-Diagonal Blements of 


Lie. A Matias 


The gh ,oft-diagonal element of tne Ay matrix is 
equal to Ufa ae where, for a crystal composed of only 


one atom type, 


eck Pasa del, (-B|g-h| 774) (4.15) 
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(See for example, Howie, 1970.) . In this expression, 
bersethe relativistic: mass correction factor for the 
incident electrons, be is the volume of. fhe unit cell, 
Dapaa is the kinematical Structure iactor, and 

f (sin yaene! is the atomic scattering factor for 
electrons. The scattering factors used in this thesis 
were those calculated by Doyle and Turner (1968) using 
relativistic Hartree-Fock atomic wave functions. 
Equation 4.15 also contains the Debye-Waller correction 
B. The B values used in this thesis were 0.5l, 0.20, 
and 0.55, for gold, molybdenum, and aluminum respec- 
tively] (Ibers Pen Vainshtein, 1962). PoOred particular 
crystal structure and accelerating voltage, the values 
(ope Uren tend to increase with atomic number and decrease 


with increasing |g-h| as shown in Fig. 16. 


40 elhe Calculatvon of the A Matrix 


iM Order wovvake avbsoLlption Anco account, it) ius 
Necessary #Oncalculate the Cerms Uji /2K Che lhiewmn, 
matrix. The values of Mie which were used were 
obtained from the ratios Dele given by Humphreys and 
Hirsch (1968). The values of U,/Ug for gold, molybdenum, 


and aluminum are plotted in Fig. 17. 
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Values ot U,saspastUnct 1LOngOiaG sion 
100 kv Bfecpnone in aluminum and gold 
(FCC crystal structure) and in molyb- 
denum (BCC crystal structure). These 
materials have atomic numbers 13, 79, 
and 42 respectively. 
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Values of Ug/Ug for 100 kV electrons obtained 
from the data presented by Humphreys and Hirsch 
(1968). The dashed lines are extrapolations 

of their results and were used to calculate 
Ug/Ug for large values of g. The ratio Us 
was pale to be zero once these extrapolations 
intersected the g axis. 
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anOee Matrix UlegonallZation ald tie Calculatloln "or 


Diffracted Beam Intensity 


In comparing the kinematical, two-beam, and 
mahy—beam theories, the A matrix was diagonalized by 
using the subroutine EIGRS in the IMSL subroutine 
package (IMSL, 1974). This routine was used because 
it was found to be very accurate and ereguinreds less 
computing time than other available routines. 

Houatdon J. 35was then used to calculate the ditfracted 
beam intensity. The intensity was generally obtained 
for thickness increments of 1-10 A depending on how 
Bap ly ethics iicensiey (varied wiles thivekness. » eins order 
to improve the accuracy to which the extinction distance 
could: be calculated, IMSL (1974) smoothing routine 
ICSIVE was used to increase the number of points in 
the regions of the intensity peaks. The extinction 
distance was then obtained from the smoothed curve. 
The depth increment and the number of smoothed points 
close to each peak were determined in order to give 
aneaccunacy Of bétbtexr wuhan, Ul5o2 4m, Ccalculatang@eme 
extinction distance. 

In examining approximations used to take the 
imaginary part of the lattice potential into account, 
it was necessary to diagonalize both real and complex 
Matrices. imvthiias case, the SSP (1970) routine EIGEN 


was used to diagonalize the real matrices. (The IMSL 
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routine EIGRS was inconvenient because some of the 
eigenvectors of adjacent orientations were often 
multiplied by -l and this led to complications in 
plotting the parameters oe See Chapter 6.) 
The complex nee were diagonalized by using the 
PMoo (L974) subroutine EIGCC, 

All the computer programs were written in 
FORTRAN and executed on the IBM 360/67 computer of 
the Department of Computing Services at the University 


of Alberta. The calculated line graphs were plotted 


using an off-line model 925/1036 Calcomp Plotter. 
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CHAPTER 5 


A COMPARISON OF THE KINEMATICAL, TWO-BEAM 


AND MANY-BEAM THEORIES 


Sele plneroduction 


It was indicated in Section 1.3.1 that the 
validity of the kinematical and two-beam theories 
for large values of 6 will be examined in the first 
part of this thesis. This assessment will be carried 
Gu Leine thas elevates by comparing these two theories 
Louthesnearlyesexact smany—beam theory. The question 
arch immediately arises is what method should be 
Us@d sCO scarry out this comparison... .The method to be 
used in this thesis will be based on a comparison of 
intensity profiles obtained from each of the three 
theories. (These intensity profiles can be obtained 
experimentally by examining the images of wedge-shaped 
crystals.) This method was chosen because the shapes 
of these profiles are very sensitive to many-beam 
effects. (In this thesis, changes in the two-beam 
intensity profile which occur when additional reflec- 
tions are considered will be referred to as many-beam 
effects.) 

The next question which arises is what methods 


should be used to compare these profiles. The first 
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method to be used will be based on a comparison of 

the ence extinction distance in the many-beam 
profiles to the two-beam and kinematical extinction 
distances. The reason that the many-beam average 
extinction distance has been chosen is that, ifsonly 
two primary Bloch waves are excited, this parameter 

is very nearly equal to the reciprocal of the spacing 
between the dispersion surface branches for these 

two Bloch waves. This is analogous to the concept of 
extinction: distance in the ee veen theory which is 
equal to the reciprocal of the spacing between branches 
one and two of the dispersion surface.’ (A method 
Similar to this has been employed by Goringe et al. 
(1966) in comparing the two-beam and many-beam extinc- 
tion distances.) 

The comparison of the many-beam average 
extinction distance to the kinematical and two-beam 
extinction distances, however, is of limited value 
in assessing the validity of the two analytical 
theories... The reason ,for this lies in the |fact, that 
secondary Bloch waves will not have a significant 
effect on the many-beam average extinction distance. 
These Bloch waves, however, can result in significant 
variations in the many-beam extinction distance. If 


these variations are large, then the kinematical and 


30 doeiisqno> 6 a0 beesd Sd LLdw boav sd od bodiem 
mased-yosm oct ai soasteib, nokdoiitxe opstsvs SAF 


fAoitomiaxe Isotjameni®s bas mped-ows sft of aelitorg 


spsteve rissdl-yasm od Jett noase1 eff ,.se0nsJelb 

vino 2i ,tedd ei sorodo aged asd oonaseih aoltonisxe ; 

istems1sq aids ,bedioxs sis esvsw doold yusmiagq ows 

pitoege sit Io Isnoigives sdt oF Isups ylasen yiev at 

szeit 103i sxedpasiad spstave sotexsqeth sft neswisd 

to tyoondo sit of asopelans ef eidT .eavew doole ow? 

ei doidw yaosit mésd-owy odd ni oonetderh npiszoaigxe 
esronsid neewisd pnipage Sd3 to Issorgioe1. say oF Taupe - 

botdom A) .sostive aolbersqeib saris to Gwe baa sno 

.is +9 spnizod yd bevelqms asad ead eiat oF reliniz 
-nnitxe mspod-ynam bus meed-owd ofd parreqiied az (ad@T) 


(,e9nretelb aoiz 


Spsisvs msod-ynom odt to soatysqnes off a 7 : 
msed-ows bois IsvivementH sat o¢ sonaderh aokgonisxe a. 
ailsy bosimil io al ,xsvaword ,esonsdaib noljonitxs . _ 

lésistyishe ows ofit to ytibitev sat pateasezs ni a 


Bi, 5 essa} eens namaniaions .28izoorls 


ea TNA fet PO AD el nq 


a 


= Fy _ ¢ = : - y > 


two-beam theories are poor approximations even though 
there may be good agreement between the many-beam 
average extinction distance and the results obtained 
from the two analytical theories. Therefore, in 
examining the validity of the two aneletecal theories, 
the variation with thickness in the many-beam extinc- 
tion distance will also be examined. 

Many-beam effects can be due to either system- 
atic or non-systematic reflections. Accordingly, both 


these cases will be examined in this chapter. 


5.2 A Comparison of the Kinematical, Two-Beam and 
Many-Beam Theories in the Systematic Case for 


Large Values of 6 


In the first set of calculations, the extinc- 
tion distances obtained using the kinematical, two- 
beam and many-beam theories were examined as a 
function ;of the deviation,parameter .d.forpthes (10) 
dark field intensity in molybdenum. The accelerating 
voltage used was 100 kV and, depending on the value of 
6, from fourteen to sixteen reflections from the (110) 
systematic row were included in the many-beam calcula- 
tions. Molybdenum was chosen because it is of medium 
atomic number and the accelerating voltage of 100 kv 


‘is commonly used experimentally. The first order 
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dark field intensity was chosen because the contrast 
in this image is generally much better than in the 
bright field image and thus it is often used in 
experimental studies of eer an distance (see for 
example, Cann and Sheinin, 1974A). The validity of 
the two-beam and kinematical theories was assessed for 
values of 6 >3.5. Values of 6 < 3.5 have previously 


been investigated by Sheinin (1967). 


5.2.1 A Comparison of the Kinematical and Two-Beam 


Extinction Distances with the Many-Beam Average 


Extinction Distance 


An examination of the many-beam intensity pro- 
' files for the (110) dark field in molybdenum showed 
that, for 6>3.5, the profiles were of two basic types. 
The first type is similar to that shown in Fig. 18(a). 
It can be seen from this figure that the intensity 
profile is sinusoidal in character over the range of 
thickness considered. The second type of profile which 
was obtained is similar to that shown in Fig. 18(b). 
This-profilenisesinusGidal,1n Character uUp_toga, thick— 
ness of approximately 1000 A, DUG, in GOnlLLastecosuhe 
DroLtile: OfMbog eel sia), this second profile then becomes 
conplex, silly thicken parca Ole the ‘crystal, «the profile 


becomes sinusoidal in character again. It was found 
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Figs ls, Typical, (110) dark field intensity profiles 
in molybdenum: (a) 6 = 3.5, and (b) 6 = 4.99. 
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that profiles of the type shown in Fig. 18(a) occurred 
for all values of 6 except those close to 5.0 and 6.0. 

Aswindreated in Section 5.157 the valadity of 
the two-beam and kinematical theories will first be 
examined by comparing the extinction distances of the 
two analytical theories to the many-beam average 
extinction distance. The many-beam average extinction 
distance was first obtained for a thickness range of 0 
to 1000 A. This average extinction distance and the 
two-beam and kinematical extinction distances are shown 
in Fig. 19. It can be seen from this figure that the 
agreement between the three theories improves with 
increasing 6. Furthermore, the kinematical extinction 
distance is a considerably better approximation to the 
very nearly exact mMany-beam average extinction distance 
than is the two-beam extinction distance. 

In order to examine how the many-beam average 
extinction distance changed with the thickness range 
considered, calculations were also performed for a 
thickness range of 0 to 1500 A. The results are shown 
Ler LO. 2 0 It can be seen that the many-beam average 
extinction distance shown in Fig. 20 is very similar 
to that shown in Fig. 19 except for values of 6 close 
to 4.0, 5.0, and 6.0 where abrupt changes in the many- 


beam average extinction distance occur. 
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Fig. 19. The kinematical and two-beam extinction 
distances and the many-beam average 
exvinctron “distance for *the iO) sdark= 
field intensity in molybdenum, © The 
many-beam average extinction distance 
was obtained for a thickness range of 
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Fig...20.. .The kinematical and two-beam extinction 


distances and the many-beam average 
extinction distance for the (110) dark 
field intensity in molybdenum. The 
many-beam average extinction distance 
was obtained for a thickness range of 
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5.2.2 The Variation with Thickness of the Many-Beam 


Extinction Distance 


As pointed out in Section 5.1, the results 
presented in the last section are of limited value 
for assessing the validity of the kinematical and 
two-beam theories because variations in extinction 
distance with thickness in the many-beam theory are 
averaged out. Since the analytical theories will be 
poor approximations if the variation in extinction 
distance is large, it is next of interest to examine 
this variation. The percentage variation in extinc- 
tion distance was first obtained for the thickness 
zange of 706-1000 A. These results are shown in 
Fig. 21. it can be seen from this: figure that there 
is a variation in the many-beam extinction distance 
of approximately 37% at 6 = 3.5. This variation then 
decreases rapidly with increaSing 6 to approximately 
b0G at 6 == 7428, For values of 6 > 438 the: variation 
decreases much more slowly to approximately 5% at 
Oa 9s 0 3 

In order to examine how the variation dara 
extinction qadeence changes with the thickness range 
considered, calculations were performed for the 
thickness ranges 0 to 500 & and 0 to 1500 &. The 
percentage variation in extinction distance for the 


thickness range of 0 to 500 R was Very Similar to that 
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Variation in Extinction Distance: (3%) 
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The variation with thickness in the many-beam 
extinctionidistance for the (110), dark field 
intensity in molybdenum. The accelerating 
Vola Gey USed gwaSeel OOskKV sand jasvecnrySia lack 
ness range of 0 to 1000 A was considered. 


prersuce (z) 


o 
m 


qu EXpTVeETOU 


ES : 
=: 
ay 
- : 
7 
of g 7 
0 
' : 7 _ 
. 2 Shee 7 
ae hat | = 
Macd-yrem of ni agoniokds itiw noit@i Tey oft LS prt | 


bisiz daveb wenn eit 297 ihe a aed 
_pritsiolsous oat . | | 
ida. 


shown in DigG.acls) lhe results obtained for the thick-— 
ness range 0 to 1500 &, on the other hand, are quite 
different and are Shown in Fig. 22. It can be seen 
from this figure that large peaks occur for values of 


oO CLOSE) to 25. 0 ande6..0. 


5.3 Analysis of the Results for the Systematic Case 


The purpose of this section 1s to carry out 
an analysis of the results presented in Section 5.2. 
In this analysis, only the kinematical and many-beam 
theories will be compared. The reason for this can be 
seen from Fig. 19 which shows that, for large 6, there 
is good agreement between the kinematical and two-beam 
theortes., “This of course is a well-known result.(see 
for example, Whelan, 1970) and occurs because the con- 
dition |s,2K| >> [Ug | (see equation 3.18 of Section 
b4)) us tsatistred for tlarge 6, Tie. large Bae in 
comparing the kinematical and many-beam theories, 
the results for values of 6 where a systematic reflec- 
Elon, as closestoswes Bragg condition will be sconsidered 
separately from seems obtained when no uel ssiia eke 
reflection is close to its Bragg condition, he 
reason for this is that the shapes of the intensity 
profiles obtained for these two cases are quite 


different (compare Figs. 18(a) and 18(b)). 
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in Extinction Distance 


Varration 
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The variation with thickness in the many- 
beam extinction distance for the (110) 
dark field intensity in molybdenum. A 
ecrystalsthickness range of 04to 1500 A 
was considered. 
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5.3.1,.An Analysis of the Results for Values of 6 
Where no Systematic Reflection is Close to 


1ts Bragg Condition 


5.3.1(a) Numerical analysis of average extinction 


distance results 


The Tesules, presented) in Fig. 19 show that 
there is good agreement between the many-beam average 
extinction distance and the kinematical extinction 
distance when no systematic reflection is close to 
tesepragg sconditien. | inorder to explain this result; 
the number of primary and secondary Bloch waves in the 
many-beam calculation will be examined. (It should be 


(eg) 
O 


noted that, it 1s not the value of C but rather 


the value of C (i) 


et na z) which determines the 


™ exp (-21q 


contribution made by Bloch wave i to the diffracted 

beam intensity at thickness z. It was found, however, 
that Bloch waves which were unimportant at thickness 

zZ were also unimportant at the top surface. Accordingly, 
a good indication of the number of primary and secondary 
Bloch waves could be obtained by examining the excita- 
(eli) 


tion amplitudes C The magnitudes of the 


excitation amplitudes lie a LOM Ow= fon ye tl o.; 
625,. /20 and onomaLemshown iin Table le slteshould be 
noted that in this table the Bloch waves have been 


labelled using both the conventional method (see 


Section tae and an alternate method referred to in 
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Table l(a). 
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The magnitudes of the Bloch wave excitation 
amplitudes C{i)cli) for the many-beam 
calculation discuss sed in the text. The 
primary Bloch waves are underlined. The 
Bloch waves have been labelled by using both 
the conventional scheme and the alternate 
K-sphere scheme described in Fig. 23. 
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Table. l(b}. 


(i) 


The magnitudes of the Bloch wave excitation 
amplitudes C‘1/ci\1/ for the many-beam 
calculation “discussed in the text. The 
primary Bloch waves are underlined. The 
Bloch waves have been labelled by using 
both the conventional method and the 


‘alternate K-sphere scheme described in 


EerGn? 23). 
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Ewald Sphere 


Fig. 23. The K-sphere scheme for labelling the branches 


of the dispersion surface. In this method, 
which is similar to a method used by Cockayne 
(1972), branch h of the dispersion surface is 


that branch which is closest to the sphere of 
radius K centered on the reciprocal lattice 
point he, For 4.0 -.6<-5. 0 thesre bation sShuge 
between the conventional method and the K-sphere 
method is shown. The relationship for a number 
of other values of 6 can be obtained from 

Table I. eelne, Kosphere me Chod Goes not sao ly 
when a systematic reflection is close to its 
Bragg condition. 
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this thesis as the K-sphere method (see Fig. 23). 

An analysis of Table 1 shows that two primary Bloch 
waves are excited in each case. In the conventional 
Bloch wave labelling scheme, these two primary Bloch 
waves are found to change when a high order systema- 
tic reflection passes through its Bragg condition. 
For example, the primary Bloch waves change from five 
and seven to six and eight as the value of 6 passes 
through 7.0. The change occurs in such a manner, 
however, that the ine primary Bloch waves are always 
labelled o and g in the K-sphere labelling scheme. 
The results obtained in Section 5.2.1 can now be 
explained by noting that the many-beam intensity 
profiles will be two-beam like in character (i.e. not 
complex) when there are only two primary Bloch waves. 
As indicated in Section 3.2.3, the average extinction 
distance for this two-beam like profile will very 
nearly be equal to the reciprocal of the separation 
between the two primary branches of the dispersion 
surface, i.e. the branches associated with the primary 
Bloch waves. It can be seen from Fig. 24 that these 
primary branches approach the spheres about o and g as 
6 increases. Since the reciprocal of the spacing 
between these spheres is the kinematical extinction 


distance, the agreement between the many-beam average 
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. K-Sphere 
Centered 


K-Sphere 
Centered 


Primary 
Bloch Wave 


A part of the many-beam dispersion surface 
for the systematic case. When a high order 
systematic reflection is not close to its 
Bragg condition, there are only two primary 
Bloch waves. These Bloch waves are associated 
with the dispersion surface branches closest 
to the spheres of radius K centered on o and 
g. When a high order systematic reflection 
passes through its Bragg condition, the 
primary Bloch waves in the conventional 
labelling scheme change as shown in the 
diagram. 
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extinct onsaistance and the kinematical extinction 


distance improves with increasing 64. 


Seon Dye Numericalanalysis Of tthe variation in 


extinceLlon Grstance results 


The variation in the many-beam extinction 
distance with thickness can also be explained by 
examining the magnitudes of the excitation amplitudes 
in Table 1. Although there are only two primary Bloch 
Waves, several secondary Bloch waves are also excited, 
It can be seen from this table, for example, that the 
excitation amplitude of Bloch wave one is quite signi- 
ficanteacy = 35.5, and. that. this, excitation amplitude 
becomes progressively smaller with increasing 6. It 
was found that the percentage variation in extinction 
distance an Fig: 21 between values of 6 of about 3.5 to 
5.0 could be primarily attributed to this secondary 
Bloch wave. At higher values of 6, it was: found that 
the variation in'rextinctionsdistance could "be "attributed 
to secondary Bloch waves labelled 2g, -g and 3g in.the 
K-sphere labelling scheme. Of these, Bloch wave 2g 


was usually the most important (see Table l). 
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Se3.l(c) yanalytical examingtion of the primary and 


secondary Bloch waves 


The spetnele just presented successfully 
explains, trom a emer) point of view, how many- 
beam effects arise in the case where no systematic 
reflection is close to its Bragg condition. However, 
further insight into how these many-beam effects 
arise can be obtained by deriving analytical expres- 
sions for the excitation coefficients of the primary 
and secondary Bloch waves. ee exact expressions 
Lor the eigenvalues and .eGigenvectors of a large A 
“matrix cannot be obtained, an approximate approach 
based on the methods of matrix diagonalization out- 
lined in Appendix A will be considered. 

If the procedure given in Appendix A is used 
to diagonalize Lieamany=bealiA gia t his saLnen Sones ie Lgien = 
Values COrrect to: first order are Simply the elements 


Ole the: diagonal Or Ay es. e. 


> (A) 


are we es eae ee =e aes ores 


(In this approximate method for diagonalizing the A 
matrix, the eigenvalues and corresponding eigenvectors 
are labelled according to the K-sphere labelling 
scheme.) The eigenvalues of the A matrix correct to 


second order can be found from equation A.14 to be 
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2 
(ee 
fe es neal: bo (5.2) 


P g#h (2K) *(s,-s,,) 
As indicated in Appendix A, these second order eigen- 
values will be a good approximation if the correction 
terms in the summation of equation 5.2 are small 
compared to Sp: Similarly, the orthogonal matrix of 
elgenvectors C/acorrect to fursteorder, can be cai- 
culated from equations; A.9 and A.i2. This matrix of 
eigenvectors, given in Table 2, will be a good 
por osimati on if the off-diagonal elements are small 


compared to the diagonal elements, i.e. 
|2K(s,-s),) | >> |Ug_y! ; top cae a) 


Since the terms oer are non-zero, the approximate 
expressions of equations 5.2 and Table 2 will break 
doOwie G re approaches Sh (Gl-en). =POr the -systematac 
case, ‘this can only occur if a systematic reflection 
is in its Bragg’ condition (see Fig. 25). The appro- 
ximate eigenvalue and eigenvector expressions of 
equation 5.2 and Table 2 can therefore only be good 
approximations if no reflection is close to. its Bragg 
conda tions. 

It is next of interest to deduce which Bloch 


waves will be the primary Bloch waves. This can be 


done by comparing the magnitudes of the Bloch wave 


VoL 


ee ig 
($2) . —— ee 8 a igi c Mle 
(a= a) CaS tp 
-aepis sehito bioosa saant ~A kibmsqqa i beteotbnr 3A 
noisnsziep oi .ti sottemixvesaqs Boop « sd Iliw seylev 


lfisma orn $.2 soitsope to aoisanmmve oft ai sprise2 


to xittem isnopottzo sad ',yitelimte 12 OF bsssqmos 
pp ad as> ,1asbto saxia oF FoaariSeS ,o 2z0tosvaspis 
to xiztem eid? .Sf.4 bas @.4 eaotsaupe mort bezeios 


boop 6 ed ifiw ,& sidst mi aevip , 2iotDsvaSepis 
{femme sis einemels Isanopsibh-tie edz ti noltsemixorggs 


.s.i ,#inemels Isnopsetb sda ot Beitsqmioo 
i] “4 
(€.2) id & By) “4 a o< [Cyt p22 aS | 


sismixoxrggs sds \oOues=-nom sis Pe a” amtst oft sonke 
Aseid fLiiw % sidsT bone S.¢€ anostaups Yo enoiaesigqxze 
Sttamoetaye sit 20% tad & 4H? 2adoso1rdgs o* he awob 
noisselie: olstamsseaye s if avsdo yiao oso eid? (be69 
-oiggqe ofT . (28 spi Soe) cary ifbon posta a7i al ei 
1G adekeaaunes tososvispis bas sulsvnepisa stemix 

boep od sen Snctenan nap & sidsT bas S.¢ ae 


sean mgacaias on tt peer 


m Ay ae Poe 


- a y > 7 


ih 


108 


f, © Sie whe © © werels 8's 6 tele 6 6 4 2 ie 2 50010 0.6% 5% 0 win 6 6m tie © cieke ec. 0 ote 6 sterekelenchvacns tencno iene aa 


Bek oe aS. a 6 
(ls-’siuz | =) (Pts—Psh ye See) iZ (Sees : Suz) _ = 
=6 Se yee 3 z B = 
| vie Sie coe ste a- 
6 Be 5 5 
Isyz ie °syz Suz Csyz SUZ f 
= ee ae me Be Be 
ea oan a oF, Qa 
Hf FSO OSES WR. (OOO (0 910 6.9) 18 <0 0 Wage 010) 0 0 16 0 16 (eu0 (0 (8:0 12 6 80 4.06 6 eleteie (6 ote omens! ole) ems tetckeronctanseenennne 
5, gran 5 es Pes 5, 5, = 9 
(U) (5¢) (5-) (62) (5) (O) = 
°5 ps meh els O8 oe °5 
(u) (Be)* (b-) (62) (8). ae) 


"O XTajzeU ToJOSeAuehbTs weeq-Auew sy} OF uoTAeUTxOATdde AJepiIo AsATZ ou °Z eTqeL 


.2 xiztem x0t>99vnepis msed-ynsm ont ot aoistsmixozggs 9b10 dext 
! 


at _ 
: ; | a Te are 
os es Pe “ 
Qo o . : o | 4 


(ii) 
p oe" Ss e D 


(il) - 
o” 


Seale ate nt eae S 6 8 WE Oe bE 6 OE Sig O8 EFS Her 


* 2 


a 
Dis 


ee ee ee : 
= Va 


7 


‘i. 


ices te Bs 


Peet. 


Ewald Sphere 


The. conditions under which the s-values of 
two different systematic reflections are 
eGudi ames Vallee spe uOreany a eCrlecu lon 2) 
Viste bemegun la GO Lhe S=Vvaluess py Lorrany 
other systematic reflection h only if the 
reflLeeconon Sh istinvits Bragg condition: 
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excitation amplitudes ee obtained from the values 
of ay and oat given in Table 2. In assessing the 


relative magnitudes of the excitation amplitudes 


ere it is important to note that the numerators 
\ 

of the terms on and oe decrease rapidly in magni- 

tude as |h| increases (see Fig. 16). This has been 


found to be the controlling factor in determining the 


relative Magnitudes of the excitation amplitudes 


c (A) , Ch) 
O g 
‘of high order Bloch waves (i.e. Bloch waves which in 


(see Appendix B). The excitation amplitudes 


the K-sphere labelling scheme have a large |h]) 
therefore tend to be smaller in magnitude than the 
excitation amplitudes of lower order Bloch waves. 
Since the magnitudes of the off-diagonal elements 


of the C matrix in Table 2 are much ‘less than one 


LO) Pam Oy Senne 


, the excitation 
1) g 


(see Appendix B) while C 
amplitudes of Bloch waves o and g will be considerably 
larger in magnitude than the excitation amplitudes of 
the remaining Bloch waves. Bloch waves o and g are 
therefore primary Bloch waves. Since the magnitudes 
of the remaining excitation amplitudes decrease with 
increasing nly the most important secondary Bloch 
waves are the low order Bloch waves such as 2g, -g 


and: sig orgrhis mismiricagrecnenth wich the? results pre- 
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5.3.1(d) Analytical methods for comparing the 


kinematical extinction distance and the 


many-beam average extinction distance 


It can be seen from Fig. 19 that the agreement 
between the kinematical extinction distance and the 
many-beam average. extinction distance improves with 
increasing 6. This result can be explained by 
examining the first and second order eigenvalues of 
equations 5.1 and 5.2. Since there are only two 
primary Bloch waves o and g, the approximate value 
of the many-beam average extinction distance will be 


so) a) | . 


given by 1/|y It can be seen from equation 


5.2 that the magnitude of the second order correction 


‘to the eigenvalue ee will depend on terms of the 
type 2 
hapa) 
peop : (5.4a) 
h#o (2K) Sy 


Similarly, the second order corrections to the eigen- 


value y $9) are of the type 
Ceo 
pag tee Ug =H) c aun sen (5.4b) 
hg ears 
(2K) Se S)) 


On the other hand, if the eigenvalues correct to 
first order (equation 5.1) are used to obtain the 


extinction distance The ge ae eee RN ae then the 
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kinematical extinction distance i is obtained. 
The difference between the many-beam average 
extinction distance and the kinematical extinction 
distance can thererore be seen by examining the 
magnitudes of the correction terms 5.4a and 5.4b. 
In examining the magnitude of these terms, it is 
Pep eond CESS 


in the numerators decrease rapidly in magnitude with 


important to note that the terms (U 


increasing |h|. This has been found to be the 
controlling factor in determining the magnitude of 

the second order correction terms 5.4a and 5.4b. 
Because of this, the higher order terms are much less 
important and therefore can be neglected. The magni- 
tude of the Second order correction terms can therefore 
be examined by considering only terms of low order. 
However, the magnitude of these low order terms 
decreases with increasing 6 because the denominators 
increase in magnitude (see Fig. 26). Accordingly, the 
many-beam average extinction distance and the kinema=- 
tical 6xXtinceionsd1 stancesmust, approach one, another 
with increasing 6, in agreement with the results pre- 
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Fig. 26. A geometrical construction showing Chatee bor 
Barge oc and a Low order rerlection i ene 
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5.3.2 Analysis of the Results for Orientations where 


a High Order Systematic Reflection is Close to 


its Bragg Condition 


5.3.2(a) Numerical approach 


The many-beam calculations of Section 5.2 
showed that the variation in the extinction distance 
with thickness increased for values of 6 close to the 
Bragg condition of a high order systematic reflection. 


In order to see how this effect arises, the excitation 
(ci) 
g 

close to the Bragg ™*conditronot—a-high-erder-systema- 


amplitudes rear have been examined for orientations 
tic 'retiection: 7 The results Por-6§=.4.97, 70 = 550 
and 67)= 5.055 Given, in Table 3, are typical of those 
obtained. It should be noted that the conventional 
Bloch wave labelling scheme is used. It can be en 
from this table that as the value of 6 passes through 
5.0, the magnitudes of the excitation amplitudes of 
Bloch waves three and five decrease while those of four 
and six nepease. wor ¢ = 5.0.45 can jbe seen that all 
four Bloch wave excitation amplitudes are approxi-: 
mately of the same magnitude. 

Because there are four primary Bloch waves 
for 6 ~ 5.0, the diffracted beam intensity cam be 


written to a good approximation as 


6 : 
op. = ) D; eetoniy a) (Ses) 
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thas) 


Typical Bloch wave excitation amplitudes 
for values of 6 close to the Bragg condi- 
tion of a high order systematic reflection. 
The primary Bloch waves are underlined. 
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where D. 
ah oO 


and absorption has been neglected. 
The dispersion surface branches for these four pri- 
mary Bloch waves are shown in Fig. 27. It can be 
seen from this figure that Lor o1closemco 5.07 the 

Res enya ies tor Bloch waves three and four differ by 


(CRE ae) 


an amount equal to y - Since the terms 


7 : 

D3 and D, are of the same sign (see Table 3), these 
two Bloch waves will be in phase at the top surface 

of the crystal. Therefore, the contributions made by 
these Bloch waves to the diffracted beam amplitude will 
Bevout Of qphadse at a thickness given by z = |1/(2c) 7. 
Similar arguments can be used to show that Bloch waves 
Eive and Sixewitl be Out sor phase at a thickness 

2 = | 1/ (2d wheres = pe ee However it can be 
seen fromgFig. w/anehnat d i= c/2. sker”“the purposes of 
the qualitative argument to be presented here, it will 


therefore be assumed that Bloch waves five and six 


are in phase at a thickness |1/ (2c) 


EQUdGLOMn Io so Calm then De WhitLen 


us a p'exp(2niy'?)z) + (De +Dg) exp (20iy \7) 2) (SO) 
where | 

D' = D3 + D,exp (2m1icz) 5 
Lit, ice assumed tiac D4 re es “De = -De (this is 


very nearly the case as can be seen from Table SNe 
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Fig. 27. A plot of a portion of the many-beam 
dispersion surface in the systematic 
case when the reflection 5g is close 
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then the term D' will decrease in Raa aceunen 
increasing z and become very small for thicknesses 
close to |1/(2c)|. For thicknesses close to |1/(2c) | 
the form of the intensity profile will he similar to 
that obtained when one primary Bloch wave (corres- 
ponding to the De+De term) .and.other. secondary Bloch 
waves (corresponding to the D' term and secondary 
Bloch waves such as 2g in the K-sphere labelling 
scheme) are excited. The beating between the primary 
Bloch wave and the secondary Bloch waves leads to a 
large variation in the extinction distance for thick- 
nesses close to.|l/(2c)| as was in facteseen in rig. 18 (b) . 
A complex region was also observed in the diffracted 
beam intensity profiles close to 6 = 6.0 for reasons 
ET rtarato those given above. 

a a es Ommoe Close, toes. 0,, 62.0; and 9.0 
complexity was not observed in the range of crystal 
thickness considered. The reason for this is that at 
these values of 6 the separation c is considerably 
sma lleneathanicatade= po 40s O4eOn= 6 Ove Therefore yethe 
thickness at which complexity would occur is consid- 
erably larger. At these large thicknesses, the 
diffracted beam intensity is very cael and this 


complexity is not significant. 
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59.3.2(b) yAnalytical approach 


It is next of interest to examine the reason 
why more than two primary Bloch waves are excited for 
6 close to the Bragg condition of a high order system- 
atic reflection. This can be seen by examining the 


magnitudes of the Bloch wave excitation amplitudes 


c(h) (h) 
O a 


Table 2. If a high order systematic reflection h 


obtained trom the 117rst order C Matrix 1OL 


approaches its Bragg condition, i.e. Sy approaches 
ZeOveccieerar Galoer SCC eLrOle Cie tact ln malin lable 2 


that the term 


(c(h) sh 
Gas 2Ks) (5.7) 
becomes very large. In addition, the term 
=) 
(h-g) h 
C = AR I (S04 
EN SAGS) ss 
_ g,°h-g! 
also becomes very large because Sy approaches Shea 


(see Fig..25).. According ly,; Bloch; waves; heandah—-g 
Wild, beocomesprimary,.Bloch.waves.¢) Furthermore, he 
excitation amplitudes of all other secondary Bloch 
waves remain small. Clearly, therefore, for values 
of. d6.close, to the, Bragg condition) of+a-.hagh.order 

systematic reflection, the four primary Bloch waves 


GreCuel— dj wande Dawid be excated. 
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It should be noted that the above analysis 
shows: that™ayhighis order systematic reflection h can 
have a significant effect, on. the intensity I. in osa 


. g 
perfect crystal when S, * 0 and when S, > So. 


g 
(Cockayne (1972) has previously used arguments 
Similar to those just presented in order to show 
that, in the systematic case, these conditions should 
be avoided if Pea pani difitiractaon, conditions are to 
be obtained.) In the systematic case, these two con- 
ditions reduce to the one condition that no systematic 
reflection should be close to its Bragg condition. 
However, it will be seen in the next section that 
these two conditions do not reduce to one in the non- 


systematic case, and therefore both conditions are 


LMpOLtCant. 


yx 4 ne tt ect Or NON-ovystematloeRerleéctrons (on Many— 


Beam Effects 


In the many-beam calculations presented in 
Section 5. 2ggonLy systematic reflections were taken 
intor account. ~In order to’ examine the role which non- 
systematic reflections play in many-beam effects, 
calculations including these reflections were also 
carried out. The results of these calculations will 


be presented and discussed in this section. 
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5.4.1 The Results of the Calculations which Included 


Non-Systematic Reflections 


The calculations taking non-systematic reflec- 
tions into account were carried out for the (110) dark 
field intensity in molybdenum. The crystal orientations 
which were considered are indicated in Fig. 28. The 
accelerating voltage used was 100 kV and the fifty 
most important reflections were included using Cann's 
criterion of equation 3.2 (Cann, 1973). An example of 
the variation in the extinction distance which was 
found for the thickness range 0 to 1500 A is shown in 
Fig. 29. The crystal orientations considered in this 
figure were obtained by incrementing jthe)ipoint.t. from 
Gore Oo a2 et RC Om 625, eo], 87) (SCCo PIG. 626s ee Lt 
can be seen from Fig. 29 that there is a variation in 
the extinction distance of at least 7% from all orien- 
tations and that there are many peaks present in the 
curve. As the thickness range used in calculating 
the percentage variation in extinction distance was 
increased, additional peaks were found. These peaks 
indicate that there are many orientations at which 
complexity 2s obtained.’ Inthe cases where no com- 
plexity was found, a many-beam average extinction 
distance was obtained. This average extinction 
distance was found to be within one or two percent 


of the kinematical extinction distance, 
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crystal orientations considered in the 
calculations which took non-systematic 
reflections into account. The projection 
of the tie-point L was initially specified 
to have: Che coordinates (3.3255 5.5257.) 5 
This coordinate was then incremented by 
(OPO OPO LT) Suntrl ‘a 4+hinal coordinate. of 
(25 oe25 f1l2)* was’ obtained. 


ie] 
"a 
~ 


be aor 5 4] 
Pai 
«he ’ 2 
egy CPs. | ta 
BLiwsl ,+° Oi. 
eyange + ‘? 
, ut 
: P | / | 


* 

a | } 
Fat Pe a is 
‘ 


x x ‘ 
a er ee 


. 


/ 
-* 
ee hee 


>. 
7 
AS so 
a 
-6@ . wea : 
ee 


emurer ae 


—— = 


edt ai boisbLenoo. a2noitasneizo Isse . 
pitamateye-mon Aoos doidw eaoiteiuaiso 9% a7) | 


(3) 


Variation in Extinction Distance 


120 


100 


80 


60 


40 


20 


Fag’ 


Mae sc age ay OU ony Ye EPR PV Cy oye Spon Pn Oc eeaawn CPCS ch aha 


29%, 


L 


The percentage variation in extinction 
distance: for the. (10). intensatyan 
molybdenum when non-systematic reflections 
are taken into: account, The crystal 
orientation is indicated by the coordinates 
of the projection of the tie-point (see 
Fa gqnue2 ey yetamcrystals thickness range of 0 
to 1500 A was considered. 
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5.4.2 An Analysis of the Results for the Non- 


Systematic Case 


An analysis of the results of Section 5.4.1 
showed that complex profiles were obtained when more 
than two primary Bloch waves were excited. The 
additional primary Bloch waves were found to be 
excited for crystal orientations where a non-systematic 
reflection had a value of Sy close to" Zero” or “ie The 
reason why additional primary Bloch waves are excited 


at these orientations can be understood by examining 


the expression for the excitation amplitude 


Ue 
th etc gwen e See (5.9) 
OE ae: (2K) S,(Sg7Sp) 


F68 Siech@waveshwoutained= Grom tie: Caematrixwin. fable 2. 
It can be seen from equation 5.9 ree the excitation 
amplitude of Bloch wave h will become large as the 
value of Sy s6ies SeHeS PASE S) KO): e en" AS a result, 
there will be three primary Bloch waves o, g, and h 
at these orientations and this results in a complex 
intensity profile. 

In examining the reasons for the complexity 
if the®*antensity* profiles) 1t was founda that ithe 
effect of a non-systematic reflection h decreased as 
its order increased. The reason for this can be seen 


from equation 5.9. This equation shows that, fora 
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particular value of s the denominators will be 


h?’ 
constant and the numerators will be larger for the 


lower order Bloch wave. Accordingly, the excitation 


(h) , (h) 
g 


é will be larger in magnitude for 


amplitude C 
lower order reflections. 

The result that, when 2 is large, many-beam 
ebfects@ane i important i fon os Over =¢ is in agreement 
with the results reported by Cann and Sheinin (1974A, 
1974B). These Atuhdis found that when the first order 
be flection cg is close to its Bragg condition, the dark 
field intensity profile g becomes complex when a non- 
systematic reflection has an s-value elosesin one of 
the eigenvalues of a two-beam calculation including 


the reflections o and g. It can be seen from the 


expression for the two-beam eigenvalues given in 


equation 3.5 that, for large 6 this condition reduces 


to. thescondrtion given) in) thevprevious paragraph. 
The fact that a non-systematic. reflection 
Can Dekmivortant) Lor values of Sh CLOSEG) LO) 2ero sok. 


So is of significance in experimental observations 
of lattice defects under weak-beam diffraction condi- 
tions. . Under, these, dittracti10on. conditions sass 


desirable in most instances to avoid many-beam effects 


125 


(COCK aye ys LO) 2) iy ee has been assumed in the literature 


that the effect of a non-systematic reflection can be 


avoided if these reflections are far from their Bragg 
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conditions. The results presented in this thesis 
show that the effects of a non-systematic reflection 
can also be large if 345 So: Therefore, this condi- 
tion should also be avoided if weak-beam diffraction 
conditions are to be obtained. Since the value of Sq 
is large under weak-beam conditions, this means that 
non-systematic reflections can have a marked effect 
even though they are far from their Bragg conditions. 

The results discussed a hee also show that 
in theoretical calculations of weak-beam contrast it 
is necessary to take all reflections into account 
which have Sh values close to zero or Js Since this 
Lactels not reflected in Cann's strong-beam criterion 
of equation 3.2, an improved criterion needs to be 


developed. 


5.25 The BHitect of Atomic Number, Accelerating Voltage, 


and the Order of the Systematic Row on Many-Beam 


Effects 


In the last three sections, the SApeanenee of 
many-beam effects on the Seon parameter 6 was 
examined. These many-beam effects were found to 
depend on the magnitudes of the excitation amplitudes 


of the primary Bloch waves 


(eGo eGo) | 
oN Cg Sami = ae (5.10) 
ee. 
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(In order to simplify the analysis to be carried out 
in this section, the t-values of equation 4.8 miles 
Usedstor Che dragonal elenenee ODP tien nat el <a Wien 
the eigenvector matrix C is calculated correct to 
first order, the primary and secondary excitation 
amplitudes are given by equations 5.10 and 5.11. In 


these equations, Th 


of the tie-point onto the systematic row to the reflec- 


tion h. See Fig. 14.) The importance of many-beam 
effects was also found to depend on the magnitude of 
the second order corrections for the many-beam eigen- 
values y 6? and yD), the magnitude or etiese correc= 
tion terms determines the difference between the many- 
beam average Extinction distance, correct to second 
order, and the kinematical extinction distance which 
is obtained a the first order approximations for 


ek and y oe Using the t-method, the second order 


Cobrection terms 5.4) are of the form 
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is the distance from the projection 
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It will be shown in this section that the magnitudes 
of the excitation amplitudes of equations 5.10 and 

5.11 and the magnitudes of the eigenvalue correction 
terms 5.12 depend on atomic number, the accelerating 


voltage and the order of the systematic row. 


5.5.1 The Effect of Atomic Number ‘on Many-Beam Effects 


In this section, the effect of atomic ies. 

on the Bloch wave excitation amplitudes of equations 
5, 10:.and 15.4: Weand son stheset-genva lie rcorre chien berms 
5.12 will be considered. As materials of higher atomic 
number are considered, the scattering factors used to 
calculate the Fourier coefficients UL in equation 4.4 
will increase (see Hirsch et al., 1965). If two 
materials of the same crystal structure and nearly 

equal lattice parameters are compared, values of UL, 

of the higher atomic number material will tend to be 
larger than the corresponding values for the lower 
atomic number material (see Fig. 16). Since the lattice 
parameters are Similar for the two cases, the corres- 
ponding t-values will be nearly equal. For a particular 
accelerating voltage, this means that the magnitude of 
the excitation amplitudes of the secondary Bloch waves 
(equation 5.10) and those of the primary Bloch waves 


(equation 5.11) will both increase with atomic number. 
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However, if the terms ee and econ Opeequationmon er 
aad @cthessinglesterm Dog in equation. 5.10 all increase 
by approximately the same factor, then the secondary 
Bloch waves ght become relatively more important. 
This should lead to an increased variation in the 
extinction distance in the higher atomic number 
material. Furthermore, the correction terms 5.12 
for the many-beam eigenvalues will also be larger 
for the Nhigher*®atomre «number wmmaterial.« This should 
lead therefore to a greater difference between the 
Many-beam average extinction distance and the kine- 
matical extinction distance for higher atomic number 
materials. 

In order to see if these effects do in fact 
occur, calculations were carried out for the high and 
low atomic number materials gold and aluminum. These 
Matervals are) both face centered cubiciin structure 
and have lattice parameters which are very nearly 
equal. The accelerating voltage used an this) calcu 
lation was 00 kv and sixteen reflectionsrvot thew 11) 
systematic: row were considered. The percentage 
Variation in the extinction distance for the (111) 
darkiseleld! intensieyawas (Obtaineayfrorscryscals 
twenky) extanctienydistancessthick. (ihesresults are 
Show Lilelapler4e bor vo =bo..5, col=s/.5, and: Go = 8, 55 
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Table 4. A comparison of many~beam effects in 
materials of high and low atomic number. 
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eonetaeranie increase in the variation of the extinc- 
Clon si usyetevaver: with dpntalat tue ders the higher atomic 
number material. The percentage difference between 
the Many-beam average extinction distance and the 
kinematical extinction distance was also found to be 
larger for higher atomic number materials as shown 
inelajles4. = Clearly, these results indicate that many— 
beam effects increase markedly in materials of high 
atomic number and therefore the kinematical theory 


becomes a poorer approximation. 


5.5.2 The Effect of Changes in the Accelerating 


Voltage on Many-Beam Effects 


It fis next of interest to examine how the 
excitation amplitudes of the primary and secondary 
Bloch waves (equations 5.11 and 5.10) and the second 
order corrections to the eigenvalues (terms 5.12) 
change as the accelerating voltage is changed. As 
the accelerating voltage is increased, the magnitude 
of Un, increases Bevaee the relativistic correction 
term 8 used in calculating this term increases (see 
equation 4.15). The numerators of the excitation 
amplitudes of equations 5. 10° and) 5.11 theretore 
increase as the accelerating voltage is increased, 


while the denominators of these terms remain constant. 
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(it) should, be.noted that.if.s-valués. are used,for the 
diagonal of the A matrix, it is not clear how the 
denominators ARS, Gore ry of the excitation amplitudes 
change with accelerating voltage since both the s- 
values and K change. The t-values rather than the s- 
values Were: Used an. equations.o.. 10, oid and: 5.22 aa 
order to overcome this problems)....Since all. the;:terms 
U_y increase by the same factor 8, the eta ree 
amplitudes of the secondary Bloch waves will increase 


with respect to those of the primary Bloch waves o and 
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G by a Tactomis. pAccordinghy .theevarniatLonean extinc- 


tion distance ie increase with an increase in the 
accelerating ene 

The effect of increases wane the accelerating 
voltage on the many-beam average extinction distance 
can be seen by examining the first and second order 
approximations for the eigenvalues of the two primary 
Bloch waves. When the t-method is used to. calculate 
fog eee aren. approximation to the many-beam 
eigenvalues, the ph else yas Seno the eigenvalues 
oan Pre two primary Bloch waves o.and g is found to be 
Itomte |Z 2K). (This is simply the difference, between 
the elements t, and te on the diagonal of the A matrix 
and is equal to the reciprocal of the kinematical. 
extinction distance.) As the accelerating voltage 


increases,the quantity |rg-tgl/ (2k) changes by a factor 
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of 1/K. The second order corrections to the eigen- 
Gabies sor the primary Bloch waves o and g also change 
by a factor 1/K. Baweret in addition, the numerators 
of these terms increase because the terms U_h 
increase. Therefore, the correction terms become 
more significant for higher acceleratren voltages. 


Accordingly, the percentage difference between the 


-Many-beam average extinction distance and the kine- 


-.Matical extinction distance should increase as the 


Pebaie racine voltage increases. 

In order to see if these effects do in Foe 
occur, calculations were carried out for the (EEO) 
dark eee intensity in molybdenum. The results for 
accelerating voltages of 100 kv and 1000 kV are shown 
tnp Lable.5., Bee enue’ lege one Be ie (110) system- 
oe row were included and the calculations Were, 
carried out for 6 = 6.5, 7.5 and 8.5. The variation ° 
in extinction distance and the many-beam average 
extinction distance were obtained form crystals 
een yec wast ion distances thick... .It can» be,seen- 
from Table 5 that there is a considerable increase 
in the variation of the extinction distance as ne 
accelerating voltage is increased. It can also be 
seen from Table 5 that the percentage difference 
between the many-beam average extinction distance 


and the kinematical extinction distance increases 
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with accelerating voltage. Clearly then, these 
. Peetits show that many-beam effects also increase 
markedly with accelerating voles and as a result 
the kinematical theory becomes a poorer approximation. 
In order to examine the role of the relativis~ 
tic correction term 8, ea licdlabions were also carried 
out Sipe ieee no relativistic mass correction was” 
applied to the terms Uae (i.e.)°8 = none hve 
results garevsiown Aan =lablé; 6" #itecan be Seéneirom 
this table that there is good agreement Metyeen rhe 
results obtained using the different accelerating 
voltages. This result is of interest since the 
increase in many-beam effects with voltage has been 
studied by many authors (see for example Goringe et al, 
1966; Howie, 1966, 1970; Fukuhara, 1966; EnaeBO unos 
Pee Sheinin, 1974). These results have been attribu- 
ted CCOntiewrelacivistiCeiIncrease in the U's as well 
as an increase in the radius of the Ewald sphere. 
However, the results presented above show that the 
increase in many-beam effects which occurs with | 
accelerating voltage is purely a relativistic effect. 
The fact that the radius of the Ewald sphere does not 
have an influence on many-beam effects is a general 
- result first pointed out by Uyeda (1968) and can be 
seen by examining the A matrix when the t-values are 


used for the diagonal elements. 
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Table 6. 


The effect of changes in the accelerating 
voltage on many-beam effects in molybdenum 
when the terms U, are not relativistically 


corrected. 
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® variation in 
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5.5.3 The Effect of the Order of the 


Systematic Row on Many~Beam Effects 


It is lastly of interest to examine how many=- 
beam effects depend on the order of the syeeenatic row. 
If a lower order systematic row is considered (Fox 
example, a (110) row instead of a (200) row), then 
the numerators in the expressions for the excitation 
amplitudes of equations 5.10 and 5.11 and in the eecond 
order eigenvalue correction terms 5.12 will increase 
(see Fig. ro Furthermore, the denominators in these 
terms will tend to decrease in magnitude since the 
values of Th decrease. Accordingly, the excitation 
amplitudes of the secondary Bloch waves should increase 
with respect to the excitation amplitudes of the 
primary Bloch waves 0 and g. Therefore, there should 
be a larger variation in the extinction distance for 
the lower order systematic row. In addition, there 
should be a larger difference between the kinematical 
and many-beam average extinction distances, The 
reason for this lies in the fact that the queena 
order gs, eee correction terms 5.12 will be larger 
for the lower order row because the numerators increase 
and the denominators decrease. Ghetwe other hand, the 
“difference [no - tg 1/ (28) between the primary Bloch 
(0) 


wave eigenvalues y and y 9 correct. to. first order 


will-decrease for the lower order row. Accordingly, 
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the percentage difference between the many-beam 
average extinction distance and the kinematical 
extinction distance should be larger for the lower — 
order systematic row. 

In order to see if these effects do in fact 
occur, calculations were carried out: for»thes(110) 
and (200) systematic rows in molybdenum. The 
accelerating voltage used was 100 kV and sixteen. 
reflections were included. The percentage variation 
in the extinction distance for the (110) and (200) 
dark fields is shown in Table 7. The crystal thickness 
considered was twenty extinction distances and the 
calculations were *carried out for 6 = 52k 8 DAP As 
and 6 = ‘8.5. it can be seen from this table that there 
is a considerable decrease in the Heese of the 
extinction distance with thickness for 7 chew highex, 
order systematic row. The percentage difference between 
the many-beam average extinction distance and the 
kinematical extinction distance was also found to 
decrease. These results show that many-beam effects 
increase markedly when a lower order systematic row is 
considered and therefore the enero theory is a 


poorer approximation. 
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Table 7..-A comparison of many-beam effects. for high 
and low order systematic rows. 


% difference 
between the 
many-beam aver- 
age extinction 
distance and 
$ variation in the kinematical 
Systema- extinction distance at extinction 
LiceROws oa distance jat. 
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5.6 Summary and Conclusions 


inven. eee the validity of the terenecical 
theory for large values of 6 has been assessed by 
comparing it to the two-beam theory. In this thesis, 
the validity has been assessed by comparing the kine- 
matical theory to the more exact many-beam theory. 
Two methods were used to carry out this comparison. 
The first method is similar to the method pees. Cie iy 
used to assess eié validity of the kinematical theory 
and involves a comparison of the kinematical extinc- 
tion distance to the many-beam average extinction 
distance. The results of this comparison showed that, 
if no systematic reflection was close to its Bragg 
condition, the kinematical and many-beam theories were 
in good agreement. This first method, however, has 
serious limitations because the many-beam ayers 
extinction distance 1S not sensitive to the effects 
atecarseed with the presence of secondary Bloch waves. 
Accordingly, a second method was also used to compare 
the two theories. In this method the variation in 
the many-beam extinction toe. with thickness oe 
examined. This method showed that there were signi- 
ficant variations in the many-beam extinction distance 
even when no systematic reflections were close to the 


Bragg condition. This result shows that the kinematical 
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theory is a poorer approximation than it had pre- 
viously been considered to be. | 

~When these two nono deee re used to compare 
the two theories close to the Bragg condition of a 
high. order systematic reflection, many-beam effects 
were found to become very important, i.e. the intensity 
profiles became complex. The kinematical theory is 
therefore a very poor approximation for Paces dif- 
fraction conditions. The nena cavenc many-beam 
theories were also compared when non=systematic 
reflections were included in the many-beam calculation. 
It was found that complexity could occur in the many- 
beam profiles if the s-value of a non-systematic 
reflection was close in value to Zero or Sq: The 
kinematical theory is therefore a very poor appro- 
SAMacioOn at etniesemaLrtrraction Conditions. | This 128 an 
important result because it shows that a non-systematic 
reflection can have an important effect when it is close 
to its Bragg condition. However, a non-systematic 
reflection:-can also have an important effect when it 
is not in the Bragg condition provided that s, fe So: 
This second condition has not previously been mentioned 
in the literature and is important in experimental 
observations ast lattice defects under weak-beam diftirac- 


tion conditions. 
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The effects of atomic number, accelerating 
voltage, and the order of the systematic row were 
also examined. It was found that the kinematical 
theory became a poorer approximation as atomic number 
or accelerating voltage increased or if a lower order 


systematic row was considered. 
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CHAPTER 6 


EXACT AND APPROXIMATE METHODS FOR TAKING 
ABSORPTION INTO ACCOUNT IN THE DYNAMICAL 


THEORY OF ELECTRON DIFFRACTION 


652" Introduction 


The effects of absorption are taken into 
account in the dynamical theory of electron diffrac- 
tion by the addition of an imaginary part to the 
lattice potential (for a recent review see Dederichs 
(1972) or Kambe and Moliére (1970)). The approach 
widely adopted in the literature to the solution of 
the Schrodinger wave equation including this potential 
is to assume that the imaginary part is much smaller 
than, the, real part, thus permitting first order per-= 
turbation theory for the non-degenerate state to be 
used. There are, however, several limitations to this 
approach. Firstly, this form of perturbation theory 
is only applicable to diffraction conditions for which 
Bloch wave degeneracies do not occur. (In this thesis 
the term Bloch wave degeneracy will be taken to mean 
an equality in the kinetic energies of the Bloch waves.) 
During the last several years, however, there has been 
considerable interest in phenomena in which the occur- 


rence of a Bloch wave degeneracy is involved. For 
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example, the critical voltage effect (Lally et al.,1972), 
which results in the minimization of the diffracted 

beam intensity of a second order systematic reflection, 
is associated with the Hen ete) of Bloch waves two 

and three. Also, when non-systematic reflections are 
taken into account in the dynamical theory, Bloch wave 
degeneracies are found to occur under certain diffrac- 
tion conditions (Gjgnnes and Hgier, 1971; Cann and 
Sheinin, 1972, 1974A, 1974B). 

A second limitation of the standard theory is 
that, it is not applicable to diffraction conditions 
which result-in (uasSa-deqcnérate or Close lever states. 
Close level states (i.e. nearly degenerate Bloch waves 
or close dispersion surfaces) occur for a wide variety 
Sbeavuercaction cond tcione. For example, they occur in 
the aot where systematic reflections are taken into 
account and when a reflection is areas COnOtntie tes 
Bragg condition. Some of these levels become particu- 
larly close for certain critical accelerating voltages. 
Close levels also occur at certain orientations in 
calculations which include non-systematic reflections. 

A third limitation of the standard method for 
calculating absorption coefficients is that it As a 
first order approximation based on the #8 ctiiptt oh that 
the imaginary part of the lattice potential is much 


smaller than the real part. Sheinin and Andrew (1974) 
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in a recent paper have shown that this assumption, 
although valid for the lighter elements, may not be 
valid for heavier elements because the imaginary part 
of the potential increases with increasing atomic 
number. 

The limitations of the standard method clearly 
indicate that the development of alternative methods 
ae taking absorption into account is desirable. 
Approximate methods, which enable calculations 0 
be carried out under conditions of exact degeneracy, 
have been developed by Sheinin and Cann (1973) and 
S5acene and Gevers (1973B). Serneels and Gevers have 
also developed expressions based on second and higher 
order perturbation theory for the non-degenerate state 
which can be employed when the imaginary Datt ole the 
lattice potential is no longer a small perturbation. 
As yet, however, no approximate methods have been 
developed for handling the case of quasi-degenerate or 
close level states. 

An alternative to the approximate methods 
discussed above is to perform calculations which make 
no assumptions about the size of the perturbation or 
about how close the kinetic energies of the Bloch waves 
are to one another. One method, referred to in this 
thesis as the exact method, has been used by Thomas 


(1972) and Baily et al. (1972) in calculations of 
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critical voltage phenomena. This method, based on 

the Bloch wave concept, simply diagonalizes the 

Complexe tat rasa AA) which occurs an the: eigen— 
value formulation of the equations of the dynamical 
theory (see equation 6.3). Other methods nde for- 
mulated on the Bloch wave concept are also available, 
as for example, the multi-slice method of Cowley and 
Moodie (see Goodman and Moodie, 1974) and the Howie 

and Whelan (1961) differential equation approach. 

All these theories are useful from a computational 
point of view, but they do not give any insight into 
the interactions between Bloch waves which Peer: from 
the introduction of the complex lattice potential. In 
this chapter an alternative approach which permits these 
interactions to be studied will be presented. The 
theory,.is based ion generalized perturbation theory 

‘and diagonalizes, without approximation, the complex 
Matrix (A+ 1Al) of thesexact theory. « Because this 
approach is a perturbation apEroaGhy the solution to 
the equations of the dynamical theory including absorp- 
tion depend explicitly on the mixing of the unperturbed 
Bloch waves. This new formulation of the theory has 
the following advantages over the methods AF eceaneiy 
available: 

(a) It offers insight into the changes in the 


mechanisms of image contrast which result from the 
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introduction of the complex lattice potential. These 
changes involve the mixing of Bloch waves. 

(b) The relationship between the lower order 
approximations, based on perturbation techniques, and 
the exact theory can be seen simply and directly. 

(c) An approximate method, based on perturbation 


theory and capable of handling quasi-degenerate states 


can be developed for calculating diffracted beam inten- 


eas, This method, referred to in this thesis as. 
the sub-matrix method, gives considerable saving akg) 
computing Pine wihle, ate tie. sale Gite,  OLVing  Vety 
close agreement with those theories which do not make 
perturbation approximations. 

in Sections) 6.2 and)6.3 of this chapter, the 
theory will be developed in its exact form. This will 
be followed in Section 6.4 by a discussion of some 
properties of the parameters Seree aivhen represent the 
mixing between pairs of Bloch waves. The reduction of 
the exact theory to approximate methods for taking 


absorption into account will then be considered in 


the remaining sections. 


6.2 The Solution of the Equations of the Dynamical 


Theory Including Absorption Using Generalized 


Perturbation Theory 


As shown in Chapter 2, the high energy electrons 
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in the crystal can be described in the dynamical theory 


of electron diffraction by Bloch waves of the form 


ip = exp (2mik ‘2) .F) Jus 
= Sanne : r) yc 1) exp (arias r) » (6.1) 
g 


ine vPOuGrLer sCOSCLi ac 1 ents a and the Bloch wave: vectors 


Rt) can be found by solving the eigenvalue equation 


AC a= tat je (6.2) 


Wiechw Ene El rects OoLredadbsorpt1on- are, Laken, into account 
by introducing an imaginary part to the lattice poten- 


tial, the eigenvalue equation 6.2 becomes 


iain) xe ' ae ’ (n) 


(A+4A')C (623) 


where iA' is a matrix in which the elements iach are 


proportional to the Fourier coefficients of the imaginary 


part of the lattice potential iu'(r) through the relation 


LA = 1U5_p/2K- The: matrix (A + iA") is no longer 


Hermitian and thus has complex eigenvalues and eigen- 


vectors. The real eigenvalues and eigenvectors yin? and 


ct) in equation 6.2 have therefore been rewritten in 


! (n) 1 (n) . 


equation 6.3 as the complex quantities y and C 


respectively. 
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The approach to be taken in this thesis to the 
solution of equation 6.3 makes no assumptions regard- 
ing the Si dedar the perturbation iA’ or the spacing of 
the levels and is based on generalized perturbation 
theory. In this approach, the eigenvectors of a 
perturbed equation are written in terms of a linear 
combination of the eigenvectors of an unperturbed 
equation (see for example Raimes, 1961). Accordingly, 
it is assumed in solving equation 6.3, referred to in 
this thesis as the perturbed equation, that the N 
eigenvectors and the corresponding eigenvalues of the 
unperturbed equation 6.2 are known. once the eigen- 
values of this unperturbed equation are discrete and 
(i) 


Enewelgonveclons. C form a complete orthonormal set, 


1 (n) 


the eigenvector elements oF 


can be written exactly as the linear combination 


3) “pind SI) 


gi) ey amc) OT ada) 
1 . 


6E?in matrix .form 


Gon eaga (6. 4b) 
in eivi Ss Vequatlons.Geised matrix in which the columns are 
(n) 


the eigenvectors of equation 6.2. anda is a column 


vector containing the a coerficient ihe the linear 
: c aie lal : q (n) 
combination as the 1 Cleneicy sl Beseigenvecrorec 


in equation 6.4b can then be substituted into equation 
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of the perturbed equation 
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This expression can then be multiplied by oe , the 
Miversenor C,)topgi ve 


ag 


Lcjtac. $c pene a hex ae . (6.6) 


Equation 6.6 can be expressed in a more convenient form 


by defining the parameters 


piie)) cee arene Ju'juG= (I 2K)\PS, Bee chews 
Cpe 


g g-h 
and a ty) = ent) , (Om 7) 
where a is the absorption coefficient of the Bloch 


wave as normally defined in the literature (see equation 
2.29). The use of equations 6.7 and 6.2 then allows 


equation 6.6 to be rewritten in matrix form as follows 


(iement(:) ard. 25 ig (13) 
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(6.8) 
This equation can then be abbreviated to the eigenvalue 


equation 
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Tis can be seensirom equations 6.7 ana 6.8 that the 
perturbed eigenvalue equation 6.9 can be solved only 
by first solving the unperturbed equation 6.2, 

Since no assumptions have been made about how 
close the levels are to one another or about the size 
of the perturbation, the eigenvalues ule Of. the 
perturbed equation 6.3 can be found without approxima- 
Gion by, diagonalizing the complex matrix B in equation 
6490 The corresponding eigenvectors, anus of this 


matrix can then be substituted into equation €.4b to 


' (n) 


construct without approximation the eigenvectors C 


Obeequation G23 hight sCerszornthogonalyand anny iS enoGs 


] (n) 


malized, then C will also be normalized. 


The form of the wave functions including absorp- 


tion can be found from equation 6.1 by substituting the 


' (n) 


complex eigenvector elements Cy 


(a) 


for the unperturbed 


> (n) 
and the complex wave vectors k 


a Ta) 


Z 


elements C for 


ame where k Gy ™ aK), The perturbed wave 


functions, eo can then be written in the form 


[pt (), = exp (2nik! ™ .Z) ) cf" exp (2nig-) - (6,10) 
g 


Clearly these wave functions are also Bloch waves since 
equation 6.10 has the form of a plane wave modulated by 


a function ao (Ms =m) of) exp (2mig-F) which has the 


fel 


{@,9) 


eld tedd¢ 8.3 bas \.3 enotdsaups mori? nose 8d “nso 31 
yiao beyloe ed aad @.3 nottsupe eulsvnspis bediuszeq 
8.98 solssups bedavsusgqny sit privioce sexi yd 

wod juwods sham ased svsi enoistqmueas on sonte 
ssia ett tuods to rorsonE ene ot 9xs elovel sdt seolfs 


(qa) $. 
1 


efit to aoulavnaspis eft ,»noistsdautiregq edt to 


~smixorgas tuodtiw hive? ed mao £.9 noltsups bediusieq 
noitsups ot @ xixésm xsiqmoo $03 paisilsnopstb yd nord 
2isia Yo _ ae 5B ,atotvevnspis paiibnoqes110>2 eit .2@.3 

ot dé.d soitsops otat Setusiszeadua sd asi? asd xitism 


ft) ¥, , F . 
— 9 2uetgosvnepie edt aoitsmixoiggs juodsiw Joursenos 


“ton ai ia bas {saoportso 2i D2 21 .t.a noitsups to 


-besilsmyon sd oals Iliw inde non> .Sositem 
-qxroeds phibuloni enotgonut evaw siz to mioi edt 


aiid paitutisedue yd {.% dotdswps moxt bavot ed aso noth 


(rt) oh 


hsdivtyeqnay aft 262 as adasmels 1032S9vaspis xelqgmoo 


(11) ws 7 (i) 
x02 A exodosv evew xeiqaoo edt bis 3 atmaneite 
evew bedawdasq sAT hy A + (a) + = (A 94 A oxen... MOG 


mso3 odd nt ‘agdsiaw od nods ee ral ae 
7 bed io 


A =i 
- ae 


a a 7 


ie, ity 
yc a= <<  Ud 
das 4: 
A : - i 
= a 
' mn 


ted da, . thgg be 


Penlodveityeomeacnes ati ce. Tt Astor Mnterest to 


1 (n) 


express the Bloch functions lu > in terms of the 
parameters characteristic of the unperturbed Bloch 
waves. This can be done by substituting equation 6.4a 


into equation 6.10 to obtain 
ieee exp (2mik' ™ .Z) Jur ™ 


= exp (2nik! (™ .Z) ) ciset Gave : (62113) 


This equation indicates that the perturbed Bloch func- 


jut ™), are in fact linear combinations of the 


tLONns 


unperturbed Bloch functions |u‘*’>, 


This can be inter- 
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preted as indicating that the perturbation (the imaginary 


part of the lattice potential) results ina He ng OL 
the Bloch a cher ecterat cS of the system without 
absorption. It can further be seen fron equation 
Chee clla te Lie off-diagonal elements, a ey of the 8B 
matrix represent the mixing which occurs between a 
particular pair of unperturbed Sige nines! Loand a7; 
since these elements depend explicitly on the Bloch 
functions |u‘*) > and ju J), The appearance of these 


q its) elements allows the importance of the mixing of 


Bloch waves to be assessed. In Section 6.4, some 


relevant properties of these elements will be discussed. 


The importance of this mixing on the diffracted beam 


intensity will then be considered in Chapter 7. 
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6.3 ‘Calculation of Diffracted Beam Intensities 


The calculation of diffracted beam intensities 
in the dynamical theory involves a determination of the 
extent to which the various Bloch waves in the crystal 
are excited. In the standard method, this is done 
relatively simply since the Brendes, Ofethenc 
matrix leads to the result that the excitation co- 
efficient of a particular Bloch wave, eee is equal 
to Crag This result, however, iS not valid in the 
general case since the matrix of eigenvectors is no 
longer orthogonal. The purpose of ae section is to 
outline the) procedure involved in calculating diffracted 
beam intensities in this case. | | 

The method for calculating diffracted beam 
intensities from the Bloch wave functions, Ip (MS 
of the perturbed system is analogous to the method used 
in the standard Bloch wave formulation of the dynamical 
theory (see Section 2.1.2). The starting point for the 
Cavculatilon, is the total wave function for the electrons 
in the crystal 
v=J x) jor Mey x) Foot) exp (ami (Et) 4G) 0X) 

aati | ie ae | 


C6..12) 


(n) 


where X is the excitation coefficient for the wave 


i ee) 


function > of equation 6.10. If equation. 6.12 


is then Ported Ito incident and dritracted beams and 
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the resulting terms are multiplied by the phase term 
+ 
exp (-20iK-r), the expression for the diffracted beam 


intensity is found to be 


Alth A 
yh =a[ XC! exp (eri! 2) (", (6.13) 
) : 


(n) 


The eigenvalues y' in this equation can be found by 


diagonalizing the matrix B, and the eigenvector elements 


on! (n) 
a 


Sie eee beiaeiboes INU sequal Tonle lie exc1 ca 


can be found by substituting the eigenvectors of 
tion coefficients can be determined from the boundary 
Conditions, atythe top surface of the crystal, ive. 


Oo lay > = 0 (g # 0) oe) 


where ae is the amplitude of the directly transmitted 
beam and o, are the amplitudes of the diffracted beams. 


Equation 6.13 together with 6.14 gives 
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where C' is a matrix having the elements Cue in the 
th ee) , 
g row and n column; X-1s a Column vector contain— 
(n) th 


ing the excitation coefficients X°© " in the n~” row, 


and u is a column. vector containing the diffracted 
beam amplitudes 1, 0, 0,;>%.. of equation 6.14.)) The 


(n) ; 


excitation coefficients X required in equation 


6.13, can then be found by obtaining numerical solutions 
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to the nonhomogeneous system of linear equations 6.15. 
In contrast to the standard theory, these coefficients 
are, in general, complex. 

Although the excitation coefficients can be 
found from equation 6.15, an alternate form of this 
equation can be obtained by substituting in equation 
6.15 the matrix form of equation 6.4, C' = Ca, where a 
is a matrix which has as its columns the column vectors 


cies Using @cnevorthogonality of the Cemateix and mul- 


Gipily ing Boy Cae from the left then gives 


akg ao (6.16) 
where Gis a column vector containing, the elements oe 
in the jae row. The excitation coefficients required 


in equation Patey ee then be found by solving the non- 
homogeneous system of linear equations 6.16. This 
alternate method for calculating the excitation co- 
efficients leads to convenient simplifications in the 
approximate diffracted beam intensity expressions to be 
derived in: Section 6.6 of this chapter. 
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6.4 Some Properties of the gq Elements of the B 


Matrix 
Before examining the approximations which can 
be made in diagonalizing the B matrix (see equation 
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some of the properties of the off-diagonal elements of 


this matrix. The first important property is that for 


centrosymmetric crystals (td) = qd mer chess 


matrix is symmetric. This can be seen from the defini- 


tion of eee in equation 6.7, Section 6.2, together 


with the fact that for centrosymmetric crystals a'=(a')? 


UierL rans poseilo is wie rand ink = om These relationships 
give 
RCAC CIA) CeACC) =O AC) = Gu a 6.17) 


where q is a matrix formed by the imaginary elements of 


ise 


Another important property of the g matrix is 


that, under certain diffraction conditions, some of the 


off-diagonal elements are identically zero. For 
example, if a centrosymmetric crystal is oriented so 
that the systematic reflection mg is in its exact Bragg 
Gond beELOn and only a systematic row of reflections is 
considered, then the aS enveseer’ of the unperturbed 


equation 6.2 are either symmetric, with 


Chih) =k 
5 Se ae Pat (6.18a) 


or anti-symmetric with 
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MuNGLCIUMe= nO pnd eae form msevenmand in =41,3,5 5.4. for m 
odd (Metherell and Fisher, 1969). If eigenvector 
elements ary and oe! are substituted into equation 
6.7, it can be shown by using equations 6.18a,b that, 
for a centrosymmetric crystal, the ee ene q't3) is 
identically equal to zero, provided that one of the 
Bloch waves is symmetric and the other is anti-symmetric. 
A PEreLCuraTmeccnpis of this second property has pre- 
viously been noted by Serneels and Gevers (1974), who 
showed that the matrix elements describing the mixing 
between Bloch waves two and three are identically zero 
in a critical voltage calculation carried exits precisely 
at—29-in* the*Bragg® condi-tions 

(Gah OBIE ge property of the q (tJ) elements can be. 
obtained by comparing these elements for values of the 
deviation parameter 6 which are equal in magnitude but 
opposite in sign (see Fig. 30). For the systematic row 


Of neflections shown in this figure, the A matrix for 


6 = y/g can be written 
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Ewald Sphere Ewald Sphere 


SAGE 


The Ewald sphere construction for deviation 
parameters 6 = y/g and 6 = -y/g which are 
equal in magnitude but opposite in sign. 

The points L(y) and L(-y) are the perpendi- 
cular projections of the tie-points T(y) and 
T(-y) .onto the systematic row. ~ 


soitsiveb 10% aoitoustened ererge biswa sit 06 SQte 


brs doirtw p\y- = 6 bas p\y = 4 exetemsisg 
apis ni etia dud ebydinpam at Late Ne ; 


~ | y-)d bes (yy) etatog siT 
bane (y)T TE, 16 anoisosfoxrg talv> Oe 


dsmezeye ord ono (y~)T fi. ioe 


Loo 


The elements on the diagonal of this matrix are the 
t-value elements of equation 4.8, i.e. 


(ty (y))* 


t,(y) = - —43— (6.20) 


where th ly) is the distance from the point ty yto the 
reciprocalalattice pointeh (see Pag. 30). (For simpli- 
city, the t-method approximation for calculating the 
diagonal elements Of thetsjmatrax willbe used ine this 
section. This will not affect the arguments presented, 
however, Since these arguments will be based on the 
SsIgciveclouseOLalne 7 amalld k,mandn Wiley elgenvecvors 


found using the t-method are exact. See Section 4.3.4.) 
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Since t(y) is the distance from the point L(y) to the 
reciprocal lattice polnteahyseltacan be seen from 


equation 6.20 and) Fig. s0nchat 


 (aduren shephhavly ee (6.22) 
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Furthermore, for centrosymmetric crystals, a = 


Uaers 


= 8) 


Accordingly the A matrices of equations 6.19 and 6.21 


are identical. The real, normalized eigenvectors of 


these matrices will therefore be equal, except for an 


auberrarvamulcuplacabive, factor ofitl,. 1.4. 
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Therefore, from equation 6.23 and the definition of . 
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The hoes elements are therefore equal in magnitude for 


values of the deviation parameter 6 which are equal in 


magnitude but opposite in sign. With the proper choice 


Of the arbitrary Signain equationi.o.25,.1unesG 
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elements for these values of 6 can be identically equal. 
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that,. when only systematic reflections .are considered, 
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Iq (43) | is a periodic function of the deviation para- 
meter 6. This property can. be derived by comparing 
the A matrix at 6 = y/g_ to 

LhesArmatrixvato'= (y+ 20)/9g 


where $ is a reciprocal lattice vector (see Fig. 31). 
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Since Ty fy) is thevdistance strom the goointe Ly) cto 
the reciprocal lattice point h, it can be seen from 
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Eicgt ol... ine Ewaldespneresconstruct10n,,forrdeviation 
parameters 6 = y/g and 6 = (y+ 2)/g where 
is a reciprocal lattice vector. The points 


Liy) and-Ly +°2)-are “the*perpendicular pro} 
yections ot the tie=points Ti(y) and Tiy +22) 
onto the systematic row. 
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From, equation 6.2/7 and the. definition .of q (13) in 


equation 6.7 it can then be seen that 
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The magnitude of the elements Rey is thus a periodic 
EUNGeLON VOL On Wi Lora 5 RO aie toe Z ir If f= g 
where |g| is the distance from o to g in Fig. 31, then 
the period is 2.0 in units of 6. Metherell and Fisher 
(1969) have previously shown this result for the diag- 
onal elements ae) = oh OL SCheS Be Matrix. s. Liieiee 
argument is based on Bloch's pheoren and can be exten- 
ded to prove the result of equation 6.28. 
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magnitudes of the terms gq vary with 6, thirteen 


beam calculations for the (220) systematic row in gold 
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have been carried out. The results are shown in 
Pig. 3208 Cangbe seen from this figure that the 
values of }q (43) | are symmetrical about 6 = 0.0 and 
that they are a periodic function, of 6 with a period 
equal CO ze UR RUE Chermore,TOr Oo "=an Where ’n Us.cn 
CiL3 ) 


integer, some of the values of q are zero because 


Bloch waves of opposite symmetry are involved. 


6.5 The Relationship Between the Exact Approach and 
the Low Order Perturbation Theory Approaches 


in eaction 6,.2,, thelLexact solution to the 
equations of the dynamical theory of electron difitrac— 
tion including absorption was formulated in 6erme Oba a 
matrix B whose off-diagonal elements, fey ae were seen 
to represent the mixing which occurs between pairs of 
Bloch waves. In this section the assumptions concern- 
Lag, Che importance of Bloch wave mixing which must be 
made in order to reduce the exact theory to the low 
order perturbation theories previously presented in the 
literature will be considered. The mathematical pro- 


cedure which will be used is discussed in detail tin 


Appendix A. 
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6.5.1 The First Order Approximation to; the Solution 


of the Perturbed Equation 


The low order perturbation theory approaches 
discussed in Section 6.1 can be developed from the 
exact theory by calculating approximate expressions 
for the eigenvalues and .cigenvectors oOfsthersamatri xs 
in order to calculate theseigenvaluess,of.thees ama trim 
correct to first order, it is assumed that different 
Bloch waves do not mix and, therefore, that all the 
off-didgonalyelements of thesBamatrixjareszexno.. qeunder 
this first order japproximatilonauhe sBematarix then: takes 


the diagonal form 


Pi roel) O O ehbcere 
Bx 0 eo ec O BS (6229) 
O O SSE Taig ee ORLA ory 


See a) 


The perturbed eigenvalues y of this approximate B 


matrix are then equal to the standard theory eigen- 


values ay i (n) 


theory for the non-degenerate state. The corresponding 


eigenvectors aoe of the approximate B matrix can be 
written as columns of the unit matrix av (see equation 
A.5 of Appendix A). Reference to equation 6.4b then 
shows that fhe perturbed eigenvector matrix C' can be 


approximated by the unperturbed eigenvector matrix C 


q found using first order perturbation 
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and that the Bloch wave functions ptm) 


; PeNen oecome 
those used in the standard first order perturbation 


approach, namely 
pt Ms & Gain) : (6.30) 


ee therefore can be seen that if the assumption is fade 
that the imaginary part of the lattice poten: a does 

not result in mixing between any of the Bloch waves, 

the exact solution reduces to the incomplete -first | 
order perturbation theory solution for the non-degenerate 
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6.5.2 The Second Order Approximation to the Solution 


of the Perturbed Equation 


In: the epee order approximate solution to the 
“matrix B, the effect or Bloch eye mixing is taken into 
account in an approximate manner by assuming that the 
mixing between two particular ayaa se LOAIiCl ge Wael. 
only affect Becee eng siles and eigenvectors of Bloch 
tied i and j.- Within this second order approximation, 
then, the only off-diagonal elements Gr tne matrix B 
which are taken into account in the calculation of the 
corrected eigenvalue and yarns of ja particular 
Bloch wave i are ig le? and iq 3) where j=1,N (j# i) 


and N is the number of beams taken into account. 
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In.order to calculate the perturbed eigenvalue 


1 (n) 


Y of “Eheanatrax By CcOLrect, Lo, second, order ,~it.1s 
first convenient to cast the perturbed eigenvalue 


equation, 6.79 in the form 


(iB - {y'})a™ = 0 | (6.31) 


where {y'} is a diagonal matrix containing the diagonal 


elements yn) 


the first order eigenvalues on the diagonal of the 
Matrax Beare non-degenerate, ,and,that.the.second, order 


corrections to the eigenvalues are small, then, to a 


good approximation the diagonal terms on SP Oe, 


(m) ,. sqm) _ Obes Hele 


can. be set equal Oe Y. iq for m #n. 


Furthermore, if it is assumed, as mentioned above, that 


only; the. off-diagonal; terms.in, the. B.matrax which 
involve Bloch wave n directly will have an important 
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effect on y then, alli other off-diagonal elements 
can be neglected and can be considered to be zero... 


Under these second order approximations, equation O.ok 
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As shown in Appendix A, the eigenvector a can be 


th 


obtained by solving all but the n equation repre- 


sented by the matrix. equation 6.32 to give 


p (1) my epee Ee (ea) (6533a) 
J ol Cie OO Gy, 
, . 4m) | | 
where Bane = ar and pe cd ha | (6.33b) 
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Since it has been assumed that the first order eigen- 
values are non-degenerate, the denominator in equation 


6.33a will be large and it can therefore be assumed 


that oe << 1) (j#-n).. -Theretore the’column vector 


(n) 


as whoseselements are a mecane Dem Cons ldeked: co: be 


normalized and the eigenvector a can be taken to be 
equal to ee The eigenvector elements of eu! can 
then be found from equation 6.4b to be 
ct) _ op(™) 
(6.34) 
or Gat CP 


: ; th : ae 
where P iS 9a matrix’ whose —n COLUMN CON ta iis sue veigen— 


(n) : 


VeClLOn P It can be seen from equations 6.33 and 
6.34 that the corrections to the first order eigen- 
vectors obtained from second order perturbation theory 


are complex. If, however, it is assumed that 
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then these corrections reduce to the purely imaginary 
first order corrections derived by Wilkens et al. (1973) 
and Serneels and Gevers (1973B) . Thestact that, in the 
present analysis, this first order result is forth- 
coming from second order perturbation theory is con- 
Sistent with the fact that, in perturbation. theory, 

the calculation of the perturbed energies to a given 
Order requires the knowledge of the wave functions to 
the next. lower order. 

It is next of interest to obtain the second 
order corrections to the first order eigenvalues. This 
can be done, as shown in Appendix A, by substituting 
the expressions for ae obtained from equations 6.33 


into the nt equation. Ot Oo 2eLOmclve 


nee 4 (vives (a) 


Y + ig ~ ) Gy, 3:G)) 


J ois ae aap 


It should be noted that this second order approximate 
solution gives both a real and an imaginary correction 
to the first order eigenvalues ee iq(™., This 
differs from the purely real correction derived by 
Serneels and Gevers (1973B) using second order per- 
turbation theory for the non-degenerate state. However, 


if the assumption in equation 6.35 is again made, 6.36 


gives the purely real correction derived by these authors. 
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It is also of interest to consider the calcu- 
lation of the excitation coefficients required to 
calculate diffracted beam intensity in the second 
order theory. in the standard theory this task is 


Simplified since these coefficients can be shown to 


be the terms in the first column of pen the inverse 


of the C matrix. Since C is orthogonal, the excita- 
(2) 


tion coefficients are simply equal to C 


In the general case when the eigenvector matrix C' is 


complex, the excitation coefficients are still equal 
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orthogonal and (cc!) "tg (rye THe Gad. CcuractOn- or “tine 
excitation coefficients in this case therefore requires 


that the first column of (eye 


be obtained by other 
means. It is possible in the second order approxima- 
tion, however, to calculate approximate expressions for 
the excitation coefficients by deriving an approximation 


to (cease If it is assumed that Raed etor sla eai, 
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the elements of Pe. (ihits. Can, De Seon MmOsurcastt yan, 
multiplying P by pl, using the fact that oe = -p{J) 
and examining the approximations required to make this 
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(c')+ = (ep) t* = priggtt zs pic? = (c1)7. (GRSZ) 


(n) 


The excitation coefficient, X , can therefore be 


¢ (n) 


ss for the second order case ina 


approximated by C 


manner similar to the standard theory. . Wilkens et.al. 


(1973) have previously given expressions for (c') 7, 


but in their sécond order theory (eovis was not shown 


to be approximately equal to Ce: 


6.6 The Sub-Matrix Approach 


In the second order approximation of the previous 
section, mixing between different Bloch waves was taken 
into account in an approximate way by assuming that the 
second order corrections to the eigenvalues were small. 


It can be seen from equation 6.36 however, that these 


: Let) ene 
corrections will become large as q' J) increases or as 


(5) Ge) Seti ae) 


the term y peaks | becomes small. Since 


the 5 eal terms tend to increase as the atomic number 


increases (see Section 7.2) and certain of the terms 
9) 


Sa) ie eee) can become small, for example, 


ay He] 
near critical voltage diffraction conditions, it is of 
interest to develop an Satay eee method for diagonali- 
zing the B matrix which can be used when the low order 
pereurDArion soluetens £or the non-degenerate state do 


not apply. Such an approximation, termed the sub- 


matrix approach, will be developed in this section. 
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In the sub-matrix approach it is assumed that 
‘the mixing between some Bloch waves is important and 
Should be taken into account, while the mixing between 
others is not important and can be neglected. This 
approach is similar to that used in quasi-degenerate 
perturbation theory where some of the unperturbed 
levels are assumed to be sufficiently close that mixing 
between these close levels must be considered. As the 
Simplest example of this method, let it be assumed that 
the mixing between only two of the Bloch waves is 
important. (The procedure can easily be generalized 

to take into account the mixing between more than two 
Billochs waves. )multy. for example, itevsthe mixing between 
Bloch waves two and three that is considered to be 


important “then the B matrix takes the approximate form 
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The perturbed eigenvalues pee (neF*2,3)2 Of ‘this 


matrix are then found to be those of the standard 


1 (2) (3) 


theory, while the perturbed eigenvalues y and y' 


can be found by diagonalizing the sub-matrix 
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The G@2i.genvectors a of this sub-matrix can then be 
incorporated into the matrix a which has. as columns 
the eigenvectors Ors Che sBonacrie in, equationv6, 35. 


This produces an a matrix of the form 


il O eC Oper 
(2) (3) 
O ay a5 OORr. 
= (2) (3) 
az |O 3 a, Omar « (6.40) 
O O O LA sexes 


The perturbed Bloch wave functions in the sub-matrix 
approach can now be obtained from equation 6.10. For 
those Bloch eats which are not considered to mix, 
the Bloch wave functions can be seen to be those of 
the standard theory, while for those which are con- 


SLdered to mix, in thistcase 2yand 3, 


[ot ™) > = exp (amiien Fy fal Jal?) > 4 al fy G)o3, (n= 2,3) 


(6.41) 


Thus Bloch waves 2 and 3 in the perturbed state are 
formed from a linear combination of unperturbed Bloch 


Functions. 


(6.39). 
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In order to calculate the diffracted igyayerr ARR 
Sity using the sub-matrix approach, it is first necessary 
to calculate the Bloch wave excitation coefficients. 
The general method for calculating these coefficients 
has been given in Section 6.3 where it was shown that 
x (0) 


, tilenexcitativon coefficient = of Block wave n;. can 


be found from the expression 


ax = G@ (6.42). 
li this -equation, fa 45 a matrix whose ee column is the 


(n) 


cugenvector a 


(n) 


Om rahiGwi Machine eroerd OO Nunn VeCCor 


in the are row, and G is a column 


A in the nth row. ~“Tne’a matrix 


containing X 
vector containing C 


of equation 6.40 can therefore be substituted into 6.42 


to give 
dt O O (G ea eesues x 1) c (1) 
fe) 
O ee a Oo xse,s x (4) c 4) 
2 2 Oo 
Ee) (3) (3) = est) 
O a3 a Oiseetars Xx = Ce . (6.43) 
O O O oe rons x 4) c (4) 
fe) 
, é (Te eect) 
It can be seen by inspection that X = Cs LOU ee eso, 


and that x (2) and x3) can be found from the reduced set 


of simultaneous equations 
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(2) (3) {2y) (2) 
ao a5 x ie 
= . (6.44) 
ai?) 4{3)| [x (3) o (3) 


The expression for diffracted beam intensity, for the 


present example, can therefore be written as 


2 _ (n). + (n) ee) 
sl = I = 
" Io, ue x Cy sep l2nty z) 
Mee iG WiC! exoloniy in 2) i (62 45) 
ne Dew! an9, 
where se (n # 2,3) are the eigenvalues, vate AS ae 


(n) 


of the standard theory and y' (ni= 2275) are, the 


eigenvalues of the sub-matrix in expression 6.39. 


6.7 A Criterion for Determining which Bloch Wave Mixing 
Terms are’ Important’ in ‘the Sub-Matrix Approach 


in order Co be able to most conveniently apply 
the sub-matrix approach, it is clearly desirable to 
have some criterion for determining which Pi cetera 
mixing terms should be taken intoj;account. The method 
to be presented in this thesis is based on the fact that 


the contribution of a particular Bloch wave to the 


i diffracted beam is determined by 


intensity of the gt 
EA ea eel 


the excitation amplitude Xo 


In the standard 


theory where Bloch wave mixing is neglected, this ampli- 
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general, result in a change in this excitation amplitude 
which can be determined to first order, by equations 
6.33 and 6.34. It can be seen from these equations that 
the mixing of Bloch wave i with. other BLOC HWAv es 5 
leads to a first order correction to the eee component 
of eigenvector i given by 

it) coe NE Tt eee (6.46) 

g g i#4 3) g . 
Since it was also shown in Section TEE tat, COrrect 


TO ltrst Oder, the excitation coerrigient x Fh) OL 


Bloch wave i is oo the first order excitation ampli- 
tude becomes 
xy er G) = ct (4) ar G) 
g oO g 
| j#i i j#i 4 
 eeithig salen ee ee 
jAi Peed ti a 
(6.47) 

where terms of second order in Ba have been neglected. 


(It has been found numerically that, when the correc- 
tions to all the excitation amplitudes epee iel (1= 1,N) 
are Fs Ws Rec all the significant corrections arise 
from terms of first order in poe Liebe ds auld Goek 
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6.33a then, to a first approximation, the mixing between 
Bloch waves i and j leads to an imaginary correction 


to the excitation amplitude of Bloch wave i of 


ne ._(i9) : 
Abb 2 ve ae es ele el) Gly (6.48) 
ce () (i) O g fe) g 


The mixing between the same two Bloch waves will also 
lead toa Beer none ae to. the excitation amplitude 
of Bloch wave j of ek An, indication, of the impor— 
tance of mixing between Bloch waves i and j can now be 
found by comparing the magnitude of Bae to the magni- 
mee of the excitation amplitudes of the most important 
Bloch waves present. In order to facilitate this com- 
Dans on, it has been found convenient to express 


We | AUS eect percentage of the magnitude of the largest 


anim) 
g 


A term. The use of these excitation amplitude 


correction factors in determining the Bloch wave mixing 
‘terms to be taken into account in a particular sub-. 


Matis -CaLlciwlatwone wil!) be 1 Liustrated an Section /,52 


6.8 Summary 


In this chapter, generalized perturbation theory 
has been used to derive an exact Solution tostne matrix 
form of the equations of the dynamical theory of elec- 
tron diffraction including absorption. It has been shown 


that the introduction of the imaginary potential leads 
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tO a Pane of the unperturbed Bloch waves character- 
istic of the system without absorption. The first and 
second ,order ,approximate solutions to the generalized 
theory were then considered and the relationship 
between these approaches and the perturbation 
approaches previously presented in the literature 
was discussed. A new perturbation approach, termed 
the sub-matrix approach, which is capable of handling 
non-degenerate, By aecen acne and degenerate states, 
was then developed. | 

In the next chapter, the importance of Bloch 
wave mixing in calculations of diffracted beam inten- 
sity will be examined. The different methods for taking 


this mixing into account will then be .compared. 
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CHAPTER 7 


THE IMPORTANCE OF BLOCH WAVE MIXING 


IN. CALCULATIONS OF DIFFRACTED BEAM INTENSITY 


etl eel Pc eOCNGEALOE 


In the previous chapter, ane dynamical theory 
of electron diffraction including absorption was for- 
mulated in terms of a matrix B whose off-diagonal 
elements were seen to represent the mixing between 
pairs of Bloch waves. In this chapter the importance 
of these elements in calculations of diffracted beam 


intensity will be considered. 


J@2 (hese rtect of Atomic Number on the wimportance of 


Bloch Wave Mixing 


In order to illustrate the importance of Bloch 
wave mixing, two sets of calculations have been carried 
out. «In the first of these, the importanceior fulas 
mixing in materials of differing atomic number has been 
examined. This aor has been adopted Pe. ne 
off-diagonal elements of the B matrix are a function 
ene Uys ahah in turn increase with atomic number 
(Humphreys and Hirsch, 1968). Table 8 shows the first 
four rows and columns of B matrices for the low atomic 


number material aluminum and the high atomic 
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Pableeo weet nee Lrrst our rows. and, columns of typical 
bematrices™ for a* thirteen” beam: calculation 
in gold and aluminum. The values to the 
right indicate the amounts by which the real 
and imaginary parts of the eigenvalues of 
these? Bymatrices*difter trom the elements 
on the diagonal of B before diagonalization. 
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number material gold. In calculating these matrices, 
thirteen reflections of the (220) systematic row 
were included, the deviation parameter was 6 = 1.8, 
andethe tacceleratingnvolitagecwasnLog,kVs i <TheesB 
matrix elements of Table 8 include the te dpicte between 
Blochhiwavesgone to lfourntandraresveryrsimplanitosthe B 
matrices which would be obtained from a calculation 
which only included the four most important reflections 
230 to-440.) “It can be seen) from Table/8 that the off- 
diagonal Bhelementsi,aréexmore thaneantactonxoftten 
larger for gold than for aluminum. The Bloch wave 
mixing introduced by the imaginary part of the lattice 
potential is thus more significant for high atomic number 
materials than for low atomic number materials. 
EtGisnnextseofmintenestrytonconsidenathelertfrects 
of the Bloch wave mixing described above on the eigen- 
values of the matrix B. This has been done by simply 
calculating the differences between the exact and 
standard theory eigenvalues. In carrying out these 
Slichasetbaen ae is important to note that the amount 
by which the standard theory eigenvalues (vy aig Sty , 
which are equal to the diagonal elements Se the matrix 
B, are modified by the inclusion of the off-diagonal 
elements depends on the process of matrix diagonaliza- 


tion. In this process, the importance of the off- 


diagonal elements depends both on the values of these 


cel a. 


,esottiam seeds pniteluolse al SBicb ‘kakeoaae Suey " 
wox oidameteye (OSS) eft to emoidoelies nsetzidd 
a,{ = 3 asw rstdompxsq aoltsiveb eft .hbebulont sx1ew 
ad sfiT) VA OOL saw apsifov pnitsisisoos ead Bas 
neswt od enttorkm: afd ebuloani 8 sidsT to estnsmelse xiazem 
4 of3 o¢ asilole yisv sis bas wot oF sao sevsw noold 
moi¢sltioisas s moxzi beaissdo sd bisow Asinw eeotissem 
ahontos lier tassiogmi tzom 2403 sax bebylont yinao doidw 
‘Sas eit ted? 8 oldeT mox2 nese sd moo 41 ¢.088 of OSS 
nat to sotos? 6s csi? stem sts aetnemele € Isnopstb 
evaew doola od? .mwnimels x6o2 neds blop 20ot 1z9pxel 
esoitgei sat to st1#q yasnipemi eit yd beoubortal pakxtm 
1edmun oimots dpid 202 tosoitiopia stom andid ei Isisnetog 
.eistusisa rwsdaun oimots wol xot nsadt elsiissem 

ajositis oli tehienos ot geoxedni Io gxon at 31 


-aspis eff mo evods bediicasb pnixim svsew doold eds to. 


yigmie yd sfiob ased asd aidT .a xtitsm edt to esutev 


a 
a> 


Bas josxs oft meswied esonere?tib ed paiisivoias 


pase 
J 


gaedd tuo poiytiso al .esutsvaepis yxrosdy Bisbaste 

drvoms ot fect eJon ot Jnsitogmi at 42 .anoktetunlias » 

: vets 2 eRe yacedg ree ee 
ana: 


sales a: — Ae ex 


183 


elements as well as on the values and differences 
between.the diagonal elements. In order to assess 
the importance of the off-diagonal elements in 
materials of different atomic number, it is therefore 
important to compare B matrices in fen Che Corres-— 
ponding diagonal elements are as nearly equal as 
possible. Reference to Table 8 shows that the Bo 
matrices for gold and aluminum meet this requirement. 
The results shown in the right column of Table 8 


‘indicate that the difference between the sty Vos eigen- 
values yom) and the eigenvalues eicinths: fe) of the 
standard theory are considerably greater in gold than 
slyg Suite This indicates that Bloch wave mixing 
leads to larger changes in the Sonali Lows 
materials of high atomic number than for materials of 
low atomic number. 

In order to see the effect which Bloch wave 
mixing has on the diffracted beam intensity, the 
matrices of Table 8 and the methods outlinedin Section 
6.3 have been used to calculate the bright field, 
Eurst, OLder, ana Second oLfder intensities in: gold, and 
aluminum for particular values of crystal thickness. 
Each intensity was then compared to the corresponding 
“intensity found using the standard theory. The com- 


parison involved expressing the differences between 


the intensities found by the two methods as a 
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percentage of the intensity found from the standard 
first order theory. The average amounts by which the 
two theories differed for thicknesses between O A and 
3000 A is shown in Table 9. It can be seen Gud, tor 
each of the diffracted beam intensities considered, | 
the differences are larger for the high atomic number 
material. For the diffraction conditions considered, 
the differences for the first order dark field are 
larger than those for thesbright fieldjwhile.the second 
order dark field Bi frerences are the largest of all. 

As an illustration of the differences just 
discussed, the second order dark field intensity 
profulesaan Hee obtained using ;the standard and exact 
theories are shown in Fig. 33. These profiles are 
typical of all those obtained in comparing the standard 
and exact theories in that there is Careers 
Beincerance in the overall shape of the profile or in 
the peak to peak spacing. There is, however, signifi- 
Cant difference in the diffracted beam intensities 
found using the standard and exact theories. For the 
second order dark field profile in Gold showm angfig.s 3, 
this difference amounts to an average of 15%. 

In summary, the results of this section indi- 
date that Bloch wave mixing can have an important 
effect on diffracted beam intensity and that this 


effect is greater in high atomic number materials. 
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Table 9. The average percentage difference in the 
intensities for gold and aluminum found 
using the standard and exact theories. 
The diffraction conditions used are indi- 
cated in the text. 
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200 400 600 800 1000 1200 


Crystal Thickness (A) 


The exact theory (upper curve) and standard 
theory (lower curve) second order dark 
field antbensity profiles; for.godld) ates — 178. 
The diffraction conditions used were the 
same as those for Tables 8 and 9. 
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7.3 The Effect of Quasi-Degenerate Diffraction 
Conditions on the Importance of Bloch Wave 
Mixing 

7.3.1 Anomalous Absorption Effects in Diffracted 


Beam Intensity Profiles 


Re indacated in’ Section 722, the effect of 
the off-diagonal elements of. the B matrix depends on 
the differences between the diagonals elements. it 1s 
therefore of interest to investigate the effects of 
the off-diagonal elements near the critical voltage 
where, for the second order reflection near its Bragg 
condition, two of the Bloch states have nearly equal 
kinetic energies. (For a more complete deseription 
of the critical voltage effect, see Chapter 8.) | 
Under these quasi-degenerate conditions, two of the 
diagonal, clements Of the Bo matrix are nearly equal, 
In this second set of calculations, the importance 
of Bloch wave mixing near the critical voltage in the 
medium atomic number material molybdenum was investi- 
gated. The standard and exact theory profiles for the 
second order dark field at the critical voltage of 
TA KV are shown in Fig. 34) in this) caiculation, 
thirteen reflections of the (110) systematic row were 
considered and the deviation parameter was 6 = 2.004. 
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Fig. 34.: The second order dark field intensity for the 
(110) systematic row in molybdenum at the 
critical voltage of 74 kV and for 6 = 2.004. 
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which Bloch waves two and three are degenerate and 
at which these Bloch waves make important contribu- 
tions to the (220) dark field intensity (Lally et al., 
UO 2) ope Caluabe seen, from rad. 34, that. there as no 
appreciable change in the overall shape of the inten- 
Sity profiles found using the two theories and that 
wane peak to peak spacing is not noticeably altered. 
There are, however, large differences in the diffracted 
beam intensities found using the two theories, and 
these, differences increase. with, thickness. . The 
Per erencas vary from a few percent in thin parts of 
ther crystals tol00s at 1200 Ns Although some 
Giirerences were. also found for Cher bright. f2etd vand 
brs t Order intensity profiles, thes largest: ditferences 
between the two theories were for the Becond order 
dark field. 

it)  vssot interest to note an Fig..34 that Bloch 
wave mixing Risa such a way that the effects of 
anomalous absorption are enhanced (see Section 2.2.3). 
Cilae should be noted, however, that this effect was © 
found to arise because of changes in the eigenvectors 
and not from changes in the absorption coefficients 
ae a Accordingly, the contrast Cala Lin)’ (4 ae +I in)? 
in thick parts of the crystal is less for the exact» 
Eneory tian for the Standard  tnéory.. s'n1s,point, is of 


interest since Lally et al. (1972). have shown that the 
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reduction in the second order dark field intensity 
which occurs near the critical voltage lowers the 
background intensity and leads to improved visibility 
@esdefects.» The, results) of Fig. 34 snow, that if this 
reduction insthe background intensity (i.e. the 
improved contrast) is to be properly assessed in 
intensity calculations, then Bloch wave mixing should 


be taken into consideration. 


7.3.2 The Effect of Bloch Wave Mixing on Rocking Curves 


In order to determine the range of 6 over which 
the differences between the standard and exact theories 
are large, rocking curve calculations have been carried 
out in molybdenum at the critical voltage. (A rocking 
curve iS simply a plot of the diffracted beam intensity, 
at some specified crystal thickness, as a function of 
the deviation parameter 6. See Hashimoto et al., 1962 
for example.) In these calculations, thirteen reflec- 
tions of the (110) systematic row were again considered, 
the accelerating voltage was 74 kV (the critical 
voltage) and the deviation parameter was varied f yarn 
6/-= 2.0 to 6'= 251. The’ second) order dark field’ rocking 
curves for the crystal thickness of 3000 A are shown 
Lier Los. #3 (al and 35 (D)s. . ott .Can be .Seen L.Com 
Pig. o5(b) (that Ghere 1S “dood agreement between the 


standard and exact theories precisely at 6= 2.0. This 
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molybdenum close to 6 = 2.0 obtain 
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systematic, reflections,o£)the,(110 
were unCcluded any the calculation, 
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good agreement resulies fT romestnemtactit hat Bloch waves 
two and three are of opposite symmetry precisely at 
6 = 2.0 and therefore do not mix (see Section 6.4). 
For 6 >2.0 to approximately 6 = 2.02 there is a large 
difference between the intensities found using the 
two theories. For 6 > 2.02 the differences are still 
apparent but become smaller. Clear the results of 
Hig 3:5 indicate that Bloch wave mixing is fd geonye 
atthe critical oirage over a range of values of 6 


ellose ttosoe= 2.0rand\should ber taken anto-.account in 


rocking curve calculations. ~ This findingjyis important - 


Since standard theory rocking curves have been used by 
ality eteal S197 2) to a llustrate she ,eduction in, the 
second order dark field intensity which occurs at the 
ClutlcalesvOltage. 

Rocking curves were also obtained at the 
critical voltage but for values of 6 close to the 
symmetry position (6 = 0.0). Under these ditlraction 
conditions, Bloch waves two and three are again quasi- 
-degenerate. The standard and exact theory rocking 
curves obtained for the bright field at a thickness 
of 3000 Afhare plotted in Pig? 36 Site can be Vseen=rrom 
this. figure that, “for "small nonzero values of Go; the 
intensity calculated using the  Gxact theory 12s consi- 


derably less than the intensity obtained using the 
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Fig...36.+)Bright £1eld,rocking curnvessefortmolybdenum 
close to the symmetry position (6 = 0.0) 
obtained using the standard and exact theories. 
Thirteen reflections of the (110) systematic 
row were considered, the accelerating voltage 
was 74 kV (the critical voltage) and the 
crystal thickness was 3000 A. 
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Standard theory. (This is the reverse of the situation 
Close to 6 = 2.0 in Fig. 35 where the standard theory 
intensity was less than that of the exact Eeeiens 
However, good agreement was obtained precisely at 
mn), = Oy, O reconeey for thus: orientation, Bioch waves two 
and three are of opposite symmetry and therefore do 
not mix. Significant differences were also found 
between the standard and exact theories in rocking 
curves obtained from the first and second order dark 
field intensities. 

In order to examine whether it is necessary 
to be Say aS TOae to che (critical voltage. in order 4to 
obtain large differences between the standard and 
exact theories, calculations were also carried out 
ativoltages somewhat sremoved#trom the critical voltage. 
ie 2 2 0) rocking curve .for molybdenum at a thickness 
of 1000 A is plotted in Pig. 397 foro close to. 2.0 
and for an accelerating voltage of 100 kV. It can be 
seen from this figure that significant differences 
between the standard and exact theories still occur. 
As would be expected, these differences were found to 
decrease as ematis, ees were carried out for voltages 
further from the critical voltage. However, "good, agréee= 
ment between the two theories far from the critical 
voltage could only be obtained in low atomic number 


materials. 
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molybdenum close to 6 = 2.0 obtained using 
the standard and exact theories. Thirteen 
systematic reflections of the 110 row were 
considered, the accelerating voltage was 

100 kV, and the crystal thickness was 1000 ie 
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In conclusion, it can be seen from the results 
preecnte et this section that Bloch wave mixing can 
have a very significant effect on the diffracted beam 
intensity when calculations are carried out under 


quasi-degenerate diffraction conditions. 


7.4 Discussion 


The calculations described in Section 7.3 have 
shown that the importance of Bloch wave mixing is deter- 
mined to a very considerable extent by how close certain 
SrecoesOlagOlaime lements Or thes Brmatri xs dare tOTone 
another. Since the values of the diagonal elements 
for a given material are in turn determined by the 
crystal orientation and the accelerating voltage, the 
importance of including the off-diagonal elements of 
the B matrix in a given calculation also depends on 
these factors. In attempting to assess the importance 
of Bloch woe Mixing at 1s) Chnererore Caper eere LO ecu 
electron microscopist to note that there is a wide 
Variety of diffraction conditions, in addition to those 
at the critical voltage mentioned above, which give 
rise to relatively close values of the diagonal elements 
of the matrix B. A familiar example is the case where 
a high order systematic reflection een the Bragg con- 


dition. Under these circumstances, some of the Bloch 
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waves which contribute significantly to the intensity 
of certain of the diffracted beams are associated with 
branches of the dispersion surface which very nearly 
touch (i.e., the Bloch waves have nearly equal eigen- 
(1), . is) 


values y Since the absorption coefficients q 


are normally much less than the eigenvalues oy LA) this 


would imply that the values of the diagonal elements 
of Cee Bie tint ye? \e associated with these 
Bloch waves, would be nearly equal. 

Roenee example, not so familiar but perhaps 
one important, is the case where non-systematic 
reftlecurvons tare, taken into account. 91t Nas been shown 
by several authors (Gj¢gnnes and Hgier, mene Cann and 
Sheinin, 972, 19 74a, 1974B; Sheinin and Cann, 1373) 
that at orientations close to a non-systematic reflec- 
tron in, the Bragg condition, Bloch wave degeneracies 
or touching branches of the dispersion surface occur. 
Again, this means that, for the Bloch waves concerned, 
equal or nearly equal values of the real parts, GE uthe 


diagonal elements of the B matrix will be obtained. 


Sheinin and Cann (1973) have shown that at such degen- 


eracies, effects of Bloch wave mixing are considerable. 


The calculations in Section /23.2 for, the quasi— 
degenerate case near the critical voltage would indi- 
cate that the Bloch wave mixing which would occur near 


the degeneracy of the non-systematic case would also 
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give rise to Significant effects. -It is clear from 
these considerations, therefore, that there are many 
possible diffraction conditions for which the effects 


of Bloch wave mixing should be taken into account. 


7.5 A Comparison of Approximate and Exact Methods for 


Taking Absorption into Account 


In this section the standard, second order, 
sub-matrix and exact methods of taking absorption into 
account will be compared. The comparison has been 
Carried out by calculating diffracted beam intensities 
int gold Lor diffraction,conditions where no Bloch wave 
dencgetadies were involved, and also in molybdenum 
close to a Bloch wave degeneracy. 

In the gold calculations, thirteen reflections 
Br the (220) systematic row were considered, the accele- 
rating voltage was 100 kV and the deviation parameter 
was 6 = 1.8. | (These are the same diffraction conditions 
as were used for the gold calculations an Section w/, 2.) 
In order to determine the size of the sub-matrix to be 
used, the excitation amplitude correction factors 
derived in Section, 6./ were calculated s=) [hes results 
obtained tor bright=-field, first order dark tieia and 
second order dark field are shown in Table 10. It is 


clear from this table that, for these diffraction 
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conditions, the effect of Bloch wave mixing on diffrac-. 
ted beam intensity should increase with increasing 
eroer of the reflection considered. This is in good 
agreement with the differences found between the 
standard and exact theories for his enientation (see 
the gold results in Table 9). Table 10 also shows 

that in second order dark field calculations, for 
example, only mixing between Bloch waves 1, 2 and 3 
need be taken into account. Accordingly, in calcula- 
tions of second order dark field intensity, the sub- 
matrix used was analogous to the one shown in expression 
6.39, Section 6.6, except that values of eae i=1,3 
and j =1,3, were included. The results obtained from 
the calculations of diffracted beam intensity showed 
that the second order theory was in better agreement 
with the exact theory than was the standard theory. 

The best agreement, however, was found by uSing the 
sub-matrix approach. For example, at a thickness of 
450 A, it was found that the exact theory differed 

from the standard, second order, and sub-matrix theories 
by 153, 8.5% and 0.24% respectively. Although there is 
very good agreement between the sub-matrix and exact 
theories, this agreement could be improved by expanding 
the size of the sub-matrix. In the limit as the sub- 
matrix size becomes equal to that of the B matrix, eige 


sub-matrix theory becomes equivalent to the exact theory. 
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In comparing the different theories it is also 
important to compare the computing times which each 
theory requires. From a numerical point of view, there 
is no point in using an approximate theory if the exact 
theory is just as fast. The total computing time which 
each of the theories required in order to calculate 
Ene wei genvalwes -andweigenvectOrs Or the yA + 1A matrix 
and the eeeation coefficients was therefore obtained. 
Por the gold calculations. just considered, it was found 
thabte the: sub-matrix approach was’ the fastest andythat 
the times taken by the second order and the exact 
theories were larger by factors of two and three res- 
pectively. Clearly, the sub-matrix approach can lead 
to a significant saving of computing time while. still 
giving results in excellent agreement with the exact 
theory. | 

In order to compare the different theories 
close to a Bloch wave degeneracy, the standard, second 
order, sub-matrix and exact theories were used to 
Galcullate second order dark field rocking curves in 
molybdenum at the critical voltage. In these calcu- 
lations, the only large excitation amplitude correction 
factor was found to result from mixing between Bloch 
waves two and three. Accordingly, a sub-matrix was 
used which took only mixing between these Bloch waves 


into account. The rocking curves which were obtained 
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are shown in Fig. 38. In this calculation, thirteen 
mcerbecttonsserethe (110) (syectematic row were consi- 
dered, the second order dark field intensity was 
calculated at a thickness of 3000 A, and the accel- 
erating voltage was 74 kV (the critical voltage). 

It can be seen from Fig. 38 that the sub-matrix 
approach is a very good approximation for all values 
of 6 while the second order theory breaks down close 
to 6 = 2.0, i.e. close to a Bloch wave degeneracy. 
As already discussed in Section 7.3, the standard 
theory iS a poor approximation over this range of 


orientations except for 6 = 2.0. An examination of 


the computing times required by the different theories 


showed that the sub-matrix approach was much faster 


than either the second order or exact theories. These 


results obtained close to a Bloch wave degeneracy 


again illustrate that the sub-matrix approach can 


give a considerable saving in computing time and still 


give very close agreement with the exact theory. 
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Exact and 
Sub-Matrix 


Standard 


Second order dark field rocking curves 


for molybdenum close to 6 = 2.0 obtained 


using the standard, second order, sub- 
matrix, and exact theories. Thirteen 
systematic reflections of the (110) row 
were included, the accelerating voltage 
was 74 kV (the critical voltage) and the 
crystal thickness was 3000 A. 
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CHAPTER 8 


THE IMPORTANCE OF BLOCH WAVE MIXING IN. THE 


CRITICAL VOLTAGE EFFECT 


o..b- intreduction 


Im .the-Last chapter, it jwas found;that the 
effects of Bloch wave mixing are important when high 
atomic number materials are considered or when the 
crystal is oriented so that quasi-degenerate diffrac- 
ELon. Condi tyonsvare .obtained. ):Livis mextpok interest 

EG examine a particular phenomenon in which this Bloch 
eos mixing is important and must be taken into account 
in order to obtain a physical understanding of the 
Con crast Teche iens involved. The, critical ,voltage 
effect, which has received considerable attention in 
the literature in the past several years, is sucha 
phenomenon and will therefore be considered in this 


chapter. 


Sie NSN RewLew OfesthNeuC i141 cal Voltage HEfect 


Before examining the role of Bloch wave mixing 
ins tthe: critical voltage effect 442s? first of Seiterest 
to briefly review the physical nature of this effect. 
The critical voltage effect was discovered experimen- 


tally by Nagata and Fukuhara (1967) and Uyeda (1968). 
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TWese= avenors = found toate ror o =—250 and "a> certain 
incident electron energy, the intensity of a second 
order reflection became very weak. This reduction 
in intensity was also observed by Watanabe and Uyeda 
(1968) in the vanishing of second order Kikuchi lines. 
Following Metherell and Fisher’ (1969) and Lally et al. 
(1972) this effect has come to be known as the critical 
voltage effect. 

They criticaly voltage’ erfect NasPgenerated a 
Great deal or interest tm; thet literature because’ 1+ 
has been found that the voltage at which the minimiza- 
Clo inesecond, Order intensity occurs. “15 e%™= the critical 
VoleiGejeis Very —selsiliuve: tO walless@h, atomic scatter— 
PAGMbactObLay steLruccuLe Jac UOlorerai ds tie luattt ce 
parameter. The effect has therefore been used in the 
accurate determination of atomic scattering factors, 
Debye temperatures, the degree of ordering of partially 
ordered alloys, and lattice parameter changes (see for 
exeamole, Watanabe et alt, 19603" Lally etal. lo 72; 
Buble eels, LIS se OMaAS ey cel ty elo eee © Lene ee Gee 
ras also been. found? tovbe useful ineancreasing the 
Vieibrilutcy Of defect contras (ia lily, et i heepi MR AN 4 

There are two approaches which have been used 
in vthe literature to explain) therminimization jn the 
Secondeorder dark field antensity at the critical 


Vol Gagens “lhe giirst approachy 15, due,to Uyeda (1968) 
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and is based on Bethe's second order approximation 
(Bethe, 1928). In this approach, it is assumed that 
two Strong beams, O did” 2g, are Included an the A 
matrix and that the effect of other weak beams such 
as,.d Caneve teahem Into account, by: CabcuLlating a4 

modi fied value pote for: the Fourier coefficient eye 
(See Uyeda (1968) for the expression for this modified 
Fourier goer evaenck Vorla pareicular, valuejor the 


accelerating voltage, My is ound to .go to zero and 


2g 
the A matrix becomes diagonal. Since the eigenvector 
matrix C is then a unit matrix, the excitation amp li- 
(Ay Rt) 


a Cog of the two Bloch waves are zero and no 


tudes C 
diffracted beam is obtained. 

However, this first approach indicates that the 
intensity log becomes zero at the critical voltage 
whereas in fact the intensity decreases to some non- 
Zere value (ialiy jet ali, 21972)., gin order to more 
accurately describe the critical voltage effect, a 
second explanation based on Bloch wave symmetry argu- 
ments was presented by Metherell and Fisher (1969) and 
ally et al. (1972)... Thisssecond approach considers 


the contribution which each Bloch wave makes to the 


diffracted beam intensity 


cf2)6 2) ass and Gy tig MP2) |? ‘ (8.1) 


Fareddexeculyeequal toe2ze0 ci teisofoundithatypeat! the otop 


206 


aocitamixoxgge x9b10 Bbaoase e'sdted no* beasd ai bas 


gerd benw2a6 ei 4k ,dosorgqs aide ot .(8S@L ,sdtea) 
A etd af bebulonit 615 ,pS bas o ,emssd paorsa, ows 
‘oue 2zmsesd Asow refto to 250779 eft tend Bas. xixtam 
6 paiseLuolso vd trvosos odnt asaists ed nao p es 
“we? tneinitieos rsixyot end tot ss st LEV bait bom 
beliiibom eins tot aoLersiqxs edt xo? (B8CL) sbheyU- 592) 
aris to eulsv wseipoisxsq 6 1O% (,dastobtfeoo zreizvo0F 
bus o%es of op oF Bbnyot et est ,epaifov pricsreiIsoos 
1osoevaepis oft soni@ .isnopaib eemoved xirdsm 4 sad 
-ilqms nolistioxe efs ,xiutem sin 6 neds at 9 xirtsm 
ant Das orez s1x6 eevew doold ows ans to eae esbut 
-honistdo eit msed betos.xtitb 

eit tedt getsoibni dosorggs texti eis ,xsvewoH 
epasioy lapisixo edt 36 otes wamooed. ot yitenetat 
~som emog oJ eeaseipeb ytianesnai edt tost nt eseiedw 
etom of tebio oI. . (800i ,.fe to ylisd) eulsv ores 
5 ,toetie eapstlov Isoitixo edd ediztje0b yleds2uoos 
-~upis yrsemmye svew doold ao boesd noiteasiqxs bnovee 
has (@aeL) rede id bas iistedieM yd betneeerq BBW adnom 


pig = bacoas eidT «(ST eA) oe 


gine ‘ Oe 
ae aca aceasta 


207 


Surface of the crystal, Bloch waves one, two and three 
are primary Bloch waves and of these, Bloch wave one 
is the least important. Furthermore, it is found that 
the absorption coefficient of Bloch wave one is much 
larger than those of Bloch waves two and three. The 
diffracted beam intensity expression 8.1 can therefore 
be written 


(1) 


( ay) 
fe) (2g 


exp (2mi (y 1) 4ig'*))z) +a]? (8.2) 
where A represents the contribution of all Bloch waves 
except “ewo and three. The magnitude of A is very; small 
unless the crystal is very thin, in which case Bloch 
wave one can make a Significant contribution. 

The reduction in the intensity which occurs at 
the critical voltage can now be understood by examining 
the excitation amplitudes and the eigenvalues of Bloch 
waves two and three in equation 8.2. The values for 
these parameters at the exact Bragg condition of the 
(222) retlection an. copper close te. the? cri ticale volcage 
are shownninghig. 8959 lt cannbenseen fromabhis dagune 
that as the accelerating voltage is increased from 250 
to 390 kV, the excitation amplitudes of Bloch waves 
two and three remain nearly equal in magnitude and of 
opposite sign, while the eigenvalues of these Bloch waves 
approach Gach other on value. At 320 KV the eigenvalues 


are degenerate and the excitation amplitudes of Bloch 
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Fig. 39. This figure shows (a) the eigenvalues 
and (b) the excitation amplitudes for 
Bloch waves two and three close to the 
critical voltage in copper. Thirteen 
reflections of “the “(11Il) “systematic 
row were taken into account and 6 was 
taken to be 2.0. 
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waves two and three are interchanged. This interchange 
results from a change in the symmetry of Bloch waves 
Newolandathres (for example, Bloch wave two is symmetric 
below 390 kV and becomes anti-symmetric above 390 kV). 
Lally et al. (1972) have defined the voltage at which 
this change in symmetry occurs to be the critical 
voltage. With an increase in accelerating voltage 
beyond the critical voltage, the excitation amplitudes 
of Bloch waves two and three remain nearly equal in 
magnitude and of opposite sign while the difference 
between the eigenvalues increases. 

The reason for the reduction in the diffracted 


beam intensity I at the critical voltage can now be 


2g 
understood by noting that at the critical voltage 
y (2) = y 3) while the excitation amplitudes of these 
Bloch waves are nearly equal in magnitude and are of 
opposite sign. It can be’ seen from equation 8.2, 
therefore, that the contributions made by these two 
Bloch waves will interfere destructively and a minimum 
in the intensity results. An example of this reduction 
in the diffracted beam intensity is shown in Fig. 40. 
In obtaining the results presented in this figure, 
thirteen reflections of the (111) systematic row were 
included, 6 was taken to be exactly 2.0 and the inten- 


SiL.ty wasecaldaulateds for a crys Cal sthickness of 5000 ae 


It can be seen from Fig. 40 that the intensity is very 
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Piguw AO. Tne critical voltage sntensi ey minimum 
for the (222) intensity in copper. 
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much Sreducsdeclose to the ‘critical voltage of 390 kV 
Dut 1S not equal to zero. This result is in better 
agreement with experiment than Bethe's second order 
approximation. (As indicated earlier in this section, 
the approach based on Bethe's approximation shows 

that the intensity becomes zero at the critical voltage.) 
Therefore, Bethe's approximation has been restricted 
LoNGe termining approximate, values ior the critical: 
voltage (Uyeda, 1968; Watanabe et. al., 1968), while 
the second approach based on the many-beam theory has 
been used in the detailed investigations into the 
Critical voltage effect (Metherell and Fisher, 1969; 
ere yee ce cle Le) oe LOMAS yoo) 2eee OU Le lee OL 2 eels Ss 
Thomas et al., 1974; Serneels and Gevers, 1974; Gevers 
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between these degenerate Bloch waves is not important 


because Gites) 


= O (see also Section 6.4). However, 
thetdifiractdon-condi. tions wircmegive rrse co tie 

exact degeneracy cannot in general be obtained 
experimentally because of such factors as beam 
divergence, specimen bending and the lack of precision 
which is inherent in determining the crystal orienta- 
tion. Furthermore, it has been found by the author 

of this thesis that conditions equivalent to the exact 
degeneracy cannot be attained in numerical calculations 
of diffracted beam intensity at the critical voltage. 
This arises because the number of significant figures 
in computer calculations is limited to some finite value. 
It was therefore felt to be of interest to investigate 
the importance of Bloch wave mixing on the critical 
voltage effect at small departures from the exact Bragg 
condition of the second order reflection. The results 


of this investigation willbe’ présented-in this section. 
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materials thére is no intensity minimum, while in others, 
the minimum persists for a range of 6 roughly between 
1.95.5 O« Go 2.05.. AS -an illustration of these two cases, 
the results ae molybdenum and copper obtained from the 
standard, sub-matrix, and exact theories are discussed 

1G thisisection. jin the sub-matirix approach ; only mixing 
between Bloch waves two and three was included. Further- 
more, high atomic number materials were not considered. 
The reason for this was that any observed effects ae 
Bloch wave mixing could be associated with quasi- 
pedenavsre diffraction conditions and not with atomic 
number (see.Section 7.2). Only the results from the 
standard and sub-matrix approaches will be presented 

| eirice the intensities obtained from the latter method 

and from the exact method were in good agreement. Riso, 
euly- Valles or s0.2e2,0 will be. considered, [c.shouid 
be noted, however, that the results to be presented are 


Loentical te those obtained, ior los 2.20. 


Soul (a) An pexamination of, the critical vyaltage intensity 
minimum in molybdenum 


The second order reflection chosen for the first 
set of calculations was the (220) reflection in molyb- 
denum. Thirteen beam systematic calculations were 
carried out for various crystal thicknesses and deviation 


parameters ranging from 6 = 2.01 to 6 = 2.0. Typical 
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results of the variation of the intensity of a second 
order reflection with accelerating voltage for small 
values of the deviation parameter are shown in Fig. 41. 
The intensity is plotted for a thickness of 3000 A 
and a deviation parameter of 6 = 2.001. An examination 
of the curves in Fig. 41 shows that there is a large 
discrepancy between the results of calculations based 
on the standard and sub-matrix theories. The calcula- 
tion based on the standard method for taking absorption 
into account shows a distinct intensity minimum at 74 
kV. This minimum lies in the same voltage range as 
does the critical voltage defined by the symmetry 
interchange of Bloch waves two and three. Clearly, if 
only a calculation based on the standard method was 
performed, there would be a strong temptation to conclude 
that the intensity minimum in Fig. 41 was the minimum 
usually associated with the critical voltage. That 
such a conclusion would be in error can be seen from 
the intensity versus voltage curve obtained from the 
sub-matrix theory. This, curve shows in fact that no 
minimum occurs in the voltage range around 74 kV and that 
the minimum obtained from the standard method is there- 
fore spurious. 

It is next of interest to ask how the spurious 


intensity minimum obtained from the standard approach 
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changes as the degeneracy is approached. In order to 
answer this question, calculations similar to those 
Shown in jhig-s4t have been carried out for deviation 
parameters decreasing to 6 = 2.0.9 li has been found 
that the depth of the minimum remains essentially the 
Same was tne degeneracy 1S approached. However, the 
range of accelerating voltages over which the minimum 
occurs (i.e., the range of accelerating voltages over 
which the intensities obtained from the standard and 
Ssub-matrix theories are appreciably different) decreased 
as the degeneracy was approached. “For example, for a 
deviation parameter of 6 = 2.00007 the two theories 
depart significantly over a voltage range of approximately 
4 kV, while at 6 = 2.001 the corresponding voltage range 
was 30 kV. 

The question then arises as to whether the 
difference between the standard and sub-matrix theories 
persists at the exact degeneracy. In order to answer 
tines Question, , analyses based “onsboth sanaly Glical Consus 
derations and numerical calculations were carried out. 
The analytical approach was based on a Bloch wave symmetry 
argument and showed that for centrosymmetric crystals, 
q(23)= O and the sub-matrix theory reduced to the standard 
theory, i.e., the spurious minimum disappeared..- The 
numerical calculations showed that the spurious minimum 


persisted even when the deviation parameter was set equal 


als ) | a 


ot tebxo at .Berdoaotaqs 2f yosaatepsb edt 86 niches 

sayd¢ 03 telimte erottsluolseS vnoltesup eid towers 
acijeiveh ro% duo beitiso nasd ovad [6 .piat ne woe 
bot deed ead JL .0.S = & of paidsatseb eretometag 


eit vylleiynsee: aniamet muminim ofa to dtqeb odd tens 


etd ferTrswotll boriquonds af yon 1sriepab ait 26 antee 
ayainim e+ doidw isvopediov poits1sfeco® Io spnsx 
Wo aopstinpv paisgereisesss To spast sn _-3.1) e@Teos0 


Seeds att mata benlétrdo eelttensini and nginw 
bessarpeb (treretith ybdefoestlqqe 975 esixroont xingsm-dge 
6 tol ,Slqnaxs 167 sbsripso1ggs e#6w yos1snspeb edz e6 
esitoarnt4 owl sft TOO00,.4 = 5 Io 3zSstamst6q nots iveb Q 
ylevemixorggs io spits spsilov 2 19v0 vyiinseottiapiea Iasqeb 
gpfax sostiov poibmegesttoo ent 160.S = 4 96 slinw ,vAb 
-VA Of BBW - 
add tenddedw od es aeette nodt nolssoup sn? 
esrroatit xixoew-doge boas Gisbabtie offf nsewied sonsysttib 
+awens 63 webtoa al ,yasTsHopsd s5axo oat 36 aera 
-~janop Lesitylena Mod no boead aoeylsas ,noiszeup aids | 


od beiatss avew enoitsluvisy ispiasmba bas ey 7 


7 a ama ae esan rte bared dew nosorqae Les Last ie 


va - acs an phins - iene pei 


7 
% 


218 


to 2.0, but that the width of the minimum depended on 
the numerical precision used in calculating the eigen- 
values and eigenvectors. For example, the Jacobi 

method of matrix diagonalization (SSP, 1970) used in 

an IBM system 360 Fortran IV single precision routine 
(approximately 7.2 significant Pees ites produced 
a spurious minimum whose width was approximately 0.05 kv 
when the deviation parameter was set exactly equal to 
2.0. The same routine in double precision, however, 
reduced the width of the spurious minimum to lyse e kV. 
These results indicate that the spurious minimum will 
always remain in numerical calculations based on pertur- 
bation theory for the non-degenerate eoee but tor-ehagh 
precision it will become very narrow. It is clear in 
comparing the numerical results with the conclusions 
based on an analytical approach to the problem, that 
even though the deviation parameter can be set equal 

to 2.0 in numerical calculations, the effective value 

of. the, deviation. parameter,is not.exactly. 2.0. 

Pes ye the standard theory should not be used 
in-numerical calculations close to. the, critical voltage 


even when the exact Bragg condition of.the second order 


reflection is satisfied. 
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8.3.1 (b) An examination of the critical voltage intensity 
minimum in copper 

Calculations were also carried out for the (222) 
intensity in copper at small departures from the exact 
second order Bragg condition. The critical voltage minima 
obtained with the standard and sub-matrix theories for a 
crystal thickness of 5000 A are shown in Fig. 42. 
Thirteen reflections of the (111) row were included 
in this calculation and the deviation parameter 6 was 
taken toywbe 2.002) It can be seen from this figure that 
a Critil@al Wwoltage intensity minimum is cbtained from 
both theories but that the minimum is more pronounced 
for the standard theory. The range of voltages over which 
the standard and sub-matrix theories differed significant- 
ly was found to decrease as values closer to 6 = 2.0 were 
considered. . (However, it should be noted that the 
numerical results obtained using the standard and sub- 
Matrix theories again differed for a very small voltage 
range when the deviation parameter was set equal to 2.0. 
The range of voltages over which the differences occurred 
were similar to those discussed for molybdenum in the 
previous section and were found to depend on the numeri- 
cal precision of the matrix diagonalization routines.) 

When examining the critical voltage intensity 
minimum at small deviations from the exact Bragg condi- 


tion of the second order reflection, the question arises 
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Of ‘the range ofa values of 0 over which the crirticar 
voltage intensity minimum is obtained. It has already 
been seen that a minimum is Seana LOL Liven i222) 
Mcenoi tye toecopper, LOL. 2.0sieo 2.002... The variation 
in the intensity with accelerating voltage for the 
standard and sub-matrix theories at 6 = 2.02 is shown 
iNet On 40em ete Can. be seen from itias figure thatea 
‘distinct intensity minimum is still obtained. fFurther- 
more, the standard and sub-matrix theories differ 
significantly over a wider voltage range than for 

6.= 2.002'fin Fig. 42. ,As values ‘of 6 > 2.02 were 
considered, it was.found that the agreement between 

the two theories improves but the intensity minimum 


becomes less pronounced. 


8.3.2 The Standard Theory Mechanism for the Critical 


Voltage Effect 


In the previous section it was found that a 
critical voltage intensity minimum can occur for small 
departures from the exact second order Bragg condition. 
The Bloch wave: mechanism which gives rise to this 


minimum will be examined in this and the following 


section. 
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Pweasetirst important toillustrate thatthe 
Bloch wave mechanism for small departures from the exact 
second order Bragg condition is in fact different from 
the mechanism at exactly 6 = 2.0. This can be done by 
examining the excitation amplitudes of Bloch waves two 
and three obtained by using the standard theory. In 
order to help visualize the change in symmetry of these 
Bloch waves at the critical voltage and for 6 exactly 
equal to 2.07 1h 1s susetul to -piot the values ‘of 


ge2) ci) iis) as a function of accelerating 
fe) 2g fe) ag 


ands 
Voltage. sucht a plot (for the (222) retlectionyin copper 
is shown in Fig. 44(a). The change in symmetry and the 
critical voltage correspond to the vertical line in 

thas b1lguresand- occur at abotita 390 KV.” “The scritical 
voltage mechanism obtained from the standard theory for 
6 exactly equal to 2.0 has already been discussed in 
Seccion S.2.)., (Le-wasenoted in Section 8.2 >that the 


= = one at the 


intensity minimum arises because y 
critical voltage and the excitation amplitudes of these 
Bloch waves are nearly equal in magnitude and are oppo- 
site in sign.) 

In order to compare; the°critical voltage 


mechanism considered in this section with that to be 


considered in the following section, it is convenient 


Ss | | | | i ee 


ait tedt edetteuiit oe sast0qnt dextt te 
 '$o6x9 oft moxt esowwtisqeb Llame tot ie thd evden 
moxt $aenst tb Jost ai et noltibaos pps rd asb10 bnodse 

vd sfob od mao eisty .O.S = 3 diveiels +6 mains room ont 
ows asvew foold to esbostitqns Awad tied ant paininsxs 

nl  .yroeds brebaste sfl+ pniew yd pentetdo sexds bas 
geadt to yrdemmye 7 epneno eft esilsavetv gqled ot teb10 
yigosse * rot bas spstiov Isoitino sft 36 eovew fino La 

to omihew sit tolq ot Ivteev ai Ji ,0.8 OF Isups 
prisstaleoss to aoisgonut 6 286 pet a baa eo 
teqqoo ai nolgoelte1 (SSS) eft 10% Jolq 6 dowe . spssiov 
ed bts yosoumve ai epastio ofT .(6)o*® .pri leeches ai 
ai snail [spitxeav adda oF Brnogesitce eestlov Leoitizo 
[eoitixo ott Va Ot tuods te xw990 bas siwpit ein 

tot yiosHs Srabaste sit mo tt Settasis mainsroem epstiov 
nt beeevoaib need ybsorls esd O.6 ot isups yljosxe é 
eit taft &.8 noidos® at Pai esw tI) .&.8 mottose 

adg +6 Su = (37, seusged ‘goeiss momisim ysienssat 
event 30 esbusiians soiietioxe edt bas pwr Iisoitito 


State nee Pipe, aes ak Loupe yiisen exe aevew doole 


020 


Excitation Amplitude 


(a) 
Oo 
Vp) Ve) 
os) =) 
| 
| 
> 7 é 
Ss NS, | © 
ao 3 2 gO 
ecu | cas Sie 
2 ! : 
E | H 
| 
oO | (=) 
Va) | Ve) 
: a i 
'_9.50 0.00 0.50 @ 
Real 
Cis) Pee = oo Ue 
: : GEE oe ' 
PIpse4 492K): Al plowmobeCeeec5e om =n2y3)haspa Sanetion 


of accelerating voltage at the exact Bragg 
Gondi GLone Ole th Renee 2G¢e(bimeandgi(cyk 
Plots ot cliesi (i) = 2), 3)5 ein Ftheecompilex 
plane for voltagés below and above the critical 
VOLLAge. 


C28 


| pe 
E 
a 
| 
Oo 
“=o ft 
im 
i oH 
i ~: 
L a ee g 
ra] 
1 Veal fa) 
ci we oa -¢ 
(Vax PHI LOV 
(A) 
a 
; as + ' . 2 
. Haas 
| j | 
; : 
7 | < 
; oF 
| : a) 
oe ee ae ee | 
(On 
| 3 
J 
, | | Pent | 
' re 
: M 
a ae f 
i] z : 
02.0 00.9 od. 0- Od. 0 
inet 


VA ohO = VY (b) Vx 08 


pt ry (4) 
> jo (has. -.&)". o to dol ws Hats i. 
bai: mas a 2 veld SPE aids baidszols ies 


By sear iy "BESS Sots 


he , 


ott 


LO pLotetiemvalues of cf) o (2) and oo in the 
. fe) 2G fe) 2g 


complex plane. The values of these excitation ampli- 
tudes for voltages below and above the critical voltage 
are shown in Figs. 44(b) and 44(c). It can be seen 
from Fig. 44(b) that, below the critical voltage, the 
excitation amplitude of Bloch wave three is negative 
while that of Bloch wave two is positive (i.e. the two 
excitation amplitudes are out of phase by 180°). For 
voltages above the critical voltage, the signs of the 
excitation amplitudes are Pore cera This change in 
Sign occurs when Bloch waves two and three interchange 
symmetries at precisely the critical voltage. 

It is next of interest to ask how the values of 


(3) 4 (3) 
2g 


: vary with accelerating voltage 


c'2)c'2) ana e 
OTT 29g 


for a small deviation from the exact Bragg condition. 


The results obtained for 6 = 2.002 are shown in Fig. 45. 
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In Fig. 45(a) the excitation amplitudes have been plotted 


as a function. of accelerating voltage, while in Figs. 
45(b), (c) and (dad), the excitation amplitudes for 300 
kV, 390 kV and 480 kV have again been plotted in the 
complex plane. It can beuseen LHom this figures thar, 
in contrast to the case for 6 exactly equal to 2.0, the 
interchange of excitation amplitudes of Bloch waves two 
and three does not occur at a particular accelerating 


voltage. Rather, this change occurs over a range of 


accelerating voltages equal to about 50 kV at 6 = 2.002. 
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Astthe* voltage is increased from about 360 kV to 
380"kV;¥the* magnitudes of ‘both excitation amplitudes 
decrease. However, the excitation amplitude for Bloch 
wave two changes sign at a lower voltage than that 
for Bloch wave three and these excitation amplitudes 
are of the same Sign near 390 kV. With a further 
increase in the voltage, the excitation amplitude of 
Bloch wave three changes sign and the magnitudes of 
both excitation amplitudes then increase. 

The first problem in analyzing the results shown 
Lieto > Ss sevet vor defining a Critical voltage since 
sucha defination’ in terms of Bloch wave symmetry is 
ho longer possible. Under the circumstances, it is 
reasonable’ to=take™ the critical voltage to be that 
voltage which corresponds to the closest separation of 
branches two and three of the dispersion surface (i.e. 
the voltage at which Bloch waves two and three are 


closest to being degenerate). This corresponds to the 


POiIMuewiternrg., 44a) at Whicll tie curves for aire and 
eS Cross anid 1S Again equaleror about S90 KV he 
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Can *be~seen from. Lgs.—45-(a)- and=45(c)-j—theretore, uit 
at the critical voltage the excitation amplitudes of 
Bloch wavesmtwoeand \threegarenorctheisame Saignse Thas 
is .in,contrast to (the Situation cated: =en0iain Fig. 44 
where,the excitation «ampl meudesmare tof topposite -sign-. 


However, the magnitudes of the excitation amplitudes 
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at the critical voltage are very much smaller for 

6 = 2.002 than for 6 = 2.0 (compare Figs. 45 and 44) 
and as a result a.minimum in the intensity is still 
obtained (see Fig. 42). 

The results just presented illustrate that 
the critical voltage mechanism for small deviations 
from the exact second order Bragg condition is quite 
different from the mechanism at 6 exactly equal to 
2.0. (It is important to emphasize that the same 
qualitative arguments as those discussed above Hola, 
even for an infinitesimally small deviation from the 
exact Bragggcondi tion.) However, 2 was. Shown. an 
Sectrone6.3, 1. that: Bloch wave mixing 1S important and 
should be taken into account for small deviations from 
6 = 2.0. Accordingly, detailed examination of the 
critical voltage mechanism for these diffraction con- 
ditions should be carried out by using the sub-matrix 
theory: The sub-matrix mechanism for the critical 


voltage effect will be considered in the next section. 


8.3.3 The Sub-Matrix Mechanism for the Critical Voltage 
BErecc 


The mechanism for the critical voltage effect 
at) small déviations from the Exact Second order Bragg 


condition can be seen by examining the contributions 
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which the perturbed and unperturbed Bloch waves make 

to the sub-matrix diffracted beam intensity expression 
of equation 6.45. (In this thesis, the Bloch waves 
obtained from the standard theory are referred to as 
unperturbed Bloch waves. On the other hand, Bloch 

waves formed from a linear aac sac ye ter of, the Bloch 
functions of the unperturbed Bloch waves are referred 

to as perturbed Bloch waves. In the sub-matrix approach 
considered here, Bloch waves two and three are perturbed 
Bloch waves while the remainder are unperturbed Bloch 
waves.) In thie analysis, it is necessary to use a 

2x 2 sub-matrix hecause the mixing between Bloch waves 
two and three was found to be important for values of 6 
close to 2.0 (see Section 8.3.1). However, the exact 
theory is not required because good agreement was 
obtained between the sub-matrix and exact theories for 
the medium atomic number materials. 

An analysis of the magnitude of the Bloch wave 
excitation coefficients in the sub-matrix intensity 
expression 6.45 showed that, except for very thin 
crystals, Bloch waves two and three were the only 
primary Bloch waves. (Thus, in this respect, the 
situation for small deviations from 6 = 2.0 is similar 
to that for 6 exactly equal to 2.0, discussed in 
Section 8,2.) The intensity expression can therefore 


be written 
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where A represents the small contribution made by the 
unperturbed Bloch waves n # 2,3. The excitation 


amplitudes X 


faves in this equation are complex and 


are obtained by the methods given Li LoectLon bea. 


8.3. (a)#The sub-matrix mechanism for % =-2.003 


In order to see how the minimum in the second 
order intensity arises, it is necessary to examine 
the excitation amplitudes for Bloch waves two and three 
in equation 8.3. These values are complex and therefore 
the manner in which they vary with voltage will be 
plotted in the complex plane. The values of these 
elements at 6 = 2.003 for the (222) reflection in 
copper are shown in Fig. 46. (This value of 6 has been 
chosen in order to first show a typical example of the 
sub-matrix mechanism for the critical voltage. Other 
values of 6 will be considered in later sections.) 
It can be seen from Fig. 46 that, for voltages well 


below 390 kVj,, they reat, paxt.or x FBerG) 


is negative 

in sign and is considerably greater in magnitude than 
the aera part. With an increase in voltage, 7the 
magnitude of the imaginary part increases relative to 


that Gietne real part until at 390 KV the inaginary 


part is large while the real part is small. With a 
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The excitation amplitudes re ea gsh et OL 
perturbed Bloch waves two and J 

three plotted in the complex plane. 

The numbers between 375 and 405 indicate 
the accelerating voltage in kV. The 
voltage range considered was 360 to 

420 kV and the voltage increment was 

1.5 KV... dhirteen retlections of pene 
(111) systematic row in copper were taken 
into account and,o = 2,003,. 
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further increase in the accelerating voltage, the 


Pealevar tro. sat A! 


again becomes considerably 
larger in magnitude than the imaginary part, but the 
real part is now positive. The excitation amplitude 
of Bloch wave two changes in the opposite manner to 
that of Bloch wave three as the voltage.1s increased ; 
As can be seen from. Fig. 46 the interchange of the 
excitation coefficients occurs over a voltage range 
OL Vabout SOLKV.a Ths. 1s approximately whe woltage 
range over which the interchange of the excitation 
amplitudes of the standard theory occurs, and is also 
the voltage range over which the intensities obtained 
from the standard and sub-matrix theories differ 
Significantly. 

The minimum in intensity which occurs close 
to the critical voltage can now be explained by noting 
firstly that, for jill voltages, the excitation ampii= 
tudes of Bloch waves two and three are nearly equal 
in magnitude and nearly out of phase. The phase 


difference at 390 kV.for example is .155°%. Furthermore, 
yt (4) ee 


for voltages close to the critical voltage, a 
The contributionsmade by these two perturbed Bloch waves 


to the diffracted beam intensity will therefore inter- 


fere destructively and a minimum in intensity results. 
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823.3 (b) The sub-matrix mechanism for 6 < 2.003 


In Section 8.3.3(a), the character of the sub- 
matrix mechanism for the critical voltage has been 
determined by examining this mechanism at a value of 
Gp=a2.005storvthe, (111) systemabticerow an copper. 

It is next of interest to examine how this mechanism 
differs for other values of 6. Calculations have 
therefore been carried out for 2.0 < 6 < 2.003 and 

Fors Oye en00R ws Por 2. 0M 60 PTDOSTi te nespeen found 
that the sub-matrix mechanism for the critical effect 
td Smillar tO -ehat GOL. — 2. 00Sinmthat thie, excitation 
amplitudes are large in magnitude and nearly out of 
phase for all voltages. However, the shapes of the 
curves obtained by plotting the excitation amplitudes 
in the complex plane are found to change. The excitation 
amplitudes of Bloch waves two and three for 6 = 2.002 
are. shownlingbigque4ieas Ttgcan¥Ybe seen from this Pagers 
that, for voltages well below the critical voltage, 


(3) G1 (3) is negative in sign and is 


Che, tealspart, of ex 
considerably greater in magnitude than the imaginary 
part. /eWithe an increase an) voltage, tne real pare 
increases in magnitude while the imaginary part at 
first increases and then decreases. At a voltage of 
approximately 390 kV, the real part then changes sign. 


With a further increase in voltage, the real part 


decreases in magnitude while the imaginary part at Lars & 
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BiG 47% Ths Seton amplitudes X y) 
x'3)05(3) gor 6 = 2.002. The nutbers 
between 375 and 405 indicate the accelerat- 
ing voltage in-kV.° The voltage range 
considered was 375 to 405 kV and the voltage 
increment was 0.25 kV. 
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increases and then decreases. For voltages well above 
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positive and considerably greater in magnitude than 


theycrivicallvoltage,stne real part of xX is 

the imaginary part. The excitation amplitude of Bloch 
wave two changes in the opposite manner as the voltage 
is increased. This interchange of the excitation co- 
efficients occurs over a voltage range of approximately 
30 kV. It can be seen, however, that for all voltages, 
the excitation amplitudes of Bloch waves two and three 
are large in magnitude and nearly out of phase. Since 


yt 7) ey at the critical voltage, a minimum in 


aay’ 
the second order diffracted beam intensity is obtained, 
in a manner similar to that for 6 = 2.003 in Section 
Srysp.3.la.) i. 

The sub-matrix mechanism for the critical 
voltage effect was also examined for 6 = 2.0017, 2.001, 
and 2.0. These results are shown in Fig. 48. It can 
be seen from this figure that, as values of 46 closer 
to 6 = 2.0 are considered, the magnitudes of the 
imaginary components of the excitation coefficients 
are reduced. Furthermore, the change in the magnitudes 
of the real parts is smaller. It is also found that 
the voltage range over which the interchange in the 
excitation amplitudes occurs becomes smaller. For 
example, this voltage range is equal to about 30 kV 


at 6 = 2.002 and 20 kV at 6 = 2.001. It can therefore 
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Fig. 48. The excitation amplitudes for Bloch waves two 
And three “for (a) "6 =" 2001/73 4b) 6 = 2.0015 
and (c) 6 = 2.0. The voltage range considered 
was 375 to 405 kV and the voltage increment 
was 0.5 KV. 
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Descecensthateass the value Of 6 approaches 2,0, the 
sub-matrix mechanism for the critical voltage inten- 
Sity minimum reduces to the standard theory mechanism 
fOm Wonexactiy equal to 2.0) (see Fig. 44). However, it 
should be noted that, for the sub-matrix theory, the 
interchange in the excitation coefficients is still 
found to occur over a very small voltage range even 
whensthe value Of 6 12s Sechexactly ‘equal to 2.0. “This 
Lange US stound to,decredses lt Matrix diagonalization 
routines of higher precision are used. These voltage 
ranges are similar to those for which .the intensities 
obtained from the standard and sub-matrix theories were 
LOUHG sLOSdLiter Signiticantly Or GO set veractly oguald 
to zero (see Section 8.3.1l(a)). Analytically, however, 
Rica Dem Siown stiate toro exaculy (equa la tO.2..0, eee! 
is zero and the sub-matrix mechanism for the critical 
voltage effect is equivalent to the mechanism obtained 
from the Standard, theory, 1.e., the excitation =ampli-— 


tudes remain real and interchange precisely at the 


Cre ECalesvG tage. 


8.3.0 (Cc) The sub-matrix, mechanismelor Oo) ee .00s 


It is also:of interest to examine the sub- 
matrix mechanism for the critical voltage at values 
Ou vec neateretiane 2 UUs, muGalCulatvonssoL the exci Las 


Cicim anplicodes, Lor =——2.02 have therefore been 
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Carried out and the results are shown in Fig. 49. 

It can be seen from this figure that the interchange 

of the excitation amplitudes occurs over a voltage 

range Of. about 200 kV,which is considerably larger 

than the 50 kV obtained fori6 = 2.003 in Fig. 46. 
Furthermore, the imaginary components of the excita- 
tion coefficients are considerably smaller in magnitude. 
It can therefore be seen that, as the deviation of the 
second order reflection from its exact Bragg condition 
is increased, the sub-matrix mechanism for the critical 
voltage reduces to the mechanism obtained for the 
standard theory. (As indicated in Section 8.3.2, a 
critical voltage intensity minimum is obtained for the 
standard theory for diffraction conditions where 

y '2) % y $3) and the excitation amplitudes are small in 


magnitude and are in phase.) As noted in Section 


re | 


8.3.1(b), however, the critical voltage intensity minimum 


becomes less pronounced for values of 6 > 2.02. 


8.4 The Critical Voltage Effect in Materials of 


Ditfemiung Atomic Number 


In this section, the critical voltage effect in 
materials of differing atomic number will be examined. 
The reason for considering the effects of atomic number 
is that the off-diagonal elements of the B matrix are 


larger for materials of higher atomic number and the 
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perturbed Bloch waves two and 2g 
three plotted in the complex plane. 
The numbers between 290 and 490 
indicate the accelerating voltage in 
kV. The voltage range considered was 
240 to 540 kV and the voltage increment 
was 10 kV. Thirteen reflections of the 
(111) systematic row in copper were 
taken into account and 6 = 2.02. 
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effects of Bloch wave mixing are enhanced (see Section 
7.2). In order to assess the effects of atomic number, 
it is important to be able to distinguish between 

these effects and effects which result from the mixing 
between the quasi-degenerate Bloch waves two and three. 
This will be done by carrying out the calculations 

for 6 exactly equal to 2.0 because, for these diffrac- 


tion conditions, the value of ey 


can analytically be 
shown to be zero (see Section 6.4). Any differences 
between the standard and exact theories can then be 
associated with the other off-diagonal elements of the 
Baath il ey COeemi xing sbetweensOther oat Ts 2Of UbLoch 
waves. As noted previously, these off-diagonal elements 
are larger for high atomic number materials (see Section 
7.2), and therefore differences between the standard 

and exact theories can be associated with the effects 

of atomic number. However, it should be noted that 

for accelerating voltages very close to the critical 


fo9) is numerically found to be 


voltage, the value of q 
non-zero for 6 exactly equal to 2.0. If the standard 
and exact theories were compared for this small voltage 
range, it would not be possible to distinguish between 
effects associated with atomic number and effects 
associated with the quasi-degenerate Bloch waves two 


and three. This problem can be overcome by comparing 


the exact theory to the sub-matrix theory, which includes 
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mixing between these two Bloch waves. Differences 
between the results obtained by these two methods 
can then be associated with the effects of atomic 
number. The results obtained from the sub-matrix 
and exact theories will be presented in this section. 
In order to examine the critical voltage 
effect in materials of differing atomic number, the 
critical voltage intensity minimum was obtained for 
materials of varying atomic number including silicon, 
molybdenum and gold. In these eaicur ea one thirteen 
systematic reflections were considered and 6 was set 
exactly equal to 2.0. The second order reflections 
considered were the (220) in molybdenum and the (440) 
in silicon and gold, and the diffracted beam intensity 
was obtained for a thickness of 5000 A. These calcu- 
lations showed that there was good agreement between 
the sub-matrix and exact theories in ine low and 
medium atomic number materials. This indicates that 
the effect of the Bloch wave mixing associated with 
atomic number is not important in these materials. 
The results for the gold calculation are shown in- 
Fig. 50. It can be seen from this figure’ that’ the 
critical voltage intensity minimum is obtained for 
the sub-matrix. theory at about 740 kV. On.tne other 
hand, the minimum is obtained for the exact theory at 
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be seen that~the-effect of Bloch wave mixing can result 
in a significant change in value of the critical vol- 
tage in high atomic number materials. Lally et al. 
(1972) have previously reported that the critical 
voltage obtained for gold by using the standard and 
exact theories differs by 2%. However, these authors 
did not offer any explanation for this effect. The 
‘results presented in this section show that this 23% 
difference is associated with the effect of atomic 
number and cannot be attributed to effects associated 
with the quasi-degenerate Bloch waves. 

It is of interest to examine the Bloch wave 
mechanism which gives rise to this difference between 
the standard and exact theories. An investigation 
Showed that the 2% difference in the value of the 
critical voltages obtained from these theories arose 
because of changes in the real parts of the eigenvalues 
of Bloch waves two and three. It was also found that 
as the size of the sub-matrix was increased, better 
agreement with the exact theory was obtained. For 
example, the critical voltage minimum for a 6x6 sub- 
Matrix including mixing between the first six Bloch 
waves differed by only 0.07% from the minimum obtained 


for the exact theory. 
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8.5 Summary and Conclusions 


In this chapter the fmportance of Bloch wave 
mixing in the critical voltage effect has been examined. 
It was noted that the exact Bragg condition of the 
second order reflection cannot in general be obtained 
experimentally because of factors such as beam diver- 
gence, specimen bending, and the lack of precision 
which is inherent in determining the crystal orienta- 
tion. The critical voltage effect was therefore 
examined for small deviations from the exact Bragg 
condition of the second order reflection, and it was 
found that it is necessary to take the mixing between 
the Pies ioqorsner. te Bloch waves two and three into 
account. It was also found that it is necessary to take 
this mixing into account in numerical calculations even 
when 6 is set exactly equal to 2.0. This is because 
the effective value of 6 cannot be set equal to 2.0 with 
greater precision than that inherent in the matrix 
diagonalization routines. It was also found that the 
critical voltage mechanism presented in the literature 
does not apply at small deviations from 6 = 2.0. An 
alternate mechanism based on the sub-matrix approach 
was therefore presented. 

The effect of atomic number on the minimum in 
intensity which occurs at the Critical voltage was’ also 


examined. It was found that Bloch wave mixing can lead 
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to a Significant change in the voltage at which this 
minimum occurs when high atomic number materials 
are considered. 

It can be concluded from these results that 
the effects of Bloch wave mixing should be taken into 
account when calculations of the critical voltage 
effect are carried out. In addition, however, it 
should be noted that these effects are important over 
a range of voltages on either side of the critical 
voltage and over a range of values of 6 on either side 
OL 0 =. 2.0 (See also Section 7.3). ,The results or 
this chapter are therefore also of importance in cal- 
culations of image contrast. For example, in copper 
the ¢ryiticad. voltage: for thebi(222) reflection 1S (390 (KV, 
yet significant differences between the standard and 


exact theories were found at 300 kV for 6 = 2.1. 
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CHAPTER’ 9 


SUMMARY AND SUGGESTIONS FOR FURTHER WORK 


o.1 Summary 


The purpose of this thesis has been to assess 
the validity of some of the approximations in the 
dynamical theory of electron diffraction. The s-value 
approximation for calculating the diagonal elements of 
the A matrix was considered it Ghapters,4netitdwas 
found for the diffraction conditions considered that 
this is a reasonable approximation. This was followed 
in Chapter 5 by an assessment of the validity of the 
kinematical theory for large values of 6. These results 
showed that this theory is not as good an approximation 
under these conditions as it had previously been consi- 
dered to be. The approximations which can be used to 
take absorption into account in the dynamical theory 
were then considered in Chapter 6. This led to the 
development of a general formulation of the dynamical 
theory of e@lectron cdifiraction ancluding absorption. 
Itewas, found that;thepintroduction of the imaginary 
potential leads to a mixing of the unperturbed Bloch 
waves characteristic of the system without absorption. 
The assumptions required to reduce the general theory 


to the low order perturbation theories previously 
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published in the literature were then examined, and a 
hew approach referred to as the. sub-matrix approach 

was developed. Numerical calculations employing this 
theory were presented in Chapter 7 which showed that 
itnis-important.te take the mixing of Bloch waves into 
account in high atomic number materials or when quasi- 
degeneracy or degeneracy of the kinetic energies of 

the unperturbed Bloch waves is obtained. A comparison 
of approximate and exact methods for taking absorption 
into account showed that the sub-matrix method has a 
number of advantages over other approximate methods. 
Firstly, it is applicable eee conditions which lead 
to either non-degenerate, quasi-degenerate or degenerate 
Bloch waves and secondly, this method can give very 
close agreement with the exact theses while at the same 
time providing advantages from a computational point of 
view. It was then shown in Chapter 8 that it is 
necessary to take Bloch wave mixing into account in 
order to obtain a physical understanding of the con- 


trast mechanisms involved in the critical voltage effect. 


9.2 Suggestions for Further Work 


There are a number of areas in which the work 
presented in this thesis can be ‘extended. “The first of 
these follows from the assessment of the kinematical 


theory “Carried out “In Chapter 5.° “The results of Section 
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5.4 showed that non-systematic reflections can have 
an important effect on the weak-beam image of reflec- 
tion g not only when S, * 0 but also when Sy * Ba 
However, it has been assumed in the literature that 

a non-systematic reflection will only be important if 
lies CLOSeetOomlts pragg cOndte1on (a,c. > Sy, * Oj e 

It would therefore be of interest to carry out an 
experimental investigation of the effect of non- 
systematic reflections on weak-beam image contrast 

for diffraction conditions where S, > ae Furthermore, 
it should be noted that in carrying out theoretical 
calculations of weak-beam contrast, it is necessary 

to take all reflections into account which have Sy 
values.close to zero or Bes SincCemtiic gaCk. is OG 
reflected in Cann's strong beam criterion of equation 
3.2, an improved criterion needs to be developed. 

A second possible extension of the work pre- 
sented in this thesis is concerned with the importance 
of Bloch wave mixing in determining image contrast 
obtained under various diffraction conditions. It was 
shown in Chapters 7 and 8 that Bloch wave mixing can 
be cer important under quasi-degenerate diffraction 
conditions. Since these conditions can occur when 
non-systematic reflections are taken into account 
(Cann and Sheinin, 1974A, 1974B), it would therefore 


be of interest to examine the effects of these reflec- 


tions on Bloch wave mixing. 
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A third extension of this research would be to 
examine the importance of Bloch wave mixing in non- 
centrosymmetric crystals. For these materials, the 
PrCperties joie thesott—-divagonal (BP elements are different 
(for example, the B matrix will not be symmetric). 
Furthermore, it would be of interest to examine the 
critical voltage effect in these materials since he 


(233 ) 


value of q will not in general be zero for 0 = 2.0. 


This is in contrast to the case for centrosymmetric 


(23) 


crystals where “qd = 0 fer. 6 = 2.0 (Sée Section) 6.4) . 
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APPENDIX A 


Approximate Methods for Calculating the Eigenvalues 


and Normalized Eigenvectors of a Symmetric Matrix 


In this appendix, the first and second order 
approximate methods for calculating the eigenvalues 
and normalized eigenvectors of a symmetrical matrix 
will be considered. The approximations used are 
Similar to those outlined in texts such as Margenau 
and Murphy (1956). 

As a first order approximation to be used in 
diagonalizing, the symmetric matrix A, Jet it be 
AScune@d anata dil eLie Ori-dvagoial.elemencs of Avare 
small and can be neglected, i.e. that they can be taken 
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The corresponding eigenvectors can be found by substi- 
tuting these first order eigenvalues into the first 


order eigenvalue equation 
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O Ay 5 Y O eeee C5 = OO e (A. 3) 
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For example, when DM man iS esubsit tuted Minto stnls 
equation, it can be seen that the elements of the 
eigenvector c(t) are proportional to 
(1) 
Si 
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Since wie eLement Cie must be chosen so that the 
eigenvector ct) is normalized, the GCigenvector Co 
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In this thesis, when the same approximations 
are used to calculate both the eigenvalues and the 
eigenvectors, the eigenvectors will be considered 
correct to one lower order than the eigenvalues. MThis 
terminology has been adopted in order to be consistent 
with the terminology used in perturbation theory. In 
perturbation theory, the same approximations which yield 
the energy to a given order also yield the wave functions 
to the next lower order. The eigenvectors of equation 
A.5 will therefore be considered ome to zero order. 

In order to calculate the eigenvalues correct 
to, second’ order’ and* the eigenvectors COPLECE, COMI not 
order, it *is*necessary’ to develop a second order appro- 
ximatwdon tothe A. matrix.» ° In" the: second orderi approexi- 
Mat#on, théVvVeéitects-of- thé off-diagonal elements or A 
are taken into account in an approximate manner by 
assuming ea particular element aS (AF 4) Was “only 
have an effect on the eigenvalues and eigenvectors i and 
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ane taken intdgaccount lin calcullating ithe nee eigenvalue 


and eigenvector are A; and a Wiere tj = "L,N and 44°. 
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As an example of the use of this second order approxi- 
mation, approximate expressions for “and and gut) will 


first be calculated. These expressions will then be 
generalized. 
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The first order eigenvector eee can be obtained by 


Solving. all’ bucethe first equacion in A.7 to give 
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The elements of the eigenvector C are therefore pro- 


portional to 
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is non-degenerate, the denominator in equation A.8 


will be large. If the off-diagonal elements of the 
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a good approximation, therefore, the eigenvector © 
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elgenvector C , correct to first order, then becomes 
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where R is a column vector containing the elements 
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p(n) in the 4° row. ‘The remaining eigenvectors can be 


calculated in a similar manner and the first order matrix 


Ofselgenvyvectors Cecah be written 
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in the nti column. 


An wimportant property of thevmatrix of first 
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to make this product equal to the unit matrix.) 


R and examining the approximations required 
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The second order expression for the eigenvalue 
yi can be found by substituting the expression for 
C in equations AUT l-and ays intox=the-first—equation 
of the set of equations A.7. The eigenvalue oe 


correct to second order, then becomes 
(A 4 3,) 


When the other eigenvalues are calculated in a similar 
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APPENDIX B 


Typical Values for the Elements in the €C Matrix of 


Table 2 and for the Bloch Wave Excitation Amplitudes 
(hy) ¢h) 
Co Cg 


In this appendix it will be shown that the 
=s 
decrease in magnitude of the terms U, as |h| increases 
is the controlling factor in determining the relative 


magnitudes of the Bloch wave excitation amplitudes 
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obtained from the C matrix|ofeTablie 2.: This will be 
dohe by examintrig the’? matrix and the excitation 
amplitudes under diffraction conditions typical of 
those obtained when a high order systematic reflection 
iS not close; to ats Bragg condition. Typical values 
fOr che firs. yewo Hows Of EheaC matrix Of Table, 2 

are presented in:Table B-l. In obtaining this matrix, 
the sixureflections (O00) to (660) of the (110) 
systematic row in molybdenum were taken into account. 
The. accelerating voltage was 100 kV and 6 = 5.5. It 
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|h| is in fact the controlling factor in determining 


the relative magnitudes of the excitation amplitudes 
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APPENDIX C 


Typical Values for the Second Order Correction Terms 


In vEquUatcTON) on 


In this appendix it will be shown that the 
decrease in magnitude of the terms U, as |h| increases 
is the controlling factor in determining the magnitudes 
of the second order corrections for the eigenvalues 


yt O7 | ye) 


Of equativonws.2. Ihis will be done by 
examining eee correction terms under diffraction 
conditions which are typical of those obtained when 
there are no high order systematic reflections close 
to wher Braggycondicion. | Typical values for the second 


(0) and y 9) of 


order corrections for the eigenvalues y 
equation 5.2 are given in Table C-1. In obtaining 
this apie, the six reflections (000) to (660) of the 
(110) systematic row in molybdenum were taken into 
account. The accelerating voltage was 100 kV oe 

6 = 5.5. @tecan be seen from Table C=-1 that the 
numerators of the correction terms tend to decrease 

in magnitude as | 4 | increases and that this is in fact 


the controlling factor in determining the magnitude of 


the correction terms. 
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